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Synchronous, permanent magnet, and induction maohines 
are modeled using computer programs. The computer programs 
^incorporate an optimization algorithm which converges on 
lowest weight, volume, and inefficiency. Machine designs 
for high and low rpms are periormed, with a varying number 
of pole-pairs. The machine designs are analyzed to find the 
optimum combination of generator and motor for inclusion in 
a naval ship propulsion system. ' S: "J „ r . 

Three ships are used for the systemsstudy a a baseline 
mechanical transmission ship, a ship retaining the same sub¬ 
division as tho baseline but with the electric machinery, 
and an electric transmission ship with subdivision and 
machinery box arrangement chosen to benefit from the in¬ 
herent arrangeability of electric transmissions. ' r /.■ 

Two generator/motor combinations <•*«>used in the final 
ship analysis/ Both employ a 3600 rpm, six-pole synchronous 
generator, which turns at the shaft speed of the prime 
mover. One combination uses a 180 rpm, direct-drive, 16- 
pole synchronous motor, and the other useeo-an 1800 rpm, 
geared, 8-pole synchronous motor. Power converters are used 
in both combinations to control motor speed. 

The geared combination in the rearranged ship 
demonstrated the best endurance speed efficiency, reducing 
the endurance fuel load by 18%, while maintaining the maxi¬ 
mum and sustained speed of the baseline ship. The savings 
in ship volume translated to an additional twenty Tomahawk 
missile cells in the rearranged ship. When the fuel load 
was held at the tonnage of the baseline ship, endurance 
range inoreased as muoh as 25%. 

Permanent magnet machines were not competitive in this 
study due to their high weight and volumer evan though their 
individual machine efficiency was the highest of all types. 
Induction machines were not used as propulsion generators 
because of the inherent difficulties in control. The induc¬ 
tion machine motor candidates were not competitive because 
of off-design-point inefficiency. 
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Chapter One: Introduction 


The use of eleotrio drive as a propulsion method for 
naval ships brings to the ship design process improved ar- 
rangeability and efficiency, though electric machines may 
increase the weight of the plant. Water-cooled electric 
machines are being studied today for ship transmissions 1 ; 
these are smaller and lighter, for the same power output, 
than air-cooled machines. They promise reduced overall 
weight for the ship through more economic prime mover load¬ 
ing, as less fuel will be needed on board. A review of the 
literature has found no work comparing various types of con¬ 
ventional motors and the effect of each type on the overall 
ship system when used as a propulsion method. St. John [1] 
showed the effects of a superconducting generator/motor 
transmission on the design of a DD863 destroyer hull. Many 
simulations of electric motors and their transients have 
been done. There have been several papers written on naval 
ship integrated electric propulsion systems^. Also, muoh 
effort has been expended in the area of electric motor 
design and optimization. Herein, various kinds of conven¬ 
tional electric machines are modeled. Those machines were 
used in ship designs to find the sensitivity of the designs 
to their use. 


1. Qreene, Mole, Welch, and Seng, "Analysis of a High-Power 
Water-Cooled Sleotric Propulsion System," 8NAMK Trans., 

Vol 86, 1878, pp 140-162. 

2. Ames, "Marine AC Generation Systems," L8K Journal, 

Vol 12(1), pp 13-29. 
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1.1. Review of electric drive 


Navy ships operate in almost every salt water location 
in the world, including the Black Sea and the Indian Ocean. 
Regular deployments are made to the Mediterranean Sea, North 
and South Atlantic Ocean, and all areas of the Pacific. 

Ships transit the Suez and Panama Canals, the Saint Lawrence 
Seaway, and operate in the Arctic and Antarctic Oceans. 

These operations are made under greatly varying environmen¬ 
tal conditions, ranging from the sub-freezing temperatures 
of the high latitudes to the hot, dusty conditions of the 
Middle East. All this variety requires ships (and men) 
capable of sustained and efficient operation under any known 
condition. To that end, naval ships are tremendously com¬ 
plicated systems. 

It is impractical to equip naval ships for every con¬ 
tingency, but ship designers try to include as much as pos¬ 
sible when deciding what systems, equipment, and spares to 
put aboard a ship. Once a ship is designed, or as part of 
the design process, political acceptance of the product and 
its purchase is required. Since ships cost tax dollars, 
there are usually limitations on the size and complexity of 
the design. Still, the designers try to work within the 
given constraints and produce an acceptable and survivable 
(in both the battle and political senses) design. Usually 
this results in a ship that has very small margins in avail¬ 
able weight and volume.* 

Height (displacement) is a semi-direct measure of ship 
cost. Volume is required to place desired systems aboard a 
ship. Therefore, any design change that results in less 
required weight or volume with no decrease in ship 


3. To have small margins is to be limited in the quantity of 
additional system weight and volume that can be added to the 
ship over its lifetime. If a ship is limited in this 
fashion, the flexibility one has in backfitting new systems 
as the ship ages is significantly decreased. 
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effectiveness, is usually a welcome one. 

There are several different kinds of ship propulsion 
systems now in use. They include steam, nuclear, diesel, 
and gas turbine, and there are two principle drive systems: 
mechanical and electric drive. Below is a crude comparison 
of the various propulsion systems, for the purpose of plac¬ 
ing the thesis in perspective. 

Steam plants burning coal or oil have been in use for 
over a hundred years. They require large amounts of prime 
ship volume, in the center part of the ship.* Steam is 
produced in boilers and used to power turbines that rotata 
the shafts and propellers through reduction gears. Steam 
plants are very reliable mechanically, but are not terribly 
economic. The large size of the system components demands 
that the boilers and turbine machinery be placed in the cen¬ 
ter of the ship. This necessitates long runs of shafting 
from the center to the stern of a ship. Shafts typically 
are 18"-24" hollow steel cylinders of two to four inch 
thickness; they are heavy and their required placement and 
length makes valuable volume unavailable for other use. All 
propulsion plants except for electric drive have this ar¬ 
rangement and shafting disadvantage. Steam plants are used 
in all sizes of ships, from the 3000-ton displacement 
frigates to the 50,000-ton battleships to the 80,000-ton 
aircraft carriers. See Table 1 for a list of ship types and 
principle dimensions. 

Nuclear plants are steam plants with a different heat 
source. They also produce steam to power turbines and suf¬ 
fer the same volume disadvantages as conventional steam 
plants. They are also very heavy and very costly. Manning 
requirements are more stringent, since personnel levels are 


4. The center part of the ship is the most useful ship 
region to place and arrange systems. Ship designers try to 
keep free as much center ship volume as possible. This al¬ 
lows much greater flexibility in arranging systems that have 
large objects, such as boilers and turbines. 
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Table 1. Typical principle dimensions for various ship types 


Ship types Total 

shaft horsepower 

Length 

Tonnage 

Patrol Hydrofoil 

18,000 shp 

145 ft. 

240 

Frigate 

37,417 

441 

3880 

Destroyer 

78,555 

510 

6395 

Cruiser (non-nuclear) 

82,462 

556 

8872 

Cruiser (nuclear) 

77,460 

590 

3487 

Carrier (non-nuolear) 

280,000 

1050 

60,000 

Carrier (nuclear) 

270,000 

1100 

73,000 

Battleship 

212,000 

887 

58,000 

Submarine (nuclear, attack) 

15,000 

292 

4640 

Submarine (nuclear, strategic missiles) 

16,000 

410 

7880 


Source: Jane’s Fighting Shies 1986-86 . Jane’s Publishing 
Company, Ltd., London, edited by John Moore. 

These numbers represent the geometric mean of several 
classes of ships within a type and should not be taken to be 
those of a particular ship. Their value lies in the ap¬ 
preciation of the differences between various ship types. 

An apt weight comparison would be that a typioal forty foot 
sailboat might displace fifteen tons. 

The U. S. Navy has other ship types besides those listed 
above. They include amphibious warfare ships, submarine and 
destroyer tenders, and fleet oilers, and supply ships. 
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rigidly controlled and the operators of the propulsion plant 
must be nuclear trained. Nuclear plants are used on 
cruisers, submarines, and aircraft carriers. 

Diesel propulsion plants have high weight to volume 
ratios, making their use costly in weight dollars, but 
cheaper in volume dollars.* They are very noisy, ruling out 
their use in antisubmarine warfare ships. Although very 
reliable mechanically, they require much tinkering and 
tuning. Their specific fuel consumption* is among the 
lowest of all the plants. Medium speed diesel engines are 
not commonly manufactured in the 25,000 hp range, which 
means diesel plants cannot be used in ships requiring high 
shaft horsepower. They are typically used in smaller ships 
as cruise engines and in amphibious warfare ships as main 
propulsion. (Amphibious warfare ships typically have lower 
top speeds than frigates or destroyers.) 

Gas turbine propulsion plants seem to have many good 
points. Thsjy are reliable, quiet, relatively low weight, 
and come in power ranges that are useful in ships ranging 
from 300 ton hydrofoils to 8000 ton destroyers. They 
require large amounts of volume for intake and exhaust duct¬ 
ing, but this is not a great disadvantage. Their fuel 
economy is not as good as other plants, but this is not 
intrinsic to the gas turbine engine. It is a fault of the 
operating method; gas turbine plants have mostly been built 
with mechanical transmissions. Usually one or two engines 
are coupled to each shaft. Jf the ship is proceeding at. 
high speed, the gas turbines are operating at their full 


5. When designing a ship, total ship cost is monitored by 
the use of marginal oost factors. Every additional cubic 
foot or ton of weight added has a marginal cost associated 
with it. When a ship’s total cost is constrained, design 
changes that add cost are disoouraged or must be offset by 
the reduction of other systems' weight or volume. 

6. Specific fuel consumption is the ratio of pounds of fuel 

burned per horsepower-hour. SFC=lbf/hp-hr. 
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load design point and are relatively fuel oconomia. Qood 
fuel economy is not usually realised, however. In the main, 
the ship proceeds at a cruise speed, using one gas turbine 
for eaoh shaft, and the gas turbines operate at about half 
power. Specifia fuel consumption rises rapidly as gas tur¬ 
bine power level drops, which makes for inefficient opera¬ 
tion. 

Usually, the above propulsion plants have ueohanical 
transmissions. This means the main engines, whether they 
are diesels, gas turbines or steam turbines, are mechani¬ 
cally connected to the shafts and propellers. There is 
usually a reduction gear between the engine and shafting. 
These gears are large, very heavy, and expensive. To 
provide an idea of size, the largest, or "bull" gear in a 
typical locked-train double-reduction gear is about seven 
feet in diameter. The reduction gear must be plaued in-line 
with the shafting, thereby using more of that prime ship 
volume. Some mechanical transmissions have cross- 
connections between shafts, but this is not common. 

Electrical transmissions are characterized by prime 
movers of any type providing power to generators. The out¬ 
put electricity ia conditioned and sent to propulsion motors 
via a distribution network. Cross-connection ia done with 
switches and breakers. There can be a mechanical reduction 
gear if it is desired to operate the propulsion motors at 
higher than propeller rotational speeds. Direot-drive 
motors may also provide the desired propeller rpm, e.g., by 
controlling the field current in a DC motor. The propulsion 
motors can be very near the propeller, i.n. aft, eliminating 
the long runs of shafting associated with a mechanical 
transmission. 

Naval ship propellers are of two types, controllable or 
fixed pitch. The pitch of a propeller is the distance the 
ship moves forward in the water for one turn of the propel¬ 
ler. A fixed pitch propeller has this characteristic dis¬ 
tance the same at all times. A oontrollable pitch propellsr 
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can vary this distance by changing the angle of attaok (the 
angle at which the blade slices through the water) of the 
propeller blades, including reversing the blade so that the 
ship moves astern. Controllable pitch propellers are prae- 
tioally required for propulsion plants that have non¬ 
reversing main engines, suoh as diesel and gas turbine 
plants .i Steam plants can reverse their propellers and 
shafts by use of an astern turbine, albeit muoh more slowly 
than a ship with a controllable pitch propeller system. 

Fixed pitch propellers have a slight advantage in efficiency 
(1-3%) over the controllable ones. This is mostly due to 
the large propeller hub required for varying the blade angle 
of a controllable pitch propeller. Quiok reversal of shaft 
direction or propeller pitch means quick ship braking and/or 
ship reversal. This ship quickness is mandatory fox' an¬ 
tisubmarine operations and safe navigation. For example, 
the ability tc stop "on a dime" may be important in a 
crowded sea lane, where a small wooden sailboat has the 
right of way over a powered naval vessel. 

Electric drive seems to combine the best of all the 
propulsion plants. It has all the advantages of a conven¬ 
tional prime mover plus the advantage of electrical cross 
connection and better arrangements. In the cruise scenario 
above, the electric drive ship could have both shafts 
operating from one gas turbine engine. That engine would be 
coupled to an electric generator which would produce enough 
power to run the motors that turn each shaft. The drive 
motors would be placed at the stern of the ship, near the 
propellers, on the same level. The long runs of shafting 
would be replaced by electric cable, which is smaller, 
weighs less, and can be placed in non-prime real estate. 


7. Ships with non-reversing prime movers can also have a re¬ 
versible reduction gear with a fixed pitch propeller instead 
of a controllable reversible pitch propeller. This is new 
technology for the United States and only the latest naval 
ship design, the DDG51, has a reversible reduction gear. 
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Cable is also in many cases cheaper than shafting, espe¬ 
cially to repair. Since one gas turbine would provide power 
to both shafts, it would operate at a higher power level and 
would be therefore more fuel economic. A typical propulsion 
plant might consist of three gas turbines with three gener¬ 
ators. Most ships have an even number of prime movers be¬ 
cause mechanical shaft cross connection between shafts is 
not often used and each shaft in a mechanical transmission 
ship requires the same number of prime movers to balance 
loading at high power levels. The extra prime mover 
requires a lot of weight and volume. An advantage of 
electric drive is that it becomes possible and perhaps 
desirable to use an odd number of prime movers. Each of the 
two shafts would have one propulsion motor. The heavy 
reduction gears could be replaced by the motors, which would 
have an Infinitely variable reduction ratio. The control¬ 
lable pitch propeller system so far required by this gas 
turbine ship would be replaced by the cheaper, slightly 
smaller, and far leas complicated fixed pitch propeller, ad¬ 
ding a small efficiency gain. The hydraulic system used to 
vary blade angle would be eliminated. Ship braking and 
reversal would be accomplished by electrically controlling 
the motor rotation direction, combined with energy dissipa¬ 
tion through the use of resistor banks. 

A disadvantage of this arrangement would be the high 
weight of the propulsion motors. They would be special 
designs and have a high capital cost. Hopefully, the high 
weight of the motors would be offset by the reduction in 
shafting and fuel weight and the possible elimination of the 
reduction gears. The high cost would be made palatable by 
the savings in fuel over the life of the ship. A Life Cycle 
Cost comparison of various propulsion plants, including 
electric drive, is available In reference two. 

The change to electric drive would likely be accom¬ 
panied by an overall decrease in propulsion plant weight and 
volume. The ship could be smaller and lighter, and would 
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require less onboard fuel for the same endurance range (the 
distance the ship can travel without refueling). Since less 
volume would be required for the fuel, the ship could be 
smaller and lighter. Since the ship would then be smaller 
and lighter, less horsepower would be required to aohieve 
the same top speed. Since less horsepower would be 
required, the ship could be smaller and lighter. Thia is an 
example of the design spiral that would result in a smaller, 
lighter, oheaper, more risk-free ship. An example of thia 
type of ship improvement is givon in reference one. There 
is a limit on ship improvement, usually due to the non¬ 
propulsion systems or payload. One cannot make an ocean 
crossing missile ship the siee of a small yacht. 

So why are not all Navy ships electric drive? They are 
not largely because the technology has not existed in a 
usable, fully developed, and manageable form. Because of 
the high cost of naval ships (a small one may cost $350 
million) and the lack of experience with current electric 
drive technology, the Defense Department is reluctant to 
build large electric drive ships. There have been electric 
drive ships, including five battleships with 21 MW shaft 
output and two aircraft carriers with 135 MW shaft output. 
Over 160 escort vessels were built during the Second World 
War with turboelectric or diesel-electric drives ranging 
from about 4.5 to 9 MW. 8 A new class of ocean surveillance 
ship, the T-AQOS 19, is being built with diesel-electric 
propulsion, but it is only a 3500 shaft horsepower (shp) 
ship. 

Electric drive was replaced by conventional mechanical 
transmission plants after World War Two because of the com¬ 
petition afforded by improved gear cutting methods. Double¬ 
reduction locked-train gear transmissions became the 
standard. Since the electric drive ships all had non- 


8. Doyle, T. J. and Harrison, J. H., "Navy Superconductive 
Machinery Program," Trans. SNAME, 1978, p. 20-1. 
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superconducting, air-cooled motors and generators, they had 
higher weight and increased space requirements and suffered 
in comparison with the mechanical transmission ships.* The 
importance of the improvement in power electronics must be 
mentioned. World War Two ships did not have the advantages 
afforded by those electronics. 

Integrated electric drive propulsion must also be men¬ 
tioned. This propulsion plant is the same as any of those 
discussed above, except that ship service electrical power 
is derived from the main propulsion plane, usually by taking 
power off the reduction gears or main propulsion generator. 
Power conditioning equipment, suoh as a cycloconverter, is 
needed to "clean up" the power and change it to fixed 
frequency for use in other equipment. Variable speed con¬ 
stant frequency equipment and concepts embody the integrated 
electric drive concept. The U. S. Navy has investigated 
this in some detail. 

Some requirements of electric drive may be viewed as 
disadvantages. The power from the electric generators has 
to be conditioned to provide frequency control of the 
propulsion motors. The power conditioners add weight and 
volume to the overall system, as well as reducing the 
efficiency of the transmission. Braking resistors, used to 
dynamically and quickly slow the propulsion motors, add more 
weight and volume to the system. There also may be a high- 
frequency radiated noise signature associated with alternat¬ 
ing current systems that may be deleterious to the mission 
of the ship. 

The research done to date has not explored specific 
motor types in detail. How can it be decided whether to put 
a synchronous, inductive, permanent magnet or other motor in 


9. Ibid. 

10. Robey, Stevens, and Page, "Application of Variable Speed 
Constant Frequency Generators to Propulsion-Derived Ship 
Service,“ Naval Engineers Journal, May 1985. 
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the electric drive system? What makes the ship system 
"best"? The effect of each motor type on the ship system 
has not been analyzed. The electric transmissions used in 
the current version of the Advanced Surface 8hip Evaluation 
Tool (ASSET), a ship design computer program written for the 
United States Navy by Boeing Computer Systems, Inc., are 
generic combinations of AC and DC motors and generators, 
using rough estimates of weight, volume, and efficiency. 

Ship designs more involved than the feasibility level need 
detail on just those items. 

Motor design is a well known subject and there are many 
texts on the subject. The use of motors and pertinent tech¬ 
nologies in a ship as a propulsion method is disoussed in 
Qreene, Powell, and Gripp [3]. The advantages of electric 
drive include flexibility of arrangement, controllability, 
variable reduction ratio, reliability, and provision of ship 
service power from the main bus. Jolliff and Greene [4] go 
on to propose a specific water-cooled Advanced Integrated 
Electric Propulsion Plant (AIEPP) for a frigate/destroyer¬ 
sized ship. They discuss the essential characteristics of 
such a plant, establish the feasibility of the drive system 
and identify the method to technically demonstrate the sys¬ 
tem. Acker, Greene, and Jolliff [5] present several model¬ 
ing techniques and scaling relationships that allow estima¬ 
tion of volumes and weights of propulsion motors and gener¬ 
ators, solid state power conditioners, electrioal 
switchgear, and associated electrical propulsion systems 
components as functions of propeller shaft power. A case 
study of AIEPP is given in the paper by Kastner, Davidson, 
and Hills [6]. 

Simulation of electric motors and associated systems is 
a popular topic. Many persons have done work in this area, 
from the micro-consideration of high frequency inductance 
changes to the more macro-consideration of hunting tran¬ 
sients, etc. Smith, Stronach, and Tsao [7] model a complete 
electromechanical marine drive system while Smith, Stronach, 
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Tsao, and Goodman [8] concern themselves more with a marine 
power system, including pump drives. Nonlinearities and 
operational transients are addressed. 

1.2. Optimization 

Optimization is the process of making a system, sub¬ 
system, or idea the best it oan be. “Best" is defined by an 
"objective function," a measure of what is optimum. For ex¬ 
ample, an optimal manufacturing process may produce the max¬ 
imum number of units at the lowest cost. The objective 
function would oombine units-produced with oost in an equa¬ 
tion that could be analyzed to find the proper production 
level. The output of the objective function is a single 
scalar measure of "goodness." It may be difficult to repre¬ 
sent complicated processes with only one number. 

Optimization can be performed on a global or subor¬ 
dinate basis. The optimum motor might be the one that has 
the highest efficiency, even if that efficiency was achieved 
by designing a very large and heavy motor. The sub-system 
(the motor) has efficiency as its objective function. The 
ship in which the motor is to be placed may be optimum when 
its overall weight and volume are the lowest (ignoring oost, 
for example}. A large, heavy motor, then, may not be op¬ 
timum for the ship, even if it is very efficient. A good 
case study of motor optimization is the EPR1 report authored 
by Fuchs, et al. [9] 

"Optimization" can be an ill-defined term but there are 
fairly well defined methods of achieving it. Linear 
programming, Markov modeling, and Monte Carlo schemes are 
examples of these methods. The accessibility of high-speed 
digital computers has made multiple random excursions in a 
multi-variable space a muoh easier way of finding the 
"optimum" solution, provided an objective function and con¬ 
straints can be devised to describe the problem. This 
method of random excursions (Monte Carlo scheme) and ex- 
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amples of It are among the methoda described in references 
ten, eleven, and twelve. 

Monte Carlo schemes take their name from the action of 
the roulette wheel in the gambling casinos of Monte Carlo. 
Around and around the wheel goes, stopping on random num¬ 
bers. If the computing power is available, this is an ac¬ 
ceptable method of exploring a large variable space. It can 
be much quicker than looking at every possible permutation 
of all variables. 

The steepest-descent scheme is so named because op¬ 
timization moves down the sharpest gradient of the objective 
function. From a valid design point, random steps are taken 
in every variable and the design point is moved "downhill" 
toward the objective function over the steepest slope. This 
is different from the “drunkard’s walk," where the random 
steps Are only evaluated on whether or not the objective 
function’s output has improved, not if it is improving at 
the fastest possible rate. 

Optimization is almost always subject to constraints. 

In the previous manufacturing process, warehouse space may 
be limited, so only a certain number of units may be stored. 
This could act to limit production. For motor design, con¬ 
straints include maximum rotor tip speed, maximum current 
density, minimum power output, etc. All constraints should 
be combined with the objective function to yield a 
"constrained optimization." 

1.3. The objective function 

Optimization is not possible without an objective func¬ 
tion. It may be very difficult to devise a good objective 
function for a complicated system such as a ship. It may be 
even harder to find one for a sub-system, of that complicated 
system. There are very many characteristics that could be 
optimizing variables, and assembling them into one objective 
function with all the constraints is not easy. 
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Even with a properly defined objective function, it may 
be difficult to choose among designs that result from the 
constrained optimisation. For example, a low-volume and 
low-weight motor may have poor efficiency. A very efficient 
motor may also be large and heavy. If all three elements 
are important, which is the best motor? Deciding between 
competing designs has been the subject of various papers, 
one of which is by Schweppe and Merrill [13]. In that 
paper, the authors suggest the use of "knee curves,” saying 
that the essential characteristics of a multiple attribute 
tradeoff can be plotted on a series of x-y graphs. Uncom¬ 
petitive designs are easily discerned and discarded. The 
decision process can be limited to only those designs that 
are competitive. 

Table 2. Optimizing characteristics for ship and transmis¬ 
sion 


Weight 

Volume 

Efficiency 

Cost 

Reliability 

Maintainability 

Commonality 

Manning 


The above table lists many of the possible optimizing 
characteristics for the ship and its propulsion sub-system. 
Manning estimates are typically based on historical data and 
do not indicate that the baseline and variant ships will 
require a significantly different number of personnel. Com¬ 
monality measures the use of the same equipment in other 
ships. Since there are no other electric drive ships at the 
power levels used in this thesis, commonality is not an 
issue. It is very difficult to discern maintainability and 
reliability differences between designs that are as close as 
the machine designs of this thesis, so these two charac- 
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teristics were not made part of the objective function. 

Cost is a measure that should be part of this sort of op¬ 
timization. The cost of the ship system, quantified in the 
ship displacement, was used in the final recommendations for 
the transmission sub-system. In the case of permanent mag¬ 
net machines, the relative costs of magnet material and mag¬ 
net steel were included in the objective function. 

Weight, volume, and efficiency were made an explicit 
part of the objective function for the computer design of 
the machines. An Effective Weight was calculated for ovary 
machine design. The design with the lowest effective weight 
was the "best" within its class. The generic objective 
function Is 

Effective Weight = weight + ke*(1-efficiency) + kv*volume 

where ke and kv are weighting factors for efficiency and 
volume, respectively. The weighting factors were obtained 
from changes in ship displacement for marginal changes in 
efficiency and volume of the transmission. They were 
modified to reflect the actual designs resulting from the 
process. 
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Chapter Two. 


General Considerations 


Only 3teady-state behavior of electric machine:, was 
modeled. The modeling of dynamic behavior is very difficult 
and was not viewed as being within the design problems posed 
by this thesis. The changes in machine design necessary to 
solve dynamic instabilities, etc., are much smaller than the 
approximate nature of the algorithms used here. 

All derivation work was performed without specifying 
the number of winding turns or the number of rotor or stater 
slots. The only exception to this was the case of induction 
machines, where an arbitrary number of rotor slots was 
selected. This selection was necessary for the calculation 
of the equivalent circuit components. The number of rotor 
slots chosen, 71, was a number designed not to induce pole 
harmonics. Since no turn numbers were specified, units in¬ 
clude volts-per-turn, ampere-turns, and impedance-per-turns- 
squared. Power is measured in watts. 

All machines used as their synchronous frequency the 
maximum allowed by the particular combination of pole pairs 
and shaft rpm. Developmental work showed that the optimisa¬ 
tion algorithm converged to the highest frequencies, so the 
algorithm now starts at the highest possible frequencies. 

Up to six pole pairs were used in the higher rpm 
machines and up tc 25 pole pairs in the 180 rpm machines. 
Diminishing improvements in volume, weight, and efficiency 
show up at half these limits. 

The random number generator was taken from Kelley and 
Pohl [14], with one change. After every run of each 
program, the random number generator seed is stored. This 
means the sequence of pseudo random numbers is not repeated 
until the full range possible has been used. For the 
machines, it gives differences at every run and means the 
multidimensional variable spaoe is more fully explored. 
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2.1. Optimization method 


The chosen optimization method is a combination of the 
Monte Carlo and steepest-descent schemes. A design point is 
established by randomly selecting machine geometric 
parameters, subject to constraints. Ten random steps are 
taken around the design point, in all variable directions. 
The effective weight of each random step is evaluated and 
compared with that of the design point. The best of the 
eleven is designated the new design point. More random 
steps are taken, and the process continues until no more im¬ 
provement is seen in effective weight. At that point, the 
size of the random steps is halved, and the process repeats 
itself, with the step size continually halved (up to ten 
times). The best effective weight is the index to the best 
design. The number of original design points used in any 
particular run of a program is under user control. If ten 
original design points are desired, the algorithm will look 
at over a thousand designs. 

The purpose of having original design points ia to 
start the optimization process in different sectors of the 
multidimensional variable space. In this fashion, the op¬ 
timization process zeroes in on several local "beat" points. 

The variables that are randomly selected include stator 
current density, rotor radius, air gap dimension, stator 
slot apace factor, and rotor slot space factor. The back 
iron dimension (the iron behind the stator teeth) is sized 
to handle a saturation level of flux. The stator slot depth 
is originally sizud as a xandom fraction of the back iron 
dimension, and the rotor slot depth is originally a random 
fraction of the rotor radius. 

Only steady-state behavior was modeled in this algo¬ 
rithm. Dynamic modeling may or may not show different op¬ 
timum configurations for machine types. 
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2.2. Constraints 


The constraints placed on the optimization process are 
listed in the following table. The moat difficult con¬ 
straint to satisfy while still achieving a valid design was 
the rotor current density constraint in induction machines. 
Only a few valid designs were achieved in induction machines 
using the above algorithm, leaving some question about the 
application of the algorithm in the case of induction 
machines. Only those induction machine designs in which 
there was reasonable confidence were included in the thesis 
analysis. 


Table 3. Optimization constraints during machine design 


Minimum air gap flux density 
Maximum (saturation) flux density 
Maximum rotor radius 
Maximum rotor tip speed 
Stator and rotor space factor 
Maximum rotor slot depth 
Maximum synchronous reactance: 
synchronous machines 
permanent magnet machines 
Power factor 


1.05 tesla rms 
1.5 tesla rms 
2.0 meters 
200.0 meters/sec 
0.35 

33*5 of rotor radius 

2.0 per unit 
3.0 per unit 
0.8 


The magnet steel chosen was 26 gauge M19. Its magnetic 
properties were found in USX technical data [15]. It has 
been observed that saturation flux levels in electric 
machines occur first where the area perpendicular to the 
flux path is the smallest. If the baok iron dimension is 
made appropriately large, this saturation will first occur 
in the teeth, as is desirable. Accordingly, the back iron 
dimension was set to 


dcore 


Br r 
Bsat p 


whore Br is the radial air gap flux density, r is the rotor 
radius, Bsat is the saturation flux density, and p is the 
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This equation is 


number of pole pairs in the machine, 
derived from Gauss’ Law. 

If electric machines are to be installed in a ship, th 
78 ey obviously will need to fit into the space designated 
for them. An electric machine with a two meter rotor radius 
is at least thirteen feet in envelope diameter. This is a 
very large machine to install in a machinery space where 
volume is already at a premium. The rotor radius limit de¬ 
scends from the physical ability to fit an electric machine 
in a ship. 

The tip speed limitation represents the physical limit 
on material strength with regard to the rotor conductors. 
Rotor conductors may break free from the rotor at higher 
tangential velocities than this limit. The magnitude of the 
limit was taken from a tip-speed-limited, 3600 rpm, two-pole 
turbogenerator, and was verified against standard Navy 
design practices. Several as-built electric machines were 
analyzed and this number seemed to fit their characteristics 
well. The tip speed limit arises when choosing a rotor 
radius (given a particular frequency and number of poles), 
and is less stringent a constraint than maximum rotor 
radius. 

Thermal considerations are often extremely important in 
machine design. The heat build-up in electric machines has 
led to many cooling schemes over the years, including 
natural convection, forced air cooling, and hydrogen cool¬ 
ing. One of the latest methods is liquid coding of the 
stator and rotor conductors through cooling passages through 
the copper itself. This has been made possible by better 
de-ionizing methods for cooling fluid and better rotating 
seals for the rotor. Naturally, the cooling passages and 
insulation limit the amount of copper area in a slot cross 
section. The copper area in a typioal conductor bar was 
measured and found to be about thirty-five percent of the 
bar cross section. This number was used for the stator and 
rotor slot space factors 
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Other thermal considerations must be made for permanent 
magnet materials, which suffer a degradation in flux as tem¬ 
peratures rise. Flux loss rises slowly with increasing tem¬ 
perature until about 100®C. Above 100®C, flux loss is more 
rapid. An assumption in the design of these machines is 
that there will be sufficient cooling in the operating space 
to limit ambient temperature to about 80*C. This, combined 
with the machine liquid cooling, should keep flux loss to a 
minimum. Transmission lines were assumed to function satis¬ 
factorily at the same temperature. 

Insulation also has a thermal rating. No insulation 
class was specified in this thesis but a typical insulation 
used in electric generators by the Navy is Class F. For 
this class, a permissible rise of 100®C over an ambient tem¬ 
perature of 50®C is standard, but lesser insulation olasses 
must run cooler. If a machine must be designed with a les¬ 
ser insulation class, the consequent lower temperatures may 
result in a quieter machine and longer machine and insula¬ 
tion life. It probably will be larger than a machine with a 
greater class of insulation. The temperatures quoted above 
are at hot spots, not in the bulk of the machine.il 

Along with the reduction in copper area for liquid 
cooling, a maximum current density was imposed. The copper 
losses, in the form of heat, have to be removed by the cool¬ 
ing fluid. There is a tradeoff between the size of the 
cooling passages, the allowable current density, and the 
rating of the machine. Twelve million amperes per square 
meter equates to forty amperes in a twelve gauge copper 
wire. 

Rotor slots were constrained to no more than one third 
of the rotor radius. Some reasonable shaft diameter is 
required to transmit the mechanical torque. Stator slots 
were allowed to grow as needed to meet the stator current 


11. Private communication, D.F. Schmuckex, Naval Sea Systems 
Command, Code &6Z31, 11 Decembor 1986. 
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density limit. 

Gynchror-ous reactance limits were taken from as-built 
machines. A power factor of 0.8 was used for all machines, 
though one researcher has indicated a power factor of 1.0 
might be best for permanent magnot machines. 12 

2.3. Geometric considerations 

End turns were modeled as described in Appendix B. 

While not exactly as machines are constructed, this model 
gives reasonable results. A length allowance equal to one 
rotor diameter on either end of the active length of the 
machine was made to allow for containment of the end turns. 

Fractional slot pitches were not considered. A stator 
winding pitch of 0.8 was assumed, resulting in the elimina¬ 
tion of the fifth harmonic from the steady-state output 
waveforms. 

Once weight and volume were calculated, an extra ten 
percent was added to allow for the frame and foundation of 
the machine. The calculated weight included an additional 
three percent of the rotor weight to allow for bearings. It 
is the final envelope weight and volume that were used in 
the decision process. 

2.4. Efficiency and losses 

The general equation for efficiency is 
minpwr 

efficiency = ----■----- 

minpwr + ph + pe + i2r + i2rr 

where minpwr is the minimum mechanical power expected of the 
machine, ph is the hysteresis loss, pe is the eddy current 


12. Robey, H.N., "Permanent Magnet Machine Technology 
Assessment," DTN6RDC Report TM-27-80-87, September 1980. 
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loss, 12r is the stator copper loss, and 12rr Is the rotor 
copper loss. This formulation is for a motor, but the ef¬ 
ficiency calculated will not be significantly different if 
the machine is a generator. 

Hysteresis and eddy current losses arise from currents 
circulating within the magnet steel that forms the rotor and 
stator. They are two different mechanisms and depend on the 
metallurgy of the steel. 

Eddy currents are a result of the time-varying magnetic 
fields within the machine, and they oppose the change in 
flux density within the machine. Eddy current losses in¬ 
crease as the square of the electrical frequency of the 
machine and also as the square of the peak flux density. 

One method of lessening eddy current losses is to use thin 
laminations to build up the rotor and stator. If the var¬ 
nish used on the laminations is sufficiently insulating, the 
eddy currents are limited to azimuthal circulation. Axial 
circulation is practically zero because of the small lamina¬ 
tion thickness. 

Hysteresis losses are inherent to magnetic materials, 
and are proportional to the total volume of the material, 
the area of the hysteresis loop, and the machine electrical 
frequency. 

USX has developed equations to calculate eddy current 
and hysteresis losses. They are 


0.01445 fl f Br He 


and pc 


0.4818 n Bm» t* f* 
rho D 


In these aquations, 8 is the hysteresis loss factor (the 
ratio of the actual hysteresis losses to the area of a 
square hysteresis loop passing through Br and Ha), f is the 
frequency in Hertz, Br is the residual induction in 
kilogauss, Ho la the coercive force in oersteds, D is the 
density in grams per cubic centimeters, rho is the electri¬ 
cal resistivity of M19 in miorohm-cm, D is the anomalous 
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loss factor. The losses are In watts per pound of material. 
The numerical factors at the beginning of both equations 
were altered to reflect the use of SI units. The factors li , 
n, Br, He, rho, and D change with the type of magnetic 
material used. The thickness of the laminations, t, is 
0.014 inches for 26 gauge steel. These equations do not 
reflect some variations caused by differences in silicon 
content and differences in processing treatments leading to 
variations in grain size and crystallographic texture. 
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Chapter Three. Synchronous machines 


Synchronous machines operate because of an interaction 
between stator and rotor flux waves. The rotor flux wave is 
developed by a field winding. The stator flux wave is 
developed by an armature winding. These two waves try to 
align themselves, which is how the machine action is 
produced. In the case of a motor, the armature wave is 
"rotating" around the periphery of the stator bore because 
of the 120° separation between the three phases. A rotating 
action ensues. In the case of a generator, the rotation is 
provided by a prime mover, such as a gas turbine, and volt¬ 
age is induced in the stator phases. 

Synchronous machines operate at a steady-state shaft 
speed specified by the number of poles and the electrical, 
or synchronous, frequency. This synchronous speed is main¬ 
tained despite changes in load. This feature makes 
synchronous machines attractive for applications where speed 
control is Important. Shaft rpm = (120*frequency)/(poles). 

A derivation of the equations of synchronous machines 
and the computer modeling program are presented in 
Appendix B. 

3.1. Assumptions 

The rotors of synchronous machines may exhibit 
saliency, or may be smooth cylindrical rotors. The dif¬ 
ferences in properties and parameters among salient and 
round-rotor machines amount to only a few percent* *. The 
approximate nature of the modeling means the saliency ef¬ 
fects will not be important. Therefore, no special provi¬ 
sion for salient rotors were made. 

The ships used in this thesis have displacements of 


13. Fitzgerald, Kingsley, and Umans, Electric Machinery. 
McGraw-Hill, New York, 1983. 
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about 5000 LT ( 1 LT = 2240 lba). Their fuel load is deter¬ 
mined by the required range and transmission efficienoy at 
endurance speed. It was first estimated that a one-percent 
increase in transmission efficiency would reduce total ship 
displacement by about 89 LT. The efficiency factor (ko = 
90,000 kg/percent) corresponds to this 89 LT change in ship 
full load displacement. (89 LT x 2240 lbs/LT x 2.205 kg/lb) 
When this produced machines with efficiencies about 95X, it 
was doubled to 180,000. Obviously, this factor may be ad¬ 
justed to any level. The volume efficiency faotor (kv = 
1286.1 kg/m*) corresponds to the density of a LM-2500 gas 
turbine module. These factors were used throughout the 
thesis. 

3.2. Machine description 

Machines with shaft speeds of 1800, 2400, 3000, 3600, 
and 7200 rpm, with the number of pole pairs varying from one 
to six, were modeled. Also, 180 rpm machines using from one 
to twenty-five pole pairs were modeled. This provided a 
good coverage of the variable space. 

3.2.1. Efficiency 

Synchronous machine efficiency at full load was about 
98.5* for the higher rpm machines, while the 180 rpm 
machines hovered around 93% to 94k efficiency. The number 
of pole-pairs seemed to have little effect in the high rpm 
machines, but there was am "arch" in the efficiency curve of 
the 180 rpm machines, peaking at 97k with 36 poles. Though 
not fully understood, the 24- and 26-pole 180 rpm machines 
had very low efficiencies. Generally speaking, machine ef¬ 
ficiency was higher when rpm was higher. This was an ex¬ 
pected result. Off-design-point efficienoy was good for the 
higher rpm machines, but bad for the 180 rpm machines. (See 
the discussion in Chapter Seven.) 
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3.2.2. Weight and volume 


Weight and volume increased almost linearly with the 
number of pole-pairs in the ISO rpm machine. This is almost 
certainly a function of tip-speed limitations, as the maxi¬ 
mum tip speed in a machine is a function of the number of 
poles in the machine. When the rotor radius is limited, the 
machine must grow in length to develop enough torque. The 
weight and volume of these machines were much higher than 
for the higher rpm machines. 

The higher rpm machines saw significant decreases in 
weight and volume when the number of pole-pairs increased 
from one to two. There also was an observable increase in 
weight and volume as the number of pole-pairs further in¬ 
creased. As rpm increased, the machines became smaller and 
lighter. 
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Figure 3. High rpm synchronous machine volume 
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Figure 4. 180 rpm synchronous machine volume 


SYNCHRONOUS MOTOR VOLUME 
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(spuosnotii) 


Figure 5. High rpm synchronous machine weight 


SYNCHRONOUS MOTOR WEIGHT 
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Figure 6. 180 rpra synchronous machine weight 


SYNCHRONOUS MOTOR WEIGHT 


DIRECT-DRIVE, 180 RPM, 25775 HP 



NUMBER OF POLE PAIRS 
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3.3. Machine characteristics 


The following tables give machine characteristics for 
many cf the higher rpm and the 180 rpm synchronous machines 
designed for this thesis. 

The stator slot factor tended to increase to the limit 
of 0.75, while the rotor slot, factor moved around the value 
0.58 quite a bit. This demonstrates a partial limit on the 
depth of the stator slots (to the same dimension as the back 
iron depth). The overall diameter was limited and more 
stator slot area was needed to develop the required power. 
Similarly, the stator current density converged to its cook¬ 
ing limit. 

The longest higher rpm machine was 5.42 meters, while 
the largest overall diameter was 0.85 meters. Machines of 
this size will causa no difficulties when placed in the 
machinery spaces of most ships. The 180 rpm machines are 
larger, typically less than 6 meters long (discounting one 
16 meter machine) and 1.8 meters in diameter. They are also 
good candidates for ahlp systems. 
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Table 4. Characteristics of 1800 rpm synchronous machines 


Miibtr of poll pairs 

1 

2 

3 

a 

3 

6 

pouer, hp 

23773 

23775 

25773 

23773 

25775 

23773 

efficiency Victor 

180000 

180000 

180000 

180000 

180000 

I8C000 

voluae factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rp* 

teoo 

1800 

1800 

1800 

1800 

1800 

stator current density 

B.0CE*06 1.20£»07 

1.20E+07 i.20E*07 

1.2CE+07 

1.20E«07 

synchronous friqumcy 

30 

60 

90 

120 

150 

180 

rotor radius 

0.1720 

0.2077 

0.2603 

0.2905 

0.3267 

0.3189 

gap dimension 

0.0138 

0.0440 

0.0250 

0 0124 

0.0090 

0.0092 

back iron depth 

0.1204 

0.0727 

0.0607 

0.0508 

0.0437 

0.0372 

stator slot depth 

0.0772 

0.0722 

0.0607 

0.0508 

0.0437 

0.0372 

rotor slot depth 

0.0401 

0.1859 

0.0747 

0.1333 

0.0907 

0,0289 

stator slot factor 

0.633 

0.730 

0.750 

0.749 

0.750 

0.748 

rotor slot factor 

0.630 

0.447 

0.687 

0.639 

0.449 

0.628 

•nvilape valun 

2.781 

1.522 

1.509 

1.541 

1.697 

1.709 

envelope weight 

20835.40 

9827.96 

8766.09 

9186.03 

8661.17 

BR33.21 

hysteresis loss 

16926.49 10677.39 

17111.27 19469.83 31370.98 

43892.02 

eddy current loss 

31-33.2 

3826.7 

9198.B 

13955.6 

28107.7 

47191.5 

stator copper lass 

273834.4 374646.5 

293335.3 164289.8 173116.1 

269901.7 

full load efficiency 

0.985 

0.980 

0.9B4 

0.990 

0.9B8 

0.982 

active length 

4.670 

1.794 

1.498 

1.313 

1,348 

1.741 

full load current density 

I.19E+07 1.06E+07 

1.17E+07 4.58E+06 B.34F+06 

1.89E*07 

no load current density 

4.40E+06 6.39EH16 

7.00E+06 

2.43E*06 4.24E+06 

1.196+07 

xs/turns-squarad, p.u. 

1.99 

0.B6 

0.87 

1.11 

1.20 

0.78 

internal volts/turn, p.u. 

2.71 

1.66 

1.67 

1.89 

1.97 

1.60 

overall length 

5.42 

2.80 

2.64 

2.72 

2.69 

3.03 

overall diaieter 

0.77 

0.79 

O.Bl 

0.81 

0.83 

0.80 


Table 5. Characteristics of 2400 rpm synchronous machines 


nuaber of polo pairs 

1 

2 

j 

4 

5 

8 

poser, iip 

23775 

23775 

23775 

25775 

23775 

23773 

efficiency factor 

180000 

180000 

190000 

180000 

180000 

180000 

voluio factor 

1288.1 

1288.1 

1288.1 

1288.1 

1288.1 

1288.1 

shaft rp* 

2400 

2400 

2400 

2400 

2400 

2400 

itator currant density 

1.15E*07 

1.20E+07 

1.19E+07 

1.20E+07 

1.20EtO7 

1.07ES07 

synchronous froqutncy 

40 

BO 

120 

180 

200 

240 

rotor radius 

0.1383 

0.1988 

0.2388 

0.2890 

0.3141 

0.3J97 

pap diientlor. 

u.0273 

0.0184 

0.0089 

0.0135 

0.0109 

0.0088 

back iron depth 

0.1094 

0.0893 

0.0357 

0.0508 

0.0440 

0.0394 

stator slot depth 

0.0929 

0.0895 

0.0337 

0.0505 

0.0434 

0.0398 

rotor slot depth 

0.0418 

0.0213 

0.0088 

0.0583 

0.1047 

0.1040 

stator slat factor 

0.821 

0.750 

0.750 

0.748 

0.750 

0.750 

rotor slot factor 

0.878 

0.749 

0.750 

0.432 

0.433 

0.327 

envelope voluse 

l-i.lt 

1.112 

1.133 

1.329 

1.438 

1.833 

envelope might 

11819.57 

7204.42 

8877.13 

8733.44 

7108.43 

7727.01 

hysteresis loss 

11981.17 

14511.73 21941.41 27815.01 

32789.75 

47103.39 

eddy current luss 

2882.8 

8934.5 

•5727.2 

28391.9 

39147.9 

47525.7 

stator copper loss 

444805.4 427187.9 482231.4 232958.2 

143185.8 

114790.8 

full luad efficiency 

0.977 

0.977 

0.974 

0.985 

0.989 

0.988 

active length 

J.397 

1.709 

1.352 

1.152 

1.144 

1.223 

full load current density 

2.1 iE+0/ 3.11E+07 5.J9E*07 l.73E*07 

7.83E+08 

7.99E*04 

no load current density 

7.7'5E*08 1.2BE+07 2.11E+07 9.84E+08 

4.55E+08 

4.32E+04 

ii/turns squared, p.u. 

2.00 

1.89 

1.83 

1.01 

0.92 

1.07 

internal volts/turn, p.u. 

2.72 

2.43 

2.58 

1.80 

1.72 

1.85 

overall length 

3,13 

2.57 

2.54 

2.38 

2.44 

2.81 

overall diaoeter 

0.77 

0.71 

0.72 

0.81 

0.82 

0.B3 
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Table 6. Characteristics of 3000 rpm synchronous machines 


nuiher o< pole pairs 

1 

2 

3 

4 

5 

6 

pmer, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

100000 

130000 

180000 

180000 

180000 

voluie factor 

1288.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rp» 

5000 

3000 

3000 

3000 

3000 

3000 

stater current density 

l.20E*07 1.20E+07 1.20E+07 

1,20E*07 

1.10£*07 

1.18E+07 

synchronous frequency 

SO 

100 

150 

200 

250 

300 

rotor radius 

0.1494 

0.1890 

0.2304 

0.2432 

0.2872 

0.2974 

qap dimension 

0.0370 

0.0172 

0.0101 

0.0052 

0.0068 

0.0088 

back iron depth 

0.1048 

0.0662 

0.0538 

0.0426 

0.0402 

0.0347 

stator slot depth 

0.0882 

0.0628 

0.0537 

0.0422 

0.0402 

0.0347 

rotor slot depth 

0.0904 

0.0879 

0.0372 

0.0151 

0.0784 

0.0427 

stator slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

0.750 

rotor slot factor 

0.598 

0.494 

0.724 

0.730 

0.255 

0.750 

envelope voluie 

1.205 

0.963 

0.972 

0.991 

1.203 

1.236 

envelope Height 

8815.49 

6498.7? 

5801.99 

5725.40 

6176.39 

6143.64 

hysteresis loss 

9574.79 14072.41 21089.10 31331.81 41078.93 46170.48 

eddy current loss 

2859.8 

8405.7 

18895.4 

37429.9 

61340.0 

82749.7 

stator copper loss 

380285.7 212841.? 

183124.5 204027.7 140446.5 134484.2 

full load efficiency 

0.980 

0.988 

0.989 

0.986 

0.988 

0.986 

active length 

1.870 

1.657 

1.363 

1.590 

1.306 

1.310 

full load current density 

1.53£f07 1.13E+07 

l.3BE*07 2.42E+07 

1.40E+07 

1.06E*97 

no load current density 

6.03E+06 5.238*08 8.09E>08 9.63E+06 7.44E+06 6.96£*0u 

xs/turns-squared, p.u. 

1.80 

1.41 

1.53 

1.78 

1.11 

0.69 

internal volts/turn, p.u. 

2.53 

2.16 

2.27 

2.51 

1.89 

1.52 

overall length 

2.42 

2.48 

2.32 

2.58 

2.48 

2.53 

overall diaieter 

0.76 

0.67 

0.70 

0.67 

0.75 

0.75 
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Table 7. Characteristics of 3600 rpm synchronous machines 


nuaber of pole purs 

1 

2 

3 

4 

5 

6 

poner, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

voluee factor 

1236.1 

1286.1 

1286.1 

1286.1 

1286.1 

1236.1 

shaft rpe 

3600 

3600 

3600 

3600 

3600 

3600 

stator current density 

i.16E+07 

1.20E+07 

1.208*07 

1.19E+07 

1.20E+07 

1.20E+07 

synchronous frequency 

60 

120 

180 

240 

300 

360 

rotor radius 

0.1413 

0.1824 

0.2191 

0.2736 

0.2846 

0.2877 

gap diiension 

0.0139 

0.0120 

0.0123 

0.0178 

0.0209 

0.0124 

back iron depth 

0.0939 

0.0638 

0.0511 

0.0479 

0.0398 

0.0336 

stator slot depth 

0.0492 

0.0633 

0.0486 

0.0479 

0.0398 

0.0336 

rotor slot depth 

0.0496 

0.0541 

0.0572 

0.0298 

0.0595 

0.03B5 

stator slot factor 

0.748 

0.750 

0.750 

0.750 

0.750 

0.750 

rotor slot factor 

0.733 

0.409 

0.528 

0.679 

0.691 

0.654 

envelope voluae 

1.233 

0.814 

0.870 

1.065 

1.127 

1.102 

envelope weight 

93B5.14 

5377.73 

5216.22 

4979.26 

5177.33 

5190.73 

hysteresis loss 

14030.95 

15587.02 

22251.57 

31042.24 

36284.98 

48590.52 

eddy current loss 

5028.6 

11172.5 

23924.3 

44501.0 

65021.1 

104436.5 

stator copper loss 

190131.9 

228605.1 

165930.4 

251863.9 

196473.3 

184008.6 

full load efficiency 

0.989 

0.987 

0.989 

0.983 

0.985 

0.9B3 

active length 

3.253 

1.493 

1.372 

0.890 

0.975 

1.166 

full load current density 

9.15E*06 

1.94E+07 

1.31E+07 

2.47E*07 

1.45E+07 

1.73E*07 

no load current density 

3.3UE+06 

7.41E+06 

6.99E+06 

1.61E*07 

1.12E+07 

1.28E+07 

xs/turns-squared, p.u. 

1.99 

1.89 

1.10 

0.71 

0.43 

0.4B 

internal volts/turn, p.u. 

2.71 

2.62 

1.83 

1.54 

1.30 

1.35 

ove all length 

3.88 

2.27 

2.30 

2.06 

2.20 

2.37 

overall disaster 

0.61 

0.64 

0.66 

0.77 

0.77 

0.73 
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Table 8. Characteristics of 7200 rpm synchronous machines 


nueher of pole pairs 

1 

2 

3 

4 

5 

6 

poser, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

voluu factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

7200 

7200 

7200 

7200 

7200 

7200 

stator current density 

1.20€‘07 

1.20E*07 

1.16E+07 

I.20E+07 

1.20E+07 

1.20E+07 

synchronous frequency 

120 

240 

360 

480 

600 

720 

rotor radius 

0.1159 

0.1710 

0.1958 

0.2044 

0.2219 

0.1932 

qap diMnsion 

0.0096 

0.0093 

0.0053 

0.0093 

0.0046 

0.0055 

back iron depth 

0.0812 

0.0599 

0.0457 

0.0358 

0.0311 

0.0225 

stator slot depth 

0.0405 

0.0452 

0.0457 

0.0358 

0.0311 

0.0225 

rotor slot depth 

0.0329 

0.0189 

0.0105 

0.0714 

0.0138 

0.0292 

stator slot factor 

0.750 

0.750 

0.749 

0.750 

0.750 

0.750 

rotor slot factor 

0.511 

0.617 

0.661 

0.633 

0.468 

0.570 

envelope voluie 

0.711 

0.544 

0.525 

0.544 

0.572 

0.564 

envelope neight 

5341.21 

3476.08 

2971.92 

3136.32 

2950.31 

3427.07 

hysteresis loss 

17305.86 22097.70 28944.03 30808.73 

50465.33 

64447.34 

eddy current lass 

12404.5 

31678.5 

62239.7 

83332.6 

180863.2 

277168.2 

stator copper loss 

149625.0 159505.1 162532.1 

74543.2 

143253.5 

117812.4 

full load efficiency 

0.991 

0.989 

0.987 

0.990 

0.981 

0.977 

active length 

2.863 

1.212 

0.972 

1.081 

1.082 

1.952 

full load current density 

1.65EV07 2.71E+C7 3.92E+07 

7.81E«Q6 

3.70E*07 

1.72E*07 

no load current density 

6.14E+06 

1.16E+07 

I.46E+07 

5.00E+06 

1.99E+07 

1.28E+07 

xs/turns-squared, p.u. 

1.97 

1.59 

1.97 

0.74 

1.08 

0.48 

internal volts/turn, p.u. 

2.69 

2.33 

2.69 

1.56 

1.86 

1.34 

overall length 

3.36 

1.93 

1.78 

1.94 

1.99 

2.75 

overall diaeeter 

0.49 

0.57 

0.58 

0.57 

0.58 

0.49 
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Table 9. Characteristics of 180 rpm synchronous machines 


nuibtr of poll pairs 

1 

2 

3 

4 

5 

pour, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180040 

180000 

IBOOOO 

180000 

180000 

voluee fictor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shift rpe 

180 

180 

180 

180 

180 

stator current density 

1.1BE*07 

1.12E+07 

1.20E+07 

9.25E*06 

1.20€*07 

synchronous frequency 

3 

6 

9 

12 

IS 

rotor ridius 

0.2641 

0.3738 

0.4588 

0.5869 

0.5507 

gip dimension 

0.0331 

0.0978 

0.0723 

0.0317 

0.0198 

hick iron depth 

0.184V 

0.1308 

0.1071 

0.1027 

0.0771 

stitor slot depth 

0.05V2 

0.1308 

0.1071 

0.1027 

0.0580 

rotor slot depth 

0.0500 

0.1867 

0.1934 

0.0565 

0.1079 

stitor slot fictor 

0.741 

0.750 

0.750 

0.742 

0.750 

rotor siot fictor 

0.558 

0.571 

0.422 

0.535 

0.287 

envelope voluee 

16.305 

9.329 

9.048 

11.378 

10.306 

envelope eeight 

119110,5 

56300.0 

48896.8 

56398.0 

56B99.5 

hysteresis loss 

10465.3 

6915.1 

9833.5 

18267.4 

23999.1 

eddy current loss 

187.5 

-47.8 

52B.6 

1309.4 

2150.3 

stitor capper loss 

2.25E+06 2.06EO6 

I.76E+Q6 1.34E+06 1.1BE+06 

full load efficiency 

0.895 

0.903 

0.916 

0.934 

0.941 

active length 

14.913 

3.135 

2.590 

2.376 

3.708 

full load current density 

1.53E+07 1.25E+07 

1.336*07 1.91E+07 1.42E+07 

no load current density 

5.61E*06 6.14E<06 7.25E06 tt.93E*06 7.27E+06 

xs/turns-squared, p.u. 

2.00 

1.28 

1.04 

1.38 

1.19 

internal volts/turn, p.u. 

2.72 

2.04 

1.83 

2.14 

1.96 

overall length 

16.10 

5.02 

4.71 

4.85 

5.99 

overall diaeeter 

1.08 

1.47 

1.49 

1.65 

1.41 
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Table 10, Characteristics of more 180 rpm synchronous 
machines 


nuaber of polt pairs 

6 

7 

a 

9 

10 

poMr, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

voIum factor 

1281.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

180 

180 

1B0 

180 

180 

stator currant density 

1.20E+07 

I.20E+07 

1.20E+07 1.20E+07 

1.20E*07 

synchronous frequency 

IB 

21 

24 

27 

30 

rotor radius 

0.6171 

0.5729 

0.7152 

0.6943 

0.7296 

pap diaansian 

0.0647 

0.0335 

0.0229 

0.0138 

0.0396 

back iron depth 

0.0720 

0.0573 

0.0626 

0.0540 

0.0511 

stator slot depth 

0.0720 

0.0569 

0.0626 

0.0540 

0.0509 

rotor slat depth 

0.1652 

0.0855 

0.2046 

0.0740 

0.1420 

stator slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

rotor slot factor 

0.368 

0.750 

0.374 

0.283 

0.615 

envelope voluie 

11.693 

10.486 

12.870 

12.3B5 

14.273 

envelope weight 

48739.8 

53234.1 

52286.9 

53189.2 

56440.5 

hysteresis loss 

22385.0 

27947.6 

31793.9 

4271B.3 

42086.1 

eddy current loss 

2406.B 

3505.7 

4557.9 

6889.4 

7541.6 

stator copper loss 

1.86E+06 1 .24E+06 7.69E+05 1.14E*06 l.UE+06 

full load efficiency 

0.911 

0.938 

0.960 

0.942 

0.943 

active length 

2.225 

3.418 

2.045 

2.349 

2.365 

full load current density 

1.74E+07 l.UE+07 7.09E+06 1.98E+07 9.69E+06 

no load current density 

1.29E+07 7.98E*06 4.1BE+06 1.07E+07 

7.76E+06 

xs/turns-squared, p.u. 

0.48 

0.53 

0.89 

1.08 

0.36 

internal volts/turn, p.u. 

1.35 

1.39 

1.69 

1.86 

1.25 

overall length 

4.95 

5.84 

5.00 

5.38 

5.44 

overall diaaeter 

1.65 

1.44 

1.73 

1.63 

1.74 
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3.4. Verification 


Data on a large turbogenerator, Big Sandy Unit Two, was 
available in reference sixteen; it was used to verify the 
synchronous machine design program. Big Sandy is rated at 
907,000 kVA, power factor 0,9, at a rated voltage of 26 kV. 
When the machine parameters were input to the synchronous 
design program, it produced a machine very close to Big 
Sandy. It was judged that the design program would yield 
good results. 
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Chapter Four. Permanent Magnet Machines 


Permanent magnet machines are very similar to 
synchronous machines. The main difference lies in the 
method used to produce the field flux wave. Instead of a 
field current causing the wave, permanent magnets provide 
the flux. No exciter is used with these machines and there 
are no brushes. 

4.1. Magnet material 

Many different elements may be used to manufacture per¬ 
manent magnets. Past designs used ceramics, aluminum- 
nickel-cobalt-iron-titanium (AlNiCo), and samarium-cobalt 
(SmCo). However, ceramic magnets do not produce sufficient 
residual flux (see Appendix C for an explanation of terms), 
and any magnet based on cobalt is high in cost and may be in 
limited supply. Recently, magnets of neodymium-iron-boron 
(NdFeB), have entered the marketplace. None of the con¬ 
stituents of the NdFeB magnets are strategic materials; it 
is expected that availability and cost will improve. 

NdFeB magnets have a high Maximum Energy Product (MEP) 
that may be used advantageously by the machine designer. 

Data from Sumitomo Special Metals [17] indicates their 
NEOMAX line have MEPs as high as 37 MGOe, which is higher 
than the 30 MGOe of the SmCo magnets marketed as REC-30 by 
TDK Corporation [18]. High MEP is not the only criteria for 
magnet selection; flux stability, cost-to-performance ratio, 
ease of machine assembly, and other characteristics may en¬ 
ter the decision process. This study needed the best per¬ 
formance of its machines, so NdFeB was selected on the basis 
of its high MEP. Thermal stability was assumed to be satis¬ 
factory if the thermal considerations in Chapter Two were 
met. ' 

Cost has been mentioned above. Magnet material is sig¬ 
nificantly more expensive than magnet steel and copper. 
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Magnets in the quantity used by a large production run of 
25,000 hp machines might cost as much as $120 per pound, 1 * 
compared with the 58$ per pound of M19 steel.* 5 Obviously, 
permanent magnet machines will be more expensive, but the 
cost of magnet materials may be m«do part of the optimiza¬ 
tion process. The degree of magnet overhang, discussed in 
Appendix C, also affects cost. 

4.1.1. Magnet cost factor 

A change to the objective function was made to incor¬ 
porate the cost of magnet material relative to magnet steel. 
The ratio of the above costs was taken and the result called 
the magnet cost factor, km. The objective function was 
modified to 

Effective Wight = (wight ♦ keleagnet wight)) ♦ kell - effcy) ♦ kvlvoluie ♦ keleageet voluee)) 

An initial value for km of 170 was used, and several 
machines designed. Then a value of km = 25 was tried. The 
machines with km = 170 indeed had less magnet material in 
them, but at a cost. The change to km = 25 resulted in a 
larger machine (23.S%) with a lower stator current density, 
20% more magnet material, and about a 1.5% increase in 
machine efficiency. That 1.5% translates to a lot of fuel 
aboard a ship, so it was decided to use km = 25. The extra 
magnet material will add about $22,000 per machine. 

4.2. Assumptions 

The largest obstacle to assembling high-power permanent 
magnet machines is their inherent residual magnetism. If 

14. Estimate by Mr. Yokokura, President of Sumitomo Special 
Metals of America. 

15. Book price for 28 gauge M19 steel from Mr. Dagg of USX. 
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the magnets possess all of their properties at assembly, it 
will be extremely difficult to place the rotor (which con¬ 
tains the magnets) inside the stator. The rotor would be 
strongly attracted to the iron of the stator. Of course, 
the magnets may be magnetized after machine assembly, but it 
may be difficult to achieve MEP without elevated tempera¬ 
tures inside the machine. The assumption is made here that 
the magnets will be magnetized prior to assembly. The 
detailed design of the machine will have to include con¬ 
sideration of the Jigs and fixtures necessary for assembly. 

The only other assumption worthy of mention is that the 
load line may be modeled as described in Appendix C. 

4.3. Machine description 

The same rpm and pole-pair combinations were used as in 
the synchronous machines. The magnets on the rotor are ar¬ 
ranged in a cylindrical-wedge configuration, as shown in the 
figure below. 

Figure 7. End view of a permanent magnet machine 





P146NCT 

8INDER 


The rotor slot factor (lr) used in the other types of 
machines is called here the magnet slot factor and refers to 
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the width of magnet per pole pitch. It may vary between 25% 
and 75% of the pole pitch, the same as for rotor slots in 
other machines. 

The rotor slot depth, ds, does not exist in this 
machine. In this case, the rotor radius, r, is added to lm, 
the magnet radial dimension, to find the actual width of the 
rotating core. (In the synchronous and induction machines, 
dr was included in the rotor radius.) 

4.3.1. Efficiency 

For the higher rpm machines, efficiency within a par¬ 
ticular rpm group decreased with an increasing number cf 
poles. The most efficient machines, at about 99%, had four 
poles. This is higher than the synchronous or induction 
machines, largely because there are no rotor copper losses. 
The twelve-pole efficiencies were about 98%, which is not 
too large a spread. 

The 180 rpm machines had a fairly flat efficiency curve 
(excepting one anomaly) up to about a 26-pole machine, where 
efficiency started to vary widely. There, the conflict be¬ 
tween the number of pole-pairs and maximum rotor radius 
started to become significant. The flat efficiency was 
about 96.5%, which is less efficient than the higher rpm 
machines. 

4.3.2. Height and volume 

The higher rpm machines had a general tendency toward 
lower weight and volume as rpm increased. Within an rpm 
group, the four- and six-pole machines had the lowest weight 
and volume. The smallest machines were larger than the 
synchronous machines. 
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Figure 9. 180 rpm permanent magnet efficiency 
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Figure 10. High rpm permanent magnet maohine volume 
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Figure 12. High rpm permanent magnet machine weight 
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The weight and volume of the 180 rpm machines increased 
with the number of pole-pairs. Again, the tip speed limita¬ 
tion is linked to this increase, as with the synchronous 
machines. There was a wide variation around the general in¬ 
crease, and these maahines are large relative to all others. 
They are not competitive as ship propulsion motors because 
of their size. 

4.4. Machine characteristics 

The following tables give the maohine characteristics 
of the permanent magnet machines designed for this study. 

The stator slot factor, Is, tended to the maximum of 0.75 
for the higher rpm and 180 rpm machines. This ocourred as 
the optimization algorithm tried to minimize envelope volume 
and weight. The rotor slot factor, Ir, was a constant 0.378 
for all machines. The rotor slot factor was calculated to 
produce load-line (MEP) operating flux, as derived in Appen¬ 
dix B. With a different magnet material selection (and a 
consequent change in operating point flux), a different lr 
would have resulted. 

The magnet overhang tended toward the maximum limit, in 
an attempt to achieve the highest flux levels. Permanent 
magnet machines cannot rival the flux level produced by the 
field winding of a synchronous machine, but the optimization 
algorithm did its best. 

The per-unit synchronous reactance-per-turns-squared 
was limited to 3.0 in these machines. This reactance tended 
toward the limit, but was lower with an increasing number of 
pole-pairs in the 180 rpm machines. It was very difficult 
to achieve valid designs with xseax = 2.0, as in the 
synchronous machines. If an xsmi greater than 2.0 is unac¬ 
ceptable, these machines will be less competitive. 

The amount of magnet material varied from 50 kg to a 
few hundred kg in the higher rpm machines to 700-4000 kg in 
the 180 rpm machines. The cost of 700 kg of NdFeB is about 
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T«bl« 11. Cherecterietic* of 1800 rpm magnet machines 


nunber of polo pair* 

1 

2 

3 

4 

5 

6 

poeer, hp 

25773 

25773 

23773 

23773 

23773 

25773 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

voIum factor 

1284.1 

1234.1 

1286.1 

12B6.1 

1286.1 

1286.1 

aagnet factor 

23 

23 

23 

23 

23 

23 

diaft rpa 

1840 

1800 

1800 

1800 

1800 

1800 

itator currant density 

7.49E*06 

1.20E+07 

7.32E+06 

9.14E+06 

I.00E+06 

1.04E+07 

lynchranoui frequency 

30 

40 

90 

120 

130 

180 

rotor radiui 

0.2887 

0.3396 

0.3643 

0.3964 

0.3344 

0.6068 

gap diatntion 

0.0324 

0.4382 

0.0047 

0.0188 

0.0067 

0.0123 

back iron depth 

0.1332 

0.0870 

0.0530 

0.0742 

0.0302 

0.0482 

stator slat depth 

0.0487 

0.0742 

0.0436 

0.0333 

0.0263 

0.0277 

aagnet radial diatnsion 

0.0342 

0.0401 

0.0070 

0.0198 

0.0071 

0.0130 

stator (lot factor 

0.730 

0.730 

0.750 

0.730 

0.750 

0.646 

rotor dot factor 

0.378 

0.378 

0.378 

0.378 

0.378 

0.378 

envelope voluae 

3.348 

2.338 

3.309 

6.304 

5.847 

6.319 

onvtlopa weight 

24280.7 

12328.8 

19366.5 

20910.3 

24826.4 

20948.1 

hysteresis loss 

20381.5 

20298.5 

52315.8 

80303.3 

121493.4 

124368.0 

eddy current loss 

3488.1 

7274.8 

28124.3 

57704.9 

106837.1 

133717.2 

stator copper loss 

171130.0 

265272.5 

93334.8 

98867.0 

69323.1 

85961.4 

full load efficiency 

0.990 

0.983 

0.991 

0,988 

0.985 

0.982 

active length 

2.328 

0.876 

2.833 

0.933 

1.883 

1.237 

■agnet wight 

403.184 

262.376 

157.335 

273.763 

178.613 

246.460 

aagnet voluue 

0.055 

0.035 

0.021 

0.037 

0.024 

0.033 

xs/turns-squared, p.u. 

2.994 

2.339 

2.996 

1.163 

1.611 

0.920 

aagnet overhang 

0.001 

0.298 

0.313 

0.326 

0.331 

0.329 

air gap flux density 

0.619 

0.687 

0.642 

0.722 

0.671 

0.700 

overall length 

3.730 

2.548 

4.345 

3.473 

4.160 

3.786 

overall disaster 

1.047 

1.078 

0.939 

1.450 

1.273 

1.390 
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Table 12. Characteristics of 2400 rpm magnet machines 


nutbtr of pol* pairs 

1 

2 

3 

4 

5 

6 

pomr, hp 

23775 

25775 

25/75 

25775 

25773 

25775 

efficiency factor 

1C0000 

100000 

100000 

100000 

IBOOOO 

180000 

voIum factor 

1286.1 

1206.1 

1206.1 

12B6.1 

1286.1 

1286.1 

lagnit factor 

25 

23 

25 

25 

25 

23 

shaft rps 

2400 

2400 

2400 

2400 

2400 

2400 

itator currant dimity 

12030000 

12000000 

12000000 

10202760 

11959030 

11997423 

synchronous infancy 

40 

B0 

120 

160 

200 

240 

rotor radius 

0.3584 

0.0005 

0.4010 

0.4408 

0.3361 

0.4760 

gap diHfisicin 

0.0371 

0.1883 

0.0136 

0.0008 

0.0032 

0.0165 

back iron dipth 

0.1640 

0.0420 

0.0625 

0.0504 

0.0290 

0.0379 

stator slot dipth 

0.0687 

0.0249 

0.0341 

0.0325 

0.0234 

0.0300 

legnet radial diamsion 

0.0390 

0.1970 

0.0143 

0.0092 

0.0033 

0.0173 

stator slot factor 

0.340 

0.750 

0.750 

0.730 

0.749 

0.678 

rotor slot factor 

0.378 

0.370 

0.378 

0.378 

0.378 

0.37B 

inyelopi voluse 

4.497 

1.342 

2.727 

3.157 

2.349 

3.438 

envelops Height 

26540.1 

8797.0 

12335.3 

13310.2 

13660.4 

11877.1 

hysteresis loss 

34020.0 

18160.1 

46136.0 

6B103.8 

86724.5 

9147P.9 

eddy current loss 

0130.4 

B677.9 

33069.4 

650B7.5 

103604.3 

131141.0 

stator copper loss 

234324.1 

217110.2 

127316.4 

90617.2 

112885.8 

91721.4 

full load efficiency 

0.986 

0.907 

0.989 

0.9B8 

0.984 

0.984 

active length 

1.536 

4.392 

1.291 

1.387 

3.061 

1.129 

eagnet weight 

384.207 

10.204 

167.444 

128.106 

73.723 

217.676 

lagnet voluie 

0.052 

0.001 

0.023 

0.017 

0.010 

0.029 

xs/turns-squared, p.u. 

2.417 

0.202 

1.973 

1.965 

2.879 

0.751 

■agnet overhang 

0.010 

0.330 

0.331 

0.315 

0.338 

0.324 

air gap flut density 

0.617 

0.635 

0.676 

0.672 

0.640 

0.692 

overall length 

3.270 

5.938 

3.010 

3.222 

4.431 

3.168 

overall diantir 

1.260 

0.511 

1.024 

1.065 

0.783 

1.121 
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Table 13. Characteristics of 3000 rpm magnet machines 


nuabcr of pole pairs 

1 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

value* factor 

1284.1 

1284.1 

1284.1 

1284.1 

1286.1 

1284.1 

•agnet factor 

25 

25 

25 

25 

25 

25 

shaft rpa 

3000 

3000 

3000 

3000 

3000 

3000 

stator current density 

10639824 9557320. 

12000000 

11999547 12000000 11400442 

synchronous frequency 

50 

too 

150 

200 

250 

300 

rotor radius 

0.2480 

0.2394 

0.2239 

0.0003 

0.3116 

0.4837 

gap dimension 

0.0227 

0.0075 

0.0035 

0.1403 

0.0027 

0.0375 

back iron depth 

0.1121 

0.0510 

0.0317 

0.0175 

0.0267 

0.0428 

stator slot depth 

0.0405 

0.0349 

0.0235 

0.0074 

0.0224 

0.0312 

•agnet 'adial diaension 

0.0238 

0.0079 

0.0037 

0.1483 

0.0028 

0.0394 

stator siot factor 

0.750 

0.750 

0.750 

0.750 

0.747 

0.681 

rotor slot factor 

0.378 

0.378 

0.378 

0.378 

0.378 

0.378 

envelope voluee 

2.213 

1.834 

1.771 

2.143 

1.950 

3.643 

envelope weight 

15051.7 

12394.7 

12381.0 

B254.8 

11507.2 

9943.3 

hysteresis loss 

22492.8 

34848.2 

57339.4 

4528*.3 

91118.4 

95420.6 

eddy current loss 

4777.4 

22010.1 

51374.8 

54098.3 134046.7 171347.8 

stator capper loss 

157512.4 108787.2 132084.7 

170746,7 

97480.0 

74849.0 

full load efficiency 

0.990 

0.991 

0.988 

0.984 

0.983 

0.9B3 

active length 

2.395 

3.719 

5.490 

16.723 

3.005 

0.732 

■agnet Height 

250.152 

123.493 

91.981 

13.685 

55.765 

330.030 

eagnet voluee 

0.034 

0.017 

0.012 

0.002 

0.008 

0.046 

xs/turns-squared, o.u. 

2.934 

2.993 

2.997 

0.043 

2.991 

0.351 

•agnet overhang 

0.007 

0.004 

0.324 

0.325 

0.311 

0.336 

air gap flux density 

0.41B 

0.4V? 

0.427 

0.423 

0.637 

0.734 

overall length 

3.574 

4.738 

4.415 

18.039 

4.274 

2.975 

overall diaeeter 

0.B47 

0.449 

0.545 

0.371 

0.727 

1.191 







Table 14. Characteristics of 3600 rpm magnet machines 


nuaber or pole pairs 

1 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

160000 

180000 

180000 

180000 

180000 

180000 

volute factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

aagnet factor 

25 

25 

25 

25 

25 

23 

shaft rpt 

3600 

3600 

3600 

3600 

3600 

3600 

stator current density 

12000000 12000000 

11966619 11677983 

11998747 11975558 

synchronous frequency 

60 

120 

180 

240 

300 

360 

rotor radius 

0.3044 

0.2596 

0.3006 

0.2526 

0.3700 

0.4536 

gap dieension 

0.0228 

0.0088 

0.0121 

0.0025 

0.0063 

0.0173 

hack iron depth 

0.1454 

0.0574 

0.0162 

0.0268 

0.0332 

0.0371 

stator slot depth 

0.0332 

0.0442 

0.0411 

0.0208 

0.0223 

0.0240 

■agnet radial ditension 

0.0240 

0.0092 

0.0127 

0.0026 

0.0067 

0.0182 

stator slot factor 

0.750 

0.553 

0.750 

0.735 

0.750 

0.749 

rotor slot factor 

0.378 

0.378 

0.37B 

0.378 

0.378 

0.378 

envelope volute 

2.475 

1.602 

2.427 

1.694 

2.063 

2.801 

envelope weignt 

15468.5 

9B53.6 

7072.7 

11238.1 

9644.1 

8977.4 

hysteresis loss 

29048.1 

36290.0 

37383.1 

84962.0 

92749.9 104238.5 

eddy current loss 

10410.6 

26012.0 

40193.3 121798.6 

166203.6 224148.8 

stator capper loss 

163951.5 114166.3 

104037.5 

95474.6 

73859.8 

66008.1 

full load efficiency 

0.990 

0.991 

0.991 

0.985 

0.983 

0.980 

active length 

1.394 

2.276 

1.279 

4.322 

1.669 

0.907 

signet weight 

221.175 

113.409 

104.473 

58.252 

91.121 

175.220 

tagnet volute 

0,030 

0.015 

0.014 

0.008 

0.012 

0.024 

xs/turns-squared, p.u. 

2.960 

2.977 

2.067 

2.981 

1.490 

0.593 

tagnet overhang 

0.336 

0.324 

0.327 

0.333 

0.340 

0.339 

air gap flux density 

0.664 

0.641 

0.663 

0.630 

0.661 

0.798 

overall length 

2.799 

3.386 

2.581 

5.352 

3.202 

2.863 

uverall diaaeter 

1.012 

0.740 

0.800 

0.605 

0.864 

1.064 





Table 15. Characteristics of 7200 rpm magnet machines 


nutber of pole pairs 
power, hp 
efficiency factor 
volute factor 
aagnet factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap ditension 
back iron depth 
stator slot depth 
aagnet radial ditension 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
tagnet weight 
aagnet volute 
xs/turns-squared, p.u. 
aagnet overhang 
air gap flux density 
overall length 
overall riiaeeter 


1 

2 

3 

25775 

25775 

25775 

180500 

180000 

180000 

1288.1 

1288.1 

1288.1 

25 

25 

25 

7200 

7200 

7200 

12000000 12000000 8313050. 

120 

240 

360 

0.2266 

0.1938 

0.2622 

0.0148 

0.0059 

0.0095 

0.1045 

0.0423 

0.0392 

0.0285 

0.0284 

0.0315 

0.0153 

0.0062 

0.0099 

0.749 

0.750 

0.750 

0.378 

0.3/8 

0.37B 

1.247 

0.811 

1.113 

8188.9 

5319.9 

6052.8 

30513.8 

3B176.1 

65438.7 

21871.7 

54727.9 

140716.0 

90777.3 

71500.7 

38515.8 

0.993 

0.992 

0.987 

1.552 

2.384 

1.621 

110.278 

57.403 

87.115 

0.015 

0.008 

0.012 

2.948 

2.884 

1.258 

0.322 

0.331 

0.307 

0.848 

0.635 

0-647 

2.578 

3.207 

2.747 

0.748 

0.541 

0.685 


4 

5 

6 

25775 

25775 

25775 

180000 

130000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

7200 

7200 

7200 

12000000 12000000 12000000 

480 

600 

720 

0.2627 

0.2467 

0.2418 

0.0041 

0.0020 

0.0022 

0.0289 

0.0209 

0.0171 

0.0244 

0.0114 

0.0104 

0.0042 

0.0021 

0.0023 

C.750 

0.748 

0.750 

0.378 

0.378 

0.378 

0,93! 

1.356 

1.398 

4898.9 

8797.5 

9221.7 

72593.1 

177996.0 218032.2 

208133.7 620002.2 937689.5 

52482.8 

49628.2 

49315.1 

0.983 

0.958 

0.941 

1.542 

3.965 

4.508 

36.486 

42.995 

50.295 

0.005 

0.006 

0.007 

2.246 

1.608 

1.142 

0-328 

0.338 

0.319 

0.650 

0.631 

0.629 

2.627 

4.968 

5.493 

0.640 

0.562 

0.543 


653 


Table 16. Characteristics of 180 rpm magnet machines 


nutber of pole pairs 

1 

2 

3 

4 

5 

power, hp 

?5775 

25775 

25775 

25775 

25775 

efficiency factor 

leoooo 

180000 

180000 

180000 

180000 

volute factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

•agnet factor 

25 

25 

25 

25 

25 

shaft rpi 

180 

ISO 

180 

180 

180 

stator current density 

7.80E+06 

1.18E+07 

1.20E+07 

1.09E+07 

1.20E+07 

synchronous frequency 

3 

6 

9 

12 

15 

rotor radius 

0.4399 

0.1787 

0.9745 

0.7389 

0.8993 

gap dimension 

0.0413 

0.0036 

0.0201 

0.0166 

0.0214 

back iron depth 

0.1994 

0.0377 

0.1491 

0.0779 

0.0865 

stator slot depth 

0.0932 

0.0282 

0.0327 

0.0476 

0.0630 

aagnet radial dieension 

0,0433 

0.003B 

0.0212 

0.0174 

0.0225 

stator slot factor 

0.463 

0.531 

0.750 

0.750 

0.740 

rotor slot factor 

0.378 

0.378 

0.378 

0.378 

0.378 

envelope volute 

29.375 

36.694 

34.709 

20.521 

22.314 

envelope Height 

211948.9 284972.1 

156312.7 107071.7 

77934.3 

hysteresis loss 

20311.6 

53970.7 

47089.7 

41032.7 

36926.0 

eddy current loss 

364.0 

1934.3 

2531.5 

2941.1 

3308.5 

stator copper loss 

7.90EV05 2.54E+06 6.97E+05 7.01E+05 

7.11E+05 

full load efficiency 

0.960 

0.881 

0.963 

0.963 

0.962 

active length 

12.100 

171.774 

3.195 

4.560 

1.866 

•agnet Height 

4072.4 

2079.6 

1941.0 

1048.8 

1087.7 

■agnet volute 

0.550 

0.281 

0.262 

0.142 

0.147 

xs/turns-squared, p.u. 

2,996 

2.990 

2.981 

2.983 

2 870 

■agnet overhang 

0.003 

0.034 

0.332 

0.010 

0.339 

air nap flux density 

0.619 

0.619 

0.674 

0.618 

0,704 

overall length 

14.198 

172.518 

7.258 

7.651 

5.639 

overall diateter 

1.548 

0.496 

2.353 

1.762 

2.140 
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Table 17. Characteristics of more 180 rpm magnet machines 


miaber of pale pairs 

6 

7 

8 

9 

10 

power, lip 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

volute factor 

1298.1 

1286.1 

1286.1 

1286.1 

1286.1 

tagnet factor 

25 

25 

25 

25 

25 

shaft rpt 

180 

180 

180 

180 

180 

stator current density 

6.0IE+06 

1.20E+07 

1.20E+07 9.42E+06 

1.20E+07 

synchronous frequency 

18 

21 

24 

27 

30 

rotor radius 

0.6047 

1.2439 

0.9696 

0.9792 

0.9526 

gap ditension 

0.013? 

0.1161 

0.0089 

0.0031 

0.0079 

back iron depth 

0.0426 

0.1117 

0.0546 

0.0469 

0,0430 

stator slot depth 

0.0337 

0.0627 

0.031? 

0.0255 

0.0360 

tagnet radial intension 

0.0146 

0.1219 

0.0093 

0.0033 

0.0083 

stator slot factor 

0.536 

0.750 

0.74B 

0.575 

0.74d 

rotor slot factor 

0.378 

0.378 

0.378 

0.378 

0.376 

envelope volute 

44.578 

53.524 

2V.253 

43.668 

26.292 

envelope weight 

313139.6 

98601.8 

120336.5 234296.3 102286.2 

hysteresis loss 

186214.5 

6751B.3 

96060.7 214344.0 101073.4 

eddy current loss 

20021.3 

8469.3 

13770.9 

34568.5 

18111.9 

stator copper loss 

3.91E+05 

7.15E+05 

5.34E+05 3.91E+05 

5.31E'05 

full load efficiency 

0.970 

0.960 

0.968 

0.968 

0.967 

active length 

24.182 

0.651 

3.523 

7.417 

3.167 

tagnet weight 

3800.1 

3431.9 

932.5 

704.3 

714.1 

tagnet volute 

0.514 

0.464 

0.126 

0.095 

0.097 

ns/turns-squared, p.u. 

0.546 

0.525 

2.350 

2.990 

2.391 

tagnet overhang 

0.010 

0.32? 

0.336 

0.340 

0.325 

air gap flux density 

0.61? 

0.859 

0.66? 

0.644 

0.671 

overall length 

26.715 

6.579 

7.474 

11.360 

7.042 

overall diateter 

1.390 

3.049 

2.128 

2.109 

2.079 
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$150,000, rendering the 180 rpm machines less economic to 
build. A machine aucn as these may cost between four and 
eight million dollars. Only those machines with magnet 
costa that are reasonable with respect to the other material 
costs should be candidates for design. 

In all the permanent magnet machines, stator currant 
density went to the maximum. The overall length and overall 
diameter of the more reasonable machines would allow them to 
fit in machinery spaces aboard a ship. 

4.5. Verification 

No high-power permanent magnet machines were discovered 
during the search to find a benchmark. Because of the high 
material cost and the competition afforded by synchronous 
and induction machines, it seems none have been built. 
Several paper studies were found [13, 19, 20, 21, and 22], 
and the paramete s resulting from this computer modeling 
seem to agree with them. The machine size is what was ex¬ 
pected, given the lower air gap flux density. The ef¬ 
ficiency was higher than the synchronous and induction 
machines. All-in-all, this modeling gave good machines. 




Chapter Five. 


Induction machines 


The stator of an induction machine is the same as those 
of synchronous and permanent magnet machines. The rotor is 
significantly different. There is no independent, mechanism 
to produce a rotor flux wave. The rotor winding is shorted, 
whether it is wound or cast, so that as the stator flux wave 
passes over the rotor, currents are induoed in the winding. 
These currents produce only a small reaction flux, but it 
still tends to align with the stator flux wave. When at 
operating speeds, the rotor speed is a bit slower than the 
stator flux wave speed, and the difference in speeds is 
called slip. Typically, slip is a few percent of the stator 
frequency. The rotor currents are at slip frequency. If 
the rotor and stator speeds were the same, slip would be 
zero, there would be no tendency to align and torque would 
be zero. Then, the rotor would Jag behind the stator until 
current was Induced in the rotor winding by the passing 
stator flux wave and torque was again produced. 

If solid bars are used as the rotor winding, they are 
shorted at the ends of the rotor by end rings, to form what 
is called a "squirrel cage" rotor. If actual turns are 
used, the wiuding may be shorted through external resis¬ 
tances to affect the starting and torque-slip characteris¬ 
tics of the machine. Fitzgerald et al [23] and Alger [24] 
discuss induction machine characteristics in some detail. 

5.1. Assumptions 

A squirrel cage rotor was assumed for these machines. 
Copper was designated as the material for the rotor bars. 
However, these machines will be fed from a frequency 
changer, so only one layer of bars was ur.ed ana the effects 
of magnetic diffusion ignored in the analysis (see 
Appendix D for a derivation of the components of an induc¬ 
tion machine equivalent circuit). The number of rotor bars 
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was arbitrarily sat at 71. This quantity should not causa 
undesirable harmonics, as it will not be an integral mul¬ 
tiple of the number of poles or stator slots in any machine. 
The number and width of the rotor bars were inextricably 
entwined and could not be separated in the analysis. 

5.2. Machine description 

An attempt was made to design the same rpu and pole- 
pair machines as was performed for the synchronous and per¬ 
manent magnet models, but problems in limiting rotor current 
density allowed only a few of the machines to be designed. 
For example, no 7200 rpm machines were designed and 180 rpm 
machines could only be designed with up to twelve poles. 

All of the induction machines are listed in the tables 
starting on page . Only medium confidence should be placed 
in the induction machine designs, as there were some conver¬ 
gence difficulties in slip. (Slip is not listed in the 
tables for that reason.) 

5.2.1. Efficiency 

The higher rpm machines showed a slight increase in ef¬ 
ficiency as rpm increased. There was much movement around 
the average value of 97.5%. The movement decreased as rpm 
Increased. With only six machines, it is hard to detect a 
trend in 180 rpm machine efficiency. Apparently, efficiency 
did increase with the number of pole-pairs, with all ef¬ 
ficiencies below 90%. Developmental studies for this thesis 
showed that off-design-point efficiencies for the 180 rpm 
maohines were sometimes below 70% for the endurance speed 
condition. 
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Figure 14. High rpm induction machine efficiency 
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Fitfure 16. 180 rpm induction machine efficiency 
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Figure 16. High rpm induction machine volume 
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Figure 18. High rpm induotion machine weight 
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5.2.2. Weight and volume 


The best of the high rpm machines rivaled the 
synchronous machines in weight and volume. The worst were 
very bad. Weight and volume decreased with an increase in 
the number of poles, but not necossarily with the increases 
in rpm. For the 180 rpm machines, both weight and volume 
decreased dramatically as the number of poles went from two 
to six, with nmoh lower decreases after that. The 180 rpm 
machines were uncompetitive in the synthesis process. 

5.3. Machine characteristics 

The previously mentioned rotor current density dif¬ 
ficulty showed in the rotor slot factor, which was at the 
limit of 0.75 for almost every motor. The stator slot fac¬ 
tor gradually grew with the increase in poles, arriving at 
0.75. The length and diameter of both the 180 and higher 
rpm machines is such that they would fit in machinery 
spaces. 

5.4. Verification 

Induction machines were expected to be close to 
synchronous machines in volume, weight, and efficiency. 

They were, and this comparison served as the verification 
for the induction machine model. Because the confidence 
level in the designed machines is only medium, more worh 
would be needed to verify that these machines would have the 
advertised properties if built. 
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Table 18. Characteristics of 1800 rpm induction machines 


nueber of pole pairs 

1 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

volute factor 

1288.1 

1286.1 

1286.1 

1286.1 

1286.1 

12B6.1 

shaft rpt 

1800 

1800 

1B00 

1800 

1800 

1800 

pritary awp-turns 

A.10E+05 7.74E+05 

2.44E+05 1.91E+05 

1.86E+03 

1.79E+05 

synchronous frequency 

30 

60 

90 

120 

150 

180 

rotor radius 

0.8236 

1.0197 

0.4638 

0.4352 

0.4060 

0.3945 

pap ditension 

0.0262 

0.0039 

0.0036 

0.0029 

0.004? 

0.0051 

back iron depth 

0.1B94 

0.0173 

0.1073 

0.0093 

0.0405 

0.0460 

stator slot depth 

2.1467 

0.5423 

0.3375 

0.3490 

0.1902 

0.13B8 

rotor slot depth 

0.2051 

0.3399 

0.1546 

0.1451 

0,1353 

0.1176 

stator slot factor 

0.285 

0.641 

0.355 

0.640 

0.699 

0.74? 

rotor slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

0.748 

envelope voluee 

148.863 

36.240 

7.558 

5.B80 

3.760 

3.164 

envelope weight 

376603.7 

80669.2 

29447.0 

21499.7 

15333.1 

13395.9 

hysteresis loss 

180851.2 

19122.1 

68838.7 

39013.4 

42456.7 

48214.8 

eddy current loss 

32425.7 

6853.2 

37006.8 

27964.1 

38040.3 

51839.3 

stator copper loss 

1.03E+06 1.45E+06 

498668.3 165391.9 

264414.2 

320360.8 

full load efficiency 

0.93? 

0.928 

0.969 

0.987 

0.981 

0.977 

active length 

0.191 

0.090 

0.75? 

0.930 

1.000 

1.081 

rotor copper loss 

8427.7 

8543.0 

19627.2 

22347.5 

30471.4 

37148.9 

eaxieue torque 

185470 

185470 

1B5470 

1B5470 

1B5470 

185470 

terminal voUs/turn 

234.31 

45.07 

266.45 

233.56 

193.76 

189.53 

air gap volts/turn 

63.2? 

32.92 

126.58 

145.45 

145.9B 

153.34 

Rl/turns-squared 

?.23E-07 8.09E-07 

2.80E-06 1.50E-06 

2.54E-06 

3.32E-06 

m/turns-squared 

3.70E-04 3,B?E-05 

9.60E-04 9.53E-04 

6.B2E-04 

6.17E-04 

Ke/turns-squared 

2.51E-03 4.41E-03 

1.23E-02 1.28E-02 

6.23E-03 

5.15E-03 

12/turns-squared 

2.03E-05 3.66E-05 

3.76E-04 6.52E-04 

7.79E-04 

9.08E-04 

R2/turns-squared 

4.50E-08 9.97E-0? 

3.89E-07 5.37E-07 

6.62E-07 

B.49E-07 

overall length 

3.990 

4.184 

2.629 

2.682 

2.644 

2.680 

overall ditKter 

6.572 

3.166 

1.824 

1.393 

1.283 

1.169 
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Table 19. Characteristics of 2400 rpm induction machines 


flutter of pole pairs 

2 

3 

4 

5 

4 

potter, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

ICOOOO 

180000 

volute factor 

1286.1 

1284.1 

1284.1 

1284.1 

1284.1 

shaft rpa 

2400 

2400 

2400 

2400 

2400 

priaary aap-turns 

5.94£*05 

3./8E+05 

2.45E+05 

1.78E+05 

1.93E+05 

synchronous frequency 

80 

120 

140 

200 

240 

rotor radius 

0.7541 

0.7958 

0.4803 

0.4124 

0.3841 

gap ditens'.on 

0.0049 

0.0020 

0.0094 

0.0087 

0.0081 

back iron depth 

0.0952 

0.0013 

0.0834 

0.0575 

0.0418 

stator slot depth 

0.4149 

0.5475 

0.2755 

0.2341 

0.1541 

rotor slot depth 

0.2514 

0.1979 

0.140! 

0.1375 

0.1197 

stator slot factor 

0.454 

0.477 

0.417 

0.413 

0.750 

rotor slot factor 

0.750 

0.750 

0.439 

0.548 

0.750 

envelope volute 

23.889 

21.180 

4.201 

4.554 

2.942 

envelope weight 

40524.1 

30880.0 

19931.3 

17243.1 

11459.0 

hysteresis loss 

78227.5 

24111.4 

80424.4 

84278.8 

54014.3 

eady current loss 

37381.4 

18714.2 

74844.2 103071.8 

77433.0 

stator copper loss 

1.10E+06 334484.0 374335.4 304049.2 282547.9 

full load efficiency 

0.940 

0.980 

0.972 

0.974 

0.978 

active length 

0.142 

0.189 

0.533 

0.909 

0.852 

rotor copper loss 

8819.4 

7424.0 

17135.5 

21589.8 

27954.3 

■axltut torque 

139102.5 139102.5 

139102.5 139102.5 

139102.5 

terminal volts/turn 

103.23 

113.49 

229.40 

344.92 

212.59 

air gap volts/turn 

51.27 

72,14 

122.49 

179.75 

154.78 

Rl/turns-squared 

1.04E-04 7.B2E-07 2.0BE-04 3.23E-04 2.52E-04 

U/turns-squared 

1.50E-04 2.32E-04 7.89E-04 1.67E-03 7.40E-04 

U/turns-squared 

3.8BE-03 1.25E-02 3.37E-03 4.31E-03 3.35E-03 

12/turns-squared 

6.22E-05 1.54E-04 4.43E-04 1.07E-03 9.34E-04 

R2/turns-squared 

2.84E-08 4.44E-08 2.97E-07 7.49E-07 4.73E-07 

overall length 

3.184 

3.380 

2,492 

2.594 

2.420 

overall diateter 

2.944 

2.493 

1.497 

1.724* 

i m 
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Table 20. Characteristics of 3000 rpm induction machines 


nuaber of pole pairs 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

l&OOOO 

1BOOOO 

180000 

180000 

1BOOOO 

voluw factor 

!784.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpi 

3008 

3000 

3000 

3000 

3000 

priaary aap-turns 

3.I0E+05 

3.62E+05 2.80E+05 : 

1.74E+05 

1.37E+05 

synchronous frequency 

100 

150 

200 

250 

300 

rotor radiu; 

0.51f8 

0.6366 

0.5580 

0.3660 

0.3325 

qap dimension 

0.0138 

0.0036 

0.0070 

0.0070 

0.0039 

back iron depth 

0.1789 

0.0134 

0.0178 

0.0503 

0.0388 

stator slot depth 

0.8926 

0.3785 

0.3668 

0.2274 

0.1288 

rotor slot depth 

0.1733 

0.2122 

0.1860 

0.1220 

0.1093 

stator slot factor 

0.340 

0.607 

0.545 

0.476 

0.721 

rotor slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

envelope voluie 

22.412 

10.122 

8.000 

3.379 

2.045 

envelope weight 

72906.0 

19541.1 

15961.5 

13000.9 

9016.0 

hysteresis loss 

185978.6 

25156.1 

44235.6 

74610.2 

52848.0 

eddy current loss 

111088.5 

22539.3 

52845.5 111415.1 

94701.2 

stator copper loss 

322718.0 359218.B 221369.1 269959.2 226785.8 

full load efficiency 

0.968 

0.979 

0.993 

0-974 

0.9B0 

active length 

0.383 

0.189 

0.307 

0.818 

1.006 

rotor copper loss 

8574.1 

9089.1 

9006.9 

18292.9 

25568.1 

aamtui torque 

111282 

111282 

111282 

111282 

111282 

terainal volts/turn 

421.67 

101.44 

174.86 

361.45 

251.07 

air gap volts/turn 

119.20 

72.02 

102.80 

179.30 

200.52 

Rl/turns-squared 

1.I2E-06 

9.14E-07 9.39E-07 2 

I.91E-04 4.01E-06 

Xl/turns-squared 

1.31E-03 

1.96E-04 5.04E-04 1.7BE-03 

1.09E-03 

Ii/turns-squared 

4.#76-03 6.93E-03 3.82E-03 5.33E-03 9.00E-03 

X2/turns-squared 

1.82E-04 1.72E-C4 3.25E-04 9.97E-04 

1.54E-03 

R2/turns-squared 

1.60E-07 5.15E-08 1 

.08E-07 4.66E-07 

1.01E-O6 

overall length 

2.517 

2.750 

2.567 

2.310 

2.352 

overall diaeeter 

3.210 

2.064 

1,899 

1.301 

1.008 
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Table 21. Characteristics of 3600 rpm induction machines 


nuiber of pole pairs 

1 

2 

3 

4 

5 

6 

poeer, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency fjctor 

1SOOOO 

180000 

180000 

180000 

180000 

180000 

volute factor 

128b.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpt 

3800 

3600 

3600 

3600 

3600 

3600 

pritary aap-turns 

2.72E+05 1.70E+05 2.99E+05 

2.95E+05 2.05E+05 I.23E+05 

synchronous frequency 

60 

120 

180 

240 

300 

360 

rotor radius 

0.4964 

0.4350 

0.5150 

0.5305 

0.4336 

0.3144 

gap ditension 

0.0061 

0.0027 

0.0333 

0.008E 

0.0079 

0.0030 

back iron depth 

0.3428 

0.0133 

0.0238 

0.0238 

0.0258 

0.0367 

stator slot depth 

1.3891 

1.0052 

0.3606 

0.3236 

0.2792 

0.1181 

rotor slot depth 

0.1654 

0.1308 

0.1717 

0.1768 

0.1446 

0.1045 

stator slot factor 

0.290 

0.338 

0.585 

0.688 

0.553 

0.750 

rotor slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

0.749 

envelope volute 

41.461 

17.960 

6.532 

6.536 

4.331 

1.722 

envelope Might 

164461.5 

54127.3 

14922.3 

14074.1 

11549.4 

7593.8 

hysteresis loss 

261300.3 145818.0 

34381.6 

42814.0 

62801.7 

52473.1 

eddy current loss 

93647.7 104519.7 

36966.2 

61376.6 

112537.8 112835.4 

stator copper loss 

341655.1 

107764.3 2B8B57.8 

210564.6 173450.1 

192320.6 

full load efficiency 

0.965 

0.981 

0.981 

0.983 

0.982 

0.981 

active length 

0.393 

0.700 

0.247 

0.249 

0.400 

0.964 

rotor copper loss 

7562.3 

7220.5 

10787.1 

9954.8 

11162.4 

23436.6 

■axiiue torque 

92735 

92735 

92735 

92735 

92735 

92735 

tereinal volts/turn 

475.12 

702.61 

143.82 

142.32 

275.51 

261.98 

air gap volts/turn 

140.25 

21B.89 

91.39 

94.86 

149.83 

218.01 

Rl/turns-squared 

1.54E-06 1.25E-06 

1.07E-06 

B.08E-07 

1.37E-06 4.25E-06 

ll/turns-squared 

1.67E-03 3.9JE-03 3.70E-04 3.59E-04 

1.13E-03 1.18E-03 

le/turns-squared 

2.39E-02 4.19E-02 9.66E-03 2.80E-03 3.94E-03 I.26E-02 

12/turns-squared 

1.20E-04 4.65E-04 2.54E-04 3.05E-04 7.24E-04 

1.84E-03 

R2/turns-squared 

1.98E-07 4.57E-07 

1.03E-07 9.68E-08 

2.79E-07 

1.07E-06 

overall length 

2.403 

2.451 

2.320 

2.406 

2.247 

2.234 

overall diaeeter 

4.469 

2.913 

1.605 

1.773 

1.494 

0.944 
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Table 22 


Characteristics of 180 rpm induction machines 


nuaber of pole pairs 

1 

2 

3 

4 

5 

6 

poaer, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180U00 

180000 

160000 

180000 

160000 

voluae factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

180 

ieo 

idO 

180 

180 

180 

priaary aaps-turns 

1054077 564571.5 569333.6 

510008.9 

512534.4 

324264.5 

synchronous frequency 

3 

6 

9 

12 

15 

18 

rotor radius 

0.9956 

0.7236 

0.6896 

0.6616 

0.7411 

9.6962 

gap diaension 

0.0025 

0.0060 

0.0195 

0.0334 

0.0235 

0.0054 

back iron depth 

0.0227 

0 . 0240 

0.1568 

0.1158 

0.1038 

0.0812 

stator slot depth 

0.6714 

0.7116 

0.2536 

0.2231 

0.2200 

0.1466 

rotor slot depth 

0.3319 

0.2412 

0.2299 

0.2039 

0.1845 

0.1501 

stator slot factor 

0.590 

0.563 

0.734 

0.750 

0.694 

0.718 

rotor slot factor 

0.750 

0.750 

0.750 

0.683 

0.696 

0.701 

envelope voluae 

47.668 

36.837 

20.316 

19.450 

21.357 

17.600 

envelope Might 

217524.1 

168344.2 

98690.0 

'*4554.9 

93084.9 

87558.9 

hysteresis loss 

4551.5 

16194.6 

19356.0 

24294.6 

31609.4 

36548.6 

eddy current loss 

81.6 

580.4 

1040.6 

1741.4 

2832.1 

3929.6 

stator copper loss 

4.41E+06 '.96E+06 3.05E+06 

3.15E+06 

2.59E+06 

2.09E+06 

full load efficiency 

0.808 

0.901 

0.854 

0.850 

0.872 

0.891 

active length 

0.825 

2.047 

1.854 

2.486 

2.159 

3.090 

rotor cop-er loss 

1.0BE+05 1.29E+05 2.0BE+05 

2.23E+05 

1.94E+05 

2.14EM15 

eaxiaui torque 

1854700 

1854700 

1854700 

1854700 

1854700 

IB547r 

terainal volts/turn 

61.849 

104.401 

58. 731 

76.530 

76.097 

87.418 

air gap volts/turn 

29.510 

53.243 

45.967 

59.124 

57.510 

77.345 

Rl/turns-squared 

1.33E-06 2.05E-06 3 

1.14E-06 

4.04E-06 

3.28E-06 

6.61E-06 

Xl/turns-squared 

5.0BE-05 1.58E-04 6.01E-05 

9.01E-05 

9.33E-05 

1.12E-04 

la/turns-squared 

1.21E-02 4.61E-03 8.17E-04 

4.60E-04 

5.08E-04 

2.49E-03 

12/turns-squared 

' 1.89E-05 6.91E-05 7.91E-05 

1.37E-04 

1.27E-04 

2.40E-04 

R2/turns-squared 

1.16E-07 4.44E-07 4.33E-07 

7.40E-07 

6.20E-07 

1.15E-06 

overall length 

4.817 

4.965 

4.691 

5.266 

5.217 

5.897 

overall diaaeter 

3.384 

2.930 

2,239 

2.068 

2.177 

1.859 
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Chapter Six. Nominal ship design 


6.1. Technology sensitivity analysis 

Technology sensitivity analyses, such as this thesis, 
must be able to quantitatively compare similar technologies. 
The pertinent differences must be made apparent through ap¬ 
propriate analysis. Inherent in the analysis must be the 
consideration of the global system complexity. Naval ships 
are extremely complex and the effects of various tech¬ 
nologies can be lost in the complexity. One methodology for 
technology characterization in naval ships has been proposed 
by Goddard [25]. This method has been followed to show the 
benefits of electric drive. 

6.2. ASSET 

The Advanced Surface Ship Evaluation Tool (ASSET) was 
developed over several years to be the U. S. Navy’s premier 
ship design computer program. It has its roots in HANDE, a 
hydrofoil design program developed in the 1970’s. The ship 
design apiral is traversed in an iterative fashion until 
convergence on a number of parameters is achieved. Boeing 
Computer Services is the contractor for ASSET, under the su¬ 
pervision of the David Taylor Naval Ship Research and 
Development Center, at Carderock, Maryland. 1 # It is divided 
into large modules by ship type. These modules include 
Monohulls, Hydrofoils, and SWATH ships. 

The geometry of a particular ship is input to ASSET. 

The following were used as characteristic ship traitsi full 
load displacement, certain Ship Work Breakdown Structure 
(SWBS) weight groups, 1 7 endurance fuel load at 20 knots, 

16. Greenwood, R.W., and Fuller, A.L., "Development of a 
Common Tool for Ship Design and Technology Evaluation," 

Proc. SNAME New England Section Marine Computers 1986. 

17. Soe Appendix E for more information on SWBS. 
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draft, maximum and sustained speed for 51,550 installed 
horsepower, and transverse GM. From the changes in these 
parameters during the various computer runs, the effects of 
transmission choices were noted. 

The inputs to ASSET describe the ship that is being 
designed. The outputs cover the range of calculations pos¬ 
sible in structures, volume, space, machinery, propeller 
characteristics, resistance, powering, and weight. There 
are ASSET performance modules on coat, stability, hydros¬ 
tatics, seakeeping, manning and space but the usual syn¬ 
thesis output is of more use during a technology sensitivity 
analysis. 

The descriptions of several ships are contained in an 
ASSET data base. For a particular ship, a Current Model is 
maintained that holds all of the parameters to desoribe that 
particular ship. In ASSET Version 2.0, over 380 parameters 
are used for each ship description. User control over most 
of these parameters is possible, or control may be given to 
the executive program which will then "design" a ship sub¬ 
ject to whatever constraints the user desires. 

Some intricacies of ship design are not handled well by 
ASSET. For example, the program is not able to handle 
equipment re-arrangements easily, and almost all equipment- 
level volumes are approximated from studies of past ships. 18 
For this reason, some equipment-level weights and volumes 
need to be calculated off-line and input to the program 
through its weight adjustment facility, especially if ac¬ 
curacy in these areas is important to the study being per¬ 
formed. 

The baseline ship used in these studies has a full load 
displacement of 5485 LT, carries 272 crew members, is 425 
feet long, and has a primary mission of anti-submarine war- 

18. The Enhanced Machinery Module [27], in the process of 
being made available, will improve the situation dramati¬ 
cally. Some of the relationships from that module were used 
in calculating electric propulsion weights. 
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fare (ASW). It is armed for that purpose and has equipment 
in keeping with its size The baseline ship is described in 
more detail below. 

6.2.1. Margins 

A naval ship design has margins in weight, vertical 
center of gravity (KG), space, ship service electrical gen¬ 
eration, propulsion power, accommodations, and structural 
strength which allow for equipment, mission, and system 
growth over its projected thirty year life. Without these 
margins, the ship would be difficult to modify, because 
whatever might be added in these areas would have to be paid 
for by a removal. For example, if a 5G ton radar system 
were added, the original 40 ton radar system and 10 tons of 
fuel might be removed to leave the ship at its original 
weight. With margins, the 10 tons of fuel might not be 
removed. 

Margins are typically split into Acquisition and Serv¬ 
ice Life allowances. Acquisition margins recognize the faot 
that ship specifications change over the design cycle and 
during construction. For example, the fourth ship built may 
have a different weapons system than the first, with a dif¬ 
ferent electrical requirement. If the electrical generation 
plant had to be changed during construction to accommodate 
the new weapons system, the total cost might be prohibitive. 
If an Acquisition margin is built into the original design, 
this may not occur. A Service Life margin makes allowance 
for configuration changes over the life of the ship. 

The ASSET program uses margins when synthesizing a 
ship. The margins are under operator control. The margins 
suggested by Goddard and used in this analysis are listed in 
the table.i# 

19. Goddard, C.H., "A Methodology for-Technology Charac¬ 
terization and Evaluation for Naval Ships," S.M. HAME AND 
O.E. thesis, MIT, 1985, p. 31. 


85 




Table 23. Recommended technology assessment design margins 
for a monohull surface combatant 


Acquisition Service Life 


Weight* 

KG 

Space 

Electrical* 5 
Propulsion powers 

Accommodations 

Strength 


12.5* of SWBS 1-7 10* 

12.5* of KG 1-7 1.0 ft. 

0 (no excess volume) 0 

20* 20* 


10* total EHP prior to prelim body plan 
8* prior to self-propelled model tests 
Accom = 1.1 x ship manning at delivery 
2.24 KSI of marginal stress at delivery 
(Max primary stress for hull material) 


Notes: 

a. The service life weight margin applies only to naval 
architectural limits of the ship (reserve buoyancy, 
stability, structures), not to the final design weight. 

b. In sizing the electric plant, the calculated maximum 
electric load plus these design margins shall be met with 
one generator out of service. The remaining generators 
shall not be loaded in excess of 90*. Note that the service 
life margin is not applied to SWBS group 200 which would be 
expected to remain stable over the life of the ship. 

c. Performance requirements (sustained speed, endurance 
range) are met at delivery full load displacement. 


6.3. Phiiosophy of effort 

The nature of this technology characterization required 
that certain limits be imposed on the total effort. (If the 
Naval Sea Systems Command were to do this study, many people 
would simultaneously be employed to investigate every 
detail.) Some items were fixed, some were allowed to float 
with the design. 

The hardest item to handle is volume. There are very 
few ways of adjusting volume as easily as weight is ad¬ 
justed. One way is through the use of Marginal Volume Fac¬ 
tors, which equate a weight penalty with every increase in 
volume. (See Howell [26].} 

The differences among transmission systems appear 
primarily in the machinery spaces and fuel tanks of a ship. 
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ASSET can handle the tankage changes without aid. However, 
it does no equipment arrangement or space analysis inside 
the machinery spaces, leaving that for arrangement experts 
to do off-lino. It was decided to keep the machinery space 
volume the same in all ships of this study, no matter how 
that volume was divided into spaces. Changes in equipment 
volume will be noted in the analysis and left for the advan¬ 
tage of others in machinery space rearrangement. 

Moderate to high technical risk has been accepted in 
specifying equipment cooling. Current densities at the 
limits of cooling technology are used, assuming that liquid 
cooling of both the stator and rotor can be performed. 

Lower current densities would be required if such cooling 
were not possible, resulting in slightly larger, heavier, 
less efficient machines. This last statement was proven 
during the course of the thesis research, as the first, 
(chronologically) current densities used were two-thirds 
larger than those listed herein. 

Other risk areas include the use of advanced vacuum 
switchgear and the assumed efficiencies of reduction gears 
and power converters. These are low risk items; the 
technology is well understood and commercially available. 


Table 24. Ship design items held constant during analysis 


Endurance speed 
Endurance range* 

Machinery box volume 
Installed horsepower 

Payload weight, volume, and electrical 

Length 

Beam 

Ship electrical load 
Ship molded lines 
Manning 

Deckhouse and superstructure geometry 


20 knots 

5500 nautical miles 
109,670 ft* 

51,560 hp 
requirements 
425 ft. 

55 ft. 

2030 kW (24 hr. avg) 
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Table 25. Ship design items allowed to float during analysis 

Maximum and sustained speeds 
Endurance fuel load 
Ship full load displacement 
Ship draft 

Ship resistance and powering 

Ship arrangement outside machinery spaces 


Note: 

a. Endurance range was allowed to float as a comparison 
between two electrical transmission ships and the mechanical 
transmission baseline ship. It was held constant for the 
rearranged ship. 




Most of the ship synthesis has been left to ASSET 
though parameters from Goddard were used where possible. 
ASSET designs a reasonable, generic ship with good seakeep¬ 
ing characteristics. As shown above, some ship characteris¬ 
tics were frozen to ensure transmission comparisons were not 
performed with different ships. Leaving the ship synthesis 
to ASSET allowed concentration on the specifics of the 
propulsion plant. 

6.4. Baseline ship 

The baseline ship has a mechanical transmission, i.e., 
two power trains consisting of a gas turbine, clutch and 
coupling, reduction gear, shafting, and propeller. There 
are two machinery rooms, each containing one gas turbine. 

The gas turbine used as the model is the General Electric 
LM-2500, rated at 25,775 brake horsepower. This is a very 
common marine gas turbine. It and its predecessors are 
powering the latest classes of naval combatant, such as the 
DD963, FFG7, and DDG51. 

The locked-train double-reduction gears are reversible, 
allowing the use of fixed-pitch propellers. There is no 
mechanical cross-connect allowed between the shafts. Except 
for the power level, this is the gear system being employed 
on the DDG51 class. Gears of this sort are about one per¬ 
cent inefficient* 0 per reduction stage. Since these are 
double-reduction gears, an efficiency of 98% was used. 

An endurance speed of 20 knots has been specified. 

This is in keeping with standard fleet practices. An en¬ 
durance range of 5500 nautical miles permits ocean crossings 
without refueling. The lack of a cross-connect capability 
between the two shafts means at least two gas turbines will 


20. Inefficiency = 1 - efficiency. Information on stage in¬ 
efficiency is from a conversation with Mr. Samuel Shank, the 
author of the ASSET Enhanced Machinery Module [27], 
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be on-line during endurance cruising. This is inherently 
inefficient; some operators alleviate the inefficiency by 
declutching one shaft and free-wheeling that propeller, 
reaching desired speeds by loading up the operating gas tur¬ 
bine . 

A measure of initial static stability is the ratio 
GM/Beam. GM is the vertical distance between the center of 
gravity and the metacenter of the ship. Typical values for 
this ratio are 8-10%. A lower value (6.5%) has been ac¬ 
cepted for purposes of comparison with the variants. A 
large ship redesign effort would have been necessary to 
bring GM/B into a better range. 

The electrical generation plant consists of three gas 
nurbines, each driving a 2000 kW generator. The data used 
for the gas turbines was taken from the Detroit Allison 501R 
turbine-generator set used aboard the BD963 ships. 

Both the deckhouse and main hull are constructed of 
High Tensile Strength (yield stress = 50,000 psi) steel. 
Active stabilizing fins and a sonar dome are included in the 
design. The payload is listed in Table 26. A coarse layout 
of the machinery spaces is shown in Figure 20. 


Table 26. Payload for baseline and variant ships 

FFG7 Command and control suite 
Satellite, UHF, and HF communications 
SLQ32V3 electronic countermeasures 
NIXIE aooustic countermeasure 

SPS-49 two-dimensional air search and tracking radar 

SPS-55 surface search and traoking radar 

SQR-19 towed array surveillance system 

MK92 missile and gun fire control system 

Harpoon fire control system 

LAMPS III helicopters and support system 

JP-5 aviation fuel 

MK32 over-the-side torpedo system 

MK13 guided missile launching system 

MK75 76mm gun 

Close-In-Weapon-System 

Small arms 

Appropriate ammunition and reloads 
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6.5. Backfit ship 


If the mechanical transmission of the baseline ship 
were replaced with an electrical transmission without chang¬ 
ing the subdivision of the ship, it would be as if an older 
ship had been updated, or "backfitted," with new technology. 
This is the idea behind the first variant ship. 

Propulsion motors and generators were added into the 
propulsion plant, in the original spaces. Shafting still 
runs from the machinery spaces to the propellers. A reduc¬ 
tion gear is still necessary for higher rpm propulsion 
motors, but electrical cross-connect may improve the en¬ 
durance fuel efficiency. 

Some rearrangement within the machinery space is neces¬ 
sitated by tha backfit. The machinery spaces in the 
baseline ship are not long enough to contain the stack-up 
length of gas turbine, generator, motor, and reduction gear, 
without "folding" the power train. This may be accomplished 
by changing the design of the reduction gear or by placing 
the gas turbine and generator (which require a mechanical 
connection) side-by-side (or transversely) with the motor 
and reduction gear (another mechanical connection), using 
transmission line to electrically connect them. Since there 
are a variety of ways to rearrange the machinery box, and 
since the chosen method has no effect on the analysis, the 
rearrangement was not specified. 

The propulsion generators and motors may operate at 
different rotational speeds, and therefore different elec¬ 
trical frequencies, so power converters must be used between 
them. Power converters change the frequency of the power 
being transferred between the generators and motors, through 
the use of cycloconverters or thyristors. They add another 
inefficiency to the transmission. A reasonable estimate of 
the efficiency of an 18 MW power converter is 97%.2i 

21. Professor John Kassakian, MIT, private communication. 
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6.6. Rearranged ship 


The second variant ship takes advantage of the benefits 
of machinery space rearrangement. The same machinery space 
total volume is preserved but split into five spaces. The 
propulsion motors are placed very near the propellers, 
resulting in much shorter shafting runs and an increased 
GM/B ratio. The decrease in shaft length means a decrease 
in shaft weight and more space available outside the 
machinery box. The shafts previously ran through shaft al¬ 
leys that may be returned to other uses. This rearrangement 
is almost certainly not the optimum one and can be improved 
in the sense of space efficiency. It is an arrangement, 
however, that can demonstrate the benefits to be expected of 
a ship designed for electric drive. Power converters again 
connect the generators and motors. Transmission line forms 
the connections, at a much lower weight than shafting. 

Many other choices in large components are possible for 
this rearranged ship. For example, three propulsion gener¬ 
ators and gas turbines driving two or four propulsion motors 
might have been chosen. The number of prime movers was 
retained from the mechanical baseline ship, however, to make 
the comparison of transmissions realistic. Too many changes 
might have obscured the fundamental differences in 
efficiency, weight and volume. 

Both geared and direct-drive propulsion motors were 
used in this variant. When geared motors are used, the 
reduction gears are also placed near the propellers. A 
coarse layout of the rearranged ship is shown in Figure 21. 

6.7. Weight and volume algorithms 

Few components in this thesis are exact commercial 
models. The weight, volume, and other characteristics are 
taken from those for which data was available. The equa¬ 
tions for shafting and reduction gears were taken from the 
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ASSET theory manuals [28], while the switchgear and.braking 
resistor equations came from the ASSET Enhanced Machinery 
Module [27], which is not yet generally available. These 
equations represent much study by ship and equipment desig¬ 
ners, incorporating equipment which ig commercially avail¬ 
able. Where possible, the ASSET equations were verified 
against other studies and actual equipment [29, 30, 31]. 

For example, the machinery in the FFQ7 was used as a model 
and verification for reduction gears and shafting. 

Appendix E gives more explanation, as well as presenting a 
computer program used to generate weight and volume figures. 
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Chapter Seven. Machine design and system synthesis 
7.1. Machine matrix 

During this study, synchronous, permanent magnet, and 
induction machines were designed at shaft speeds of 180, 
1800, 2400, 3000, 3600, and 7200 rpm. The number of pole- 
pairs for the 180 rpm machines ran from one to twenty-five. 
For the higher rpm machines, from one to six pole-pairs were 
used. If every machine could be a generator or a motor, a 
160x165 machine matrix results. From these 27,225 possible 
combinations, two were chosen and input to the t.wo variant 
ships for synthesis in ASSET. 

These particular rpms were chosen partially because of 
the choice of the gas turbine. The LM-2500 operates at a 
full-load speed of 3600 rpm, making multiples and "nice" 
fractions of that speed desirable. A 3600 rpm, two-pole 
machine has a synchronous frequency of 60 herts, the stand¬ 
ard in the United States. A four-polo machine at 1800 rpm 
is also a 60 hertz machine. If "nice" frequencies result 
from rpm choices, results may not be obscured. The low rpm 
machine is tied to propeller rpm. For the baseline ship, 
maximum speed propeller rpm is 170 rpm. A ten rpm allowance 
for "battle override" gives a requirement for 180 rpm. 

The reduction in the number of maohine combinations is 
a bit more difficult to explain. First, it was observed 
that reduction gears add greatly to the weight and volume of 
the transmission and detract from its efficiency. Second, 
induction generators are notoriously difficult to control. 

It was then decided that generators would not be induction 
machines and any generators used would operate at the same 
shaft speed as the gas turbine, eliminating a possible 
reduction gear. The matrix then measured 12x165 and had 
1980 combinations. 

From this point on, the decision theory espoused in 
Schweppe and Merrill [13] was used, specifically using “knee 
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sets" to eliminate uncompetitive designs. The considera¬ 
tion, or Figure-of-Merit ; was always minimum weight, volume, 
and inefficiency. These external characteristics are those 
"seen" by a ship design (since propulsion current and volt¬ 
age were specified in terms of turns — fertile ground for 
another tradeoff study). Initially, individual motors of 
all of the higher rpms were considered together to select 
the best geared propulsion motors. Since the weight and 
volume of a reduction gear varies with shaft rpm, the weight 
and volume of the reduction gear was included with that of 
the machine to select the best machine-gear system. Direct- 
drive machines were also selected. The 3600 rpm synchronous 
and permanent magnet machines were also considered 
separately as propulsion generators. 

Once the initial selection of machines was made, the 
matrix measured 5x11, or 55 combinations. Ail of the com¬ 
binations were plotted in knee curves that showed the 
volume, weight., and inefficiency of each transmission. The 
volume and weight of shafting, braking resistors, cooling 
systems, power converters, and the inefficiency of the power 
converters were common to all combinations and were not in¬ 
cluded at this level. The inefficiency of any reduction 
gears was included where appropriate. Three of the gener¬ 
ators were synchronous machines and two were permanent mag¬ 
net machines. Of the motors, two were 1800 rpm synchronous 
machines, two were 1800 rpm permanent magnet machines, three 
were 1800 rpm induction machines, and two were 180 rpm 
synchronous machines. The 1800 rpm motors clearly dominated 
the higher rpm machines, .largely because of the differences 
In reduction gear weight and volume. 

Since there was no single dominating combination, ten 
of the 55 combinations were selected. These ten were among 
the best at least twice on the knee curves. These ten com¬ 
binations were composed of only synchronous machines. 

Programs to calculate off-design-point efficiency were 
written. Each motor and generator was evaluated at the 
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power level and rpm appropriate for the sustained and en¬ 
durance speed conditions of the ship, using the delivered 
horsepower (DHP) and propeller rpm taken from the ASSET syn¬ 
thesis run on the baseline mechanical ship. Some of the 
combinations had very low endurance efficiencies. There was 
a sharp division evident between the geared (lower weight, 
much higher volume) and the direct-drive systems. 

The ten combinations, with their maximum, sustained, 
and endurance speed transmission inefficiencies, were again 
made the subject of knee curves. A simple scoring scheme 
was devised to rank the combinations according to their 
grouping on these last knee curves. If a combination was in 
the beat group on a particular knee curve, it was given two 
points for that curve If it was in the second heat group, 

it received one point. If it was in neither the best or 
second best group, it received no points. When the scores 
were totaled, two combinations stood out. One was a geared 
drive system and one a direct-drive system. These two com¬ 
binations were used in both of the variant ships and are the 
subject of the next chapter. 

7.2. Knee curves 

Figures 22 through 24 show the knee curves for the 
propulsion generators The first letter of the generator 
ID indicates whether it is a synchronous machine or a per¬ 
manent magnet machine. All of the generators are 3600 rpm 
machines. The generators selected were SB (four poles), SC 
(six poles), SD (eight poles), PB (four poles), and PC (six 
poles). Generator SA was not selected because of its poor 
showing on the volume-weight curve, even though it was com¬ 
petitive on the volume-efficiency curve. 

The 180 rpm, direct-drive propulsion motor curves are 
in Figures 25 through 27. They were not combined with 
geared motors because one of the points being investigated 
was whether or not geared motors were "better" than direct- 
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drive motors. The clumping of the machines necessitated 
other graphs on different scales to distinguish between the 
machines. Machines 1-25 are synchronous, 26-50 are per¬ 
manent magnet, and 51-56 are induction machines. The number 
within a group indicates the number of poles in the machine, 
e.g., machine 32 has (32-25)x2=14 poles. Valid designs with 
over twelve poles were not achieved for the induction 
machines. Machines 3, 6, 11, and 16 were selected for fur¬ 
ther work. These are all synchronous motors. 

The higher rpm motors were oombined with their reduc¬ 
tion gears to form system knee curves. In all of these knee 
curves, machines 1-6 are 1800 rpm, 7-12 are 2400 rpm, 13-18 
are 3000 rpm, 19-24 are 3600 rpm, and 25-30 are 7200 rpm. 
Figure 28 is the volume-efficiency curve for synchronous 
machines, showing the distinct grouping of the machine-gear 
systems due to the high volume of the reduction gears. Note 
the high values of the permanent magnet machines in the 
volume-weight curve Figure 29. The clumping of induction 
motors around the low inefficiencies is shown in the weight- 
efficiency curve of Figure 30. From these curves, the 
motors on page 97 were selected. 

The initial motor and generator combinations were made 
and plotted on more knee. curves (Figures 31 to 33). On 
those curves, the high-volume or high-weight nature of the 
combinations can be seen. Since the multiple-attribute 
decision theory embodied in knee curves does not say how to 
select between high-volume or high-weight, the best of eaoh 
were selected. Combinations 1, 2, B, 9, 12, 13, 19, 20, 30, 
and 31 were chosen. The off-design-point efficiencies were 
calculated and all of the information was plotted. 

Figure 34 is a bar-graph of the maximum, sustained, and en¬ 
durance speed transmission efficiencies of the various com¬ 
binations, including the reduction gears, if any, and power 
converters. The final combination knee curves are sum¬ 
marized in the scoring scheme of Table 27, which was ex¬ 
plained on page 96 . Combinations 12 and 20 were chosen to 
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use in the ships of the study. 

This is a good method to choose among the possible 
machines. During the course of this thesis, the above path 
was followed through several complete iterations and a few 
partial ones. As stated in Schweppe and Merrill, knee 
curves serve very well to eliminate uncompetitive options, 
allowing concentration on the better ones. 

Table 27. Final combination knee curve scores 


Combo ID 

Firsts 

Seconds 

Total 

1 

3 

1 

7 

2 

3 

0 

6 

8 

1 

0 

2 

9 

0 

1 

1 

12 

6 

1 

13 

13 

4 

0 

8 

19 

1 

1 

3 

20 

3 

4 

10 

31 

0 

2 

2 


Conclusion: test combinations 12 and 20 





Figure 22. Curve of volume-efficiency for 3600 rpm gener¬ 
ators 
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Figure 23. Curve of volume-weight for 3600 rpm generators 
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Figure 24. Curve of weight-efficiency for 3600 rpm gener 
atora 
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Figure 25. Curve of volume-efficiency for 180 rpm motors 
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Figure 27. Curve of weight-efficiency for 180 rpm motors 
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Figure 29. Curve of volume-weight for geared motors 
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Figure 30. Curve of weight-efficiency for geared motors 
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Figure 31. Volume-efficiency curve for initial PM and PG 
combinations 
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Figure 32. Volume-weight curve for initial PM and PG com 
binations 
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Figure 33. Weight-efficiency curve for initial PM and PG 
combinations 
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Figure 34. Final combination transmission efficiencies 
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Chapter Eight. 


Analysis 


This chapter presents an analysis of the ships with 
electric transmissions. Standard naval architectural 
methods have been used to observe and comment on the variant 
ships, comparing them to the mechanical baseline ship. Con¬ 
clusions and recommendations follow the analysis. 

The names used to describe the various ships imply 
their internal arrangement, equipment, and ASSET Design Mode 
Indicator (DMI). The two DMIs used were ENDURANCE, when en¬ 
durance range was held constant at 5500 NM, and FUEL WT, 
when the usable fuel weight was held constant at 996.6 LT. 22 
The ship names are as follows: 


MECH 23 BASELINE: Baseline, meohanical transmission 
ELEC 23 BASELINE 12: Backfit ship, geared motors 
ELEC 23 BASELINE 20: Backfit ship, direct-drive motors 
NEW MR ELEC 12: Rearranged ship, geared motors 
NEW MR ELEC 20: Rearranged ship, direct-drive motors 
CONSTANT FUEL ELEC 23 BASE 12: Backfit ship, geared 
motors 

CONSTANT FUEL ELEC 23 BASE 20: Backfit ship, direct- 
drive motors 


8.1. Direct effects 

The direct effects of an electric transmission are the 
changes in weight and volume of the propulsion system, as 
well as the transmission efficiency. Included are the 
weight and volume of the propulsion motors and generators, 
transmission lines, cooling systems, switchgear, power con¬ 
verter, exciter, braking resistor, any reduction gears and 
their associated gear oil, and shafting. These items are 
listed in Table 31. A positive difference from the baseline 
ship means a heavier and/or larger ship. 


21. Professor John Kassakian, MIT, private communication. 

22. Not all fuel in a ship is usable. There are nooks and 
corners of fuel tanks that are inaccessible to the fuel sys- 
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Note that the only variant that haa a lower direct 
weight effect than the baseline Is NEW MR ELEC 12. The ac¬ 
cumulation of weight increases in the others makes them 
heavier, while NEW MR ELEC 12 has lower motor and shafting 
weight than the rest. Geared drive is always lighter than 
direct-drive, largely due to the high weight of the direct- 
drive propulsion motors. With respect to volume, direct- 
drive is always smaller than geared drive, beoause of the 
reduction gears. All electric transmissions are larger than 
the mechanical baseline, but the smallest variants, within 
motor type groups, are the rearranged ships. 

Table 28 contains the maximum, sustained, and endurance 
speed transmission efficiencies of the two generator-motor 
combinations. Note that the off-design-point efficiency of 
the direct-drive combination is significantly lower than the 
geared combination, even though it does not have the added 
inefficiency of reduction gears. This is in large part due 
to the poor efficiency of the slowly rotating direct-drive 
motor at the endurance speed. 

Table 28. Transmission efficiencies 

Combo Maximum Sustained Endurance 

12 0.9307 0.9266 0.8817 geared combination 

20 0.9209 0.9093 0.7754 direct-drive combination 

Endurance efficiency with one generator driving two motors. 

Table 29. Propulsion generator efficiencies 

PG ID Maximum Sustained Endurance 

SC 0.9891 0.9870 0.9737 geared combination 

SC 0.9891 0.9872 0.9768 direct-drive combination 

Table 30. Propulsion motor efficiencies 
PM ID Maximum Sustained Endurance 

S4 0.9898 0.9876 0.9526 geared combination 

SL8 0.9598 0.9496 0.8184 direct-drive combination 
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The above tables show the efficiencies of the motor and 
generator used in two particular combinations. The 
efficiency of the motor, reduction gears, and power con¬ 
verter have a direct effect on the efficiency of the gener¬ 
ator, as they change the loading point of the generator. 
Generally, motors and generators are more efficient when 
they are loaded more closely to their design point. (The 
same is true of gas turbines.) The inefficiencies of the 
power converters and reduction gears, if any, are included 
in the transmission efficiencies. 



Table 31. Direct volume and weight effects 
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ELEC 

ELEC 

NEN 

NEN ' 

ELEC 

ELEC 


23 

23 

23 

NR 

HR 

23 

23 


BASELINE 

BASELINE 

BASELINE 

ELEC 

ELEC 

BASE 

BASE 



12 

20 

12 

20 

12 

20 

Electric Propulsion Heights: in LT 






PNs 

0 

18.08 

102.94 

18.08 

102.94 

1B.0G 

102.94 

PGs 

0 

10.27 

10.27 

10.27 

10.27 

10.27 

10.27 

Trans lines 

0 

0.18 

0.18 

0.77 

0.77 

0.18 

0.18 

Cooling sys 

0 

5.98 

5.98 

5.98 

5.98 

5.98 

5.98 

Switchgear 

0 

1.56 

1.56 

1.56 

1.56 

1.56 

1.56 

Power converter 

0 

7.16 

7.16 

7.16 

7.16 

7.16 

7.16 

Eiciters 

0 

3.50 

3.50 

3.50 

3.50 

3.50 

3.50 

Braking resistor 

0 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

Reduction gears 

78.90 

41.19 

0.00 

41.19 

0.00 

41.19 

0.00 

Shafting 

69.00 

66.84 

66.84 

41.72 

41.72 

66.84 

66.84 

N298 (op fluid) 

14.7 

14.7 

0 

14.7 

0 

14.7 

0 

Direct effect 

162.60 

179.46 

208.43 

154.93 

IBS.90 

179.46 

208.43 

Diff fe Baseline 

0 

16.86 

45.83 

-7.67 

21.30 

16.86 

45.83 

Electric Propulsion voluaesi in cubic 

feet 





PNs 

0 

108.86 

909.13 

108.86 

909.13 

108.86 

909.13 

P6s 

0 

61.49 

61.49 

61.49 

61.49 

61.49 

61.49 

Trans lines 

0 

20.9 

20.9 

6.39 

6.39 

0.33 

0.33 

Cooling sys 

0 

200 

200 

200 

200 

200 

200 

Switchgear 

0 

70.2 

70.2 

70.2 

70.2 

70.2 

70.2 

Power converter 

0 

1089.73 

1089.73 

1089.73 

1089.73 

1089.73 

1089.73 

Exciters 

0 

245.02 

245.02 

245.02 

245.02 

245.02 

245.02 

Braking resistor 

0 

1422.88 

1422.88 

1422.88 

1422.88 

1422.88 

1422.83 

Reduction gears 

27J1.04 

1425.75 

0.00 

1425.75 

0.00 

1425.75 

0.00 

Shafting 

517.71 

501.5 

501.5 

306.06 

306.06 

501.5 

501.5 

Direct effect 

3248.75 

5146.33 

4520.85 

4936.38 

4310.90 

5125.76 

4500.28 

Diff fe Baseline 

0.00 

1897.5B 

1272.10 

1687.63 

1062.15 

1877.01 

1251.53 


8.2. Indirect effects 

Indirect effects are again composed of weights and 
volumes, but these arc the ripple effects of the propulsion 
system through the ship. For example, if a transmission is 
more efficient at endurance speed, it should be expected 
that less onboard fuel would be needed to achieve the same 
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endurance range as a less efficient transmission. This is 
indeed the case. Another important indirect effect is the 
change in full load displacement, which is tied to the 
powering characteristic of the ship. 23 The following tables 
lir.t the characteristics and indirect effects of the various 
ship configurations. 

Every electric transmission had lower maximum and sus¬ 
tained speeds than the baseline ship (by about 0.43 knots), 
but also lower EHP requirements at those speeds. The lower 
EHPs are a reflection of lower drafts (less resistance). 

The lower speeds show that the transmission efficiencies of 
the variant ships are lower than the mechanical baseline 
ship. There are more components in the electrical power 
trains, hence the lower efficiencies. The speed difference 
of 0.43 knots may be regarded by some as significant; it is 
about the speed reduction to be expected by a fouled bottom. 

The endurance range of all ships except those with con¬ 
stant fuel load is 5500 NM. The fuel load in the others 
varies greatly, showing the benefit of electrical cross- 
connection. In the two constant fuel ships, there was an 
increase in the endurance range of 1350 and 1400 NM, respec¬ 
tively, for the geared and direct-drive transmissions. This 
is an indication of fuel cost savings from the electric 
transmission. If a ship refuels every three steaming days 
(receiving a third of its tank capacity), steams 100 days 
each year, and fuel is priced at about $18 per barrel, this 
represents about a $600,000 savings per ship per year. 

The initial static transverse stability of the variants 
was degraded by the change in propulsion equipment. As pre¬ 
viously stated, the 6.5% GM/Beam ratio of the mechanical 
baseline ship is not as large as desired for an actual ship, 
but provided a benchmark to measure relative changes. 

tern. Typically, 95% of the onboard fuel'is usable. 

23. For the same molded lines, ships with higher displace¬ 
ments will have greater wetted surface areas and higher 
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Table 32. General ship characteristics 


CONSTANT CONSTANT 



NECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEN 

HR 

ELEC 

12 

NEN 

HR 

ELEC 

20 

FUEL 

El EC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

LBP 

422 

425 

425 

425 

425 

425 

425 

Ship disp FI 

5482 

5297.5 

5417.1 

5234.7 

5353.4 

5484.1 

5530.2 

Diff fe Dateline 

0 

-187.5 

-47.9 

-250.3 

-131.4 

-0.9 

45.2 

Nav heal 

55 

55.04 

55.02 

55.04 

55.03 

55 

55 

Nav draft 

14.44 

14.07 

14.31 

15.95 

14. IB 

16.44 

16.53 

Depth sta 10 

38 

38 

38 

38 

38 

38 

38 

Accoeeodations 

272 

272 

272 

272 

272 

272 

272 

snt/B 

0.045 

0.059 

0.058 

0.041 

0.044 

0.063 

0.041 

6H1 

950.01 

973.41 

959.1 

979.74 

947.39 

950.07 

943.41 

Full load KB 

23.23 

23.87 

23.75 

23.84 

23.53 

23.34 

23.44 

LCB/IBP 

0.514 

0.514 

0.515 

0.513 

0.514 

0.516 

0.516 

LCF/LBP 

0.545 

0.545 

0.545 

0.544 

0.545 

0.545 

0.545 

Netted surface 

23347.7 

23010 

23230.5 

22884.3 

23115.7 

23346.6 

23425.B 

Cp 

0.4 

0.594 

0.598 

0.594 

0.597 

0.6 

0.601 

Naterplane area 

17441 

17575.9 

17424.B 

17539.8 

17602 

17640.9 

17653.2 

Hull voluee 

418504 

418504 

418504 

418504 

618504 

61B504 

618504 

Okhi voluee 

121528 

121528 

121528 

12152B 

121528 

121528 

121528 

Total voluee 

740032 

740032 

740032 

740032 

740032 

740032 

740032 

Endur fuel vol 

44345 

34592 

39447 

34347 

39396 

44344 

44344 

Nachy bos vol 

109471 

109471 

109471 

109424 

109624 

109671 

109471 




Table 33. Powering 


CONSTANT CONSTANT 



HECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEN 

I1R 

ELEC 

12 

NEN 

NR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

Poxering: 

Vmi 

29.13 

29.02 

28.83 

29.09 

2B.9 

28.82 

28.7 

EHP i Vlax 

31821 

33131 

32741 

33154 

32743 

33044 

32700 

Vsus 

27.94 

27.82 

27.4 

27.88 

27.44 

27.43 

27.48 

EHP ) Vsus 

28163 

24484 

24145 

24700 

24181 

24434 

24134 

Endurance 

5500 

5500 

5500 

5500 

5500 

4550 

4084 

EHP 8 Vend 

4051 

4432 

4772 

4557 

4498 

4851 

4905 

HPi 

51550 

51550 

31550 

51550 

51550 

51550 

51550 

IcNi 

4000 

4000 

4000 

4000 

4000 

6000 

4000 

ftvg 24 hr load 

2030 

2030 

2030 

2030 

2030 

2030 

2030 
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In the constant fuel ships, the high weight of the 
propulsion motors and propulsion motors and propulsion gen¬ 
erators, combined with a smaller amount of vartically-lower 
fuel reduced GM/B by over one-half percent, a not incon¬ 
siderable amount. For the rearranged variants, GM/B 
decreased less than the two previous ships. They also have 
less low fuel, but the propulsion generators are lower than 
in the backfit ships and the propulsion motors are very much 
lower. Only a very small decrease in GM/B was seen in the 
constant fuel ships because the constant fuel load compen¬ 
sated for the increased high weight of the electric trans¬ 
mission. The longitudinal metacentric height, GM1, in¬ 
creased for all variants, though it seems it should have 
decreased with the decrease in waterplane area and draft. 

The machinery space volume was the same for all ships. 

There were no big surprises in the area of weight. The 
structural weight (W100) encloses the same volume in every 
ship, so it was about constant. The propulsion plant weight 
(W200) varied with the type of transmission. Weight groups 
W300, W400, W500, W600, and W700 were virtually identical in 
every ship, and the variable loads were dominated by the 
change in fuel weight. The Design and Builders Margin is a 
function of the light ship weight (summation of W100 through 
W700), so the margin weight moved with the light ship 
weight. The miles-per-gallon figure of NM/LT of fuel shows 
the endurance efficiency of electrical cross-connection. 

Only a few comments need be made regarding Table 35. 

The structural weight fraction shows the changes in full 
load displacement, remembering that the W100 weights were 
ail about the same. The same may be said for the weight 
fraction of the W300 through W700 groups and payload weight. 
Higher propulsion plant weights in the variant ships drove 
up the W200 fraction, except for NEW12. The fuel weight 
fraction shows the same behavior as the miles-per-gallon 
figure. 
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Table 34. Ship weights 


CONSTANT CONSTANT 



DECK 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEN 

HR 

ELEC 

12 

NEN 

HR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

Might: 
DFN onboard 

1049.1 

865.3 

937.5 

860 

931.6 

1049.1 

1049.1 

Usable fuel *t 

996.6 

822 

890.6 

817 

885 

996.6 

996.6 

Diff ft Baseline 

0 

-183.8 

-111.6 

-189.1 

-117.5 

0 

0 

NH per IT fuel 

5.52 

6.69 

6.18 

6.73 

6.21 

6.57 

6.11 

Payload Height 

571.2 

571.2 

571.2 

571.2 

571.2 

571.2 

571.2 

hi 00 

1684.8 

1686 

1692 

1660.6 

1666.6 

1686.2 

1692.1 

K200 

343.6 

351.8 

387.1 

326.1 

361.3 

351.9 

387.1 

KJOO 

236.5 

236.5 

236.5 

236.5 

236.5 

236.5 

236.5 

N400 

302.2 

302.2 

302.2 

302.2 

302.2 

302.2 

302.2 

hSOO 

615.4 

613.5 

614.3 

613.5 

614.2 

615.4 

615.4 

N600 

426.6 

426.6 

426.6 

426.6 

426.6 

426.6 

426.6 

N700 

95.9 

95.9 

95.9 

95.9 

95.9 

95.9 

95.9 

Loads 

1317.1 

1120.9 

1193.2 

1115.6 

1187.2 

1305.1 

1304.9 

DitB aargin 

463.1 

464.1 

469.3 

457.7 

462.9 

464.3 

469.5 

Disp FL 

5485.2 

5297.5 

5417.i 

5234.7 

5353.4 

5484.1 

5530.2 


Table 35, Naval architectural analysis indices 


CONSTANT CONSTANT 
FUEL FUEL 



(TECH 

ELEC 

ELEC 

NEN 

NEN 

ELEC 

ELEC 


23 

23 

23 

NR 

NR 

23 

23 


BASELINE 

BASELINE 

BASELINE 

ELEC 

ELEC 

BASE 

BASE 



12 

20 

12 

20 

12 

20 

L/D 

11.16 

11.18 

11.18 

11.18 

11.18 

11.18 

11.18 

L/B 

7.73 

7.72 

7.72 

7.72 

7.72 

7.73 

7.73 

B/T 

3.345 

3.425 

3.373 

3.451 

3.401 

3.345 

3.327 

SN1/LBP 

2.235 

2.290 

2.257 

2.305 

2.276 

2.235 

2.220 

NIOO/Df1 

0.307 

0.318 

0.312 

0.317 

0,311 

0.307 

0.306 

N200/Df 1 

0.063 

0.066 

0.071 

0.062 

0.067 

0.064 

0.070 

N300/DT1 

0.043 

0,045 

0.044 

0.045 

0.044 

0.043 

0.043 

H400/DF1 

0.055 

0.057 

0.056 

0.058 

0.056 

0.055 

0.055 

N500/Bfl 

0.112 

0.116 

0.113 

0.117 

0.115 

0.112 

0.111 

NA00/DT1 

0.078 

0.081 

0.079 

0.081 

0.080 

0.078 

0.077 

N700/DF1 

0.017 

0.018 

0.018 

0.018 

0.01B 

0.017 

0.017 

Nfuel/Dfl 

0.182 

0.155 

0.164 

0.156 

0.165 

0.182 

0.180 

Npjyloail/Dfl 

0.104 

0.108 

0.105 

0.109 

0.107 

0.104 

0.103 

Nld/Dfl 

0,240 

0.212 

0.220 

0.213 

0.222 

0.238 

0.236 

Vib/Vtot 

0.14B 

0.148 

0.148 

0.148 

0.148 

0.148 

0.14B 

HPi/Dfl 

9.398 

9.731 

9.516 

9.848 

9.629 

9.400 

9.322 

HPitVun/BFl 

273.76 

282.39 

274.35 

286.47 

278.29 

270.91 

267.53 

kNi/Dfl 

1.094 

1.133 

1.108 

1.146 

1.121 

1.094 

1.085 

Nt/DTl 

0.050 

0.051 

0.050 

0.052 

0.051 

0.050 

0.049 

Vtot/Df1 

134.914 

139.695 

136.610 

141.370 

138.236 

134.941 

133.816 

NIOO/Vtot 

5.100 

5.103 

5.122 

5.026 

5.045 

5.104 

5.122 

N200/HPi 

14.930 

15.287 

16.821 

14.170 

15.700 

15.291 

16.821 

Vib/HPi 

2.127 

2.127 

2.127 

2.127 

2.127 

2.127 

2.127 

N300/kNi 

0.039 

0.039 

0.039 

0.039 

0.039 

0.039 

0.039 

V»b/(HPi+kNi) 

1.906 

1.906 

1.906 

1.905 

1.905 

1.906 

1.906 

N500/Vtot 

1.863 

1.857 

1.859 

1.857 

1.859 

1.863 

1.363 

NAOC/Vtot 

1.291 

1.291 

1.291 

1.291 

1.291 

1.291 

1.291 

DTl/Vtot 

16.603 

16.035 

16.397 

15.845 

16.204 

16.600 

16.739 


For this technology characterization, everything 
devolves to total values. What is the total effect on the 
3hip, once the individual pieces are put together? Table 36 
gives the bottom line. The electric propulsion plants are 
larger and heavier than their mechanical drive cousin; 
however, the extra weight and volume are more than compen¬ 
sated by the savings in fuel weight and volume. If a ship 
is designed from the beginning to be an "optimized" electric 
drive ship, over 6300 cubic feet of volume and 250 LT may be 
saved. The savings might be used for other systems, to 
reduce the overall size and cost of the ship (maybe allowing 
a larger buy, since 30 ships times 250 LT is a 7500 LT 
ship), or to extend the naval architectural limits of the 
ship design. 

If a ship is baokfitted with this technology, it is un¬ 
likely that tank volume can be recovered. However, the dual 
benefits of increased time-on-station and better fuel 
economy are realized. In this case, the choice between 
geared or direct-drive systems can be made by selecting the 
system with the most leverage, i.e., if the ship is volume- 
limited, use the lower volume direct-drive system (since the 
shafts are already in place). 

To put the volume and weight savings in perspective, 
note that 6300 cubic feet and 250 LT translates to twenty 
Tomahawk missile cells. The ship would be volume limited, 
with about 200 LT of weight savings still unused. This is a 
significant addition to the firepower of any ship, and the 
unused weight allows for ship growth. 
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Table 36, Total differences 


CONSTANT CONSTANT 



RECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

MEN 

NR 

ELEC 

12 

NEN 

NR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

Total voluMit 
Fuel voluu 

44365 

36592 

39647 

36367 

39396 

44364 

44364 

fuel diff 

0 

-7773 

-4718 

-7998 

-4969 

-1 

-1 

Prpln volute 

3246.75 

5146.33 

4520.85 

4936.38 

4310.90 

5125.76 

4500.28 

Prpln diTT 

0.00 

1897.58 

1272.10 

1687.63 

1062.15 

1877.01 

1251.53 

Total 

47613.75 

41738.33 

44167.85 

41303.38 

43706.90 

49489.76 

48864.28 

Total dift 

0.40 

-5875.42 

-3445.90 

-6310.37 

-3906.85 

1876.01 

1250.53 


Total eeightii 

Dfl difference 

0.00 

-187.50 

-67.90 

-250.30 

-131.60 

-0.90 

45.20 

prpln diff 

0.00 

16.86 

45.83 

-7.67 

21.30 

16.86 

45.83 

fuel diff 

0.00 

-183.80 

-111.60 

•189.10 

-117.50 

0.00 

0.00 

other 

0.00 

,-20.56 

-2.13 

-53.53 

-35.40 

-17.76 

-0.63 
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8.3. Closure 


8.3.1. Conclusions 

This thesis has demonstrated the usefulness of electric 
drive transmissions in reducing ship weight and volume. 
Electric drive transmissions are better than mechanical 
drive transmissions on a ship basis. They provide, besides 
the weight and volume advantages, substantial arrangement 
flexibility and the opportunity to use new technologies in 
the ship design arena. Tho technical risk associated with 
these different technologies is minimal, as there is much 
industrial experience with electric machines, advanced 
switchgear, and the like. If the weight and volume reduc¬ 
tions are reinvested in the ship design through more optimum 
arrangements and subdivision, a substantially more efficient 
ship may be realized. Such a ship could successfully com¬ 
pete with the best of current ships. 

Small, light, high-power motors can be designed to a 
fair degree of detail with a computer optimization scheme if 
a meaningful objective function can be devised. For a ship 
system, the objective function should contain measures for 
volume, weight, efficiency, and relative cost (if a par¬ 
ticular material is significantly more expensive than other 
used). A steepest-descent scheme can be combined with a 
Monte Carlo scheme to quickly converge on the objective 
function. 

The use of electric drive, and its consequent electri¬ 
cal cross-connect, can reduce the endurance fuel load by as 
much as 17 5%. When usod in combination with an improved 
machinery arrangement and subdivision, that percentage car. 
rise to 181;. If the fuel load stays constant, the endurance 
range may Increase as much as 25%. 

On both an equipment weight basis arid a ship weight 
basis, systems composed of a direct-drive propulsion gener¬ 
ator (with the same shaft rpm as tho prime mover) and a 
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geared propulsion motor are better than those systems using 
no gears. Regarding volume, a non-geared system has lower 
equipment volume but higher ship volume due to the lower en¬ 
durance fuel efficiency. 

Geared propulsion motors have better off-design-point 
efficiencies than those in direct-drive systems, primarily 
due to their higher rpm. A reduction in output power (in a 
motor) of 75% means only a few percent reduction in 
efficiency for a geared motor, while the same power reduc¬ 
tion means a 20% or more efficiency reduction for a direct- 
drive motor. 

Permanent magnet machines do not appear attractive for 
ship propulsion systems. They are both heavier and larger 
than candidate systems using synchronous, and, to a smaller 
extent, induction machines. Their low air gap flux density 
is the main detractor. Current permanent magnet materials 
cannot develop the energy product to compete with other al¬ 
ternatives, even though the NdFeB magnets are now in the 
marketplace. Induction machines may be useful as propulsion 
motors, but in this thesis they did not appear so. There¬ 
fore, ship propulsion generators and motors should only be 
synchronous machines. 

8.3.2. Recommendations 

The same modeling approach used in this thesis should 
be taken with variable reluctance machines (VRM). Although 
no VRMs have been built at ship propulsion power levels, it 
is not inconceivable that they could serve in such a 
capacity. 

The induction machine model uBed here needs refinement, 
especially in the area of limiting maximum rotor current 
density. All of the machines need an analysis of their 
transient and dynamic characteristics. 

The recent advent of liquid hydrogen temperature super¬ 
conducting materials may signal an era where conventional 
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machines are overshadowed by the smaller, higher-flux 
machines possible with superconducting technology. However, 
if these new materials fail to provide the required current 
density, design of satisfactory machines may not be pos¬ 
sible . 

Integrated electric ship service propulsion plants may 
be beneficial additions to electric drive technology. Their 
influence on the systems suggested here may be an area of 
interest for future ship designs. 



Appendix A. Definitions of raaahine variables and oonstants. 


bd 

BETA 

br 

BR 

BR1 

BREM 

BRNGS 

BSAT 

CP 

CRHO 

cw 

D 

dcore 

DCU 

DMAG 

doa 

dr 

ds 

eaf 

effcy 


ew 

f req 

g 


magnet operating point flux density, Tesla. Used 
in permanent magnet machines. 

hysteresis loss factor of M19 magnetic steel. The 
figure 2.5 was used. 

air gap flux density, Tesla. Used only in per¬ 
manent magnet machine model, 
air gap flux density, 1.05 rms Tesla, 
residual induction, 7.89 kilogauss, of M19 mag¬ 
netic steel,at BSAT and T1. 
remanence flux of NdFeB, 1.21 Tesla, 
weight percentage of rotor shaft bearings, 1.03, 
or three percent of rotor weight, 
max flux density anywhere, 1.5 rms Tesla, 
stator ooil pitch. The figure 0.8 was used, 
copper electrical resistivity, 1.724E-8 ohm- 
meters . 

copper weight, kg. 

density, 7.65e3 kg/m3, of M19 magnetic steel. 

back iron depth, meters. 

copper density, 8968.0 kg/meter*. 

density of NdFeB, 7.4e3 kg/meter*. 

overall machine diameter, meters. 

depth of rotor slots, meters. 

depth of stator slots, meters. 

p.u. internal voltage, used in syn only. 

efficiency, defined as 

output power 

effcy = - 

output power + ph + pe + i2r + i2rr 

effective weight of machine, a combination of 
weight, volume, and efficiency. Used as the ob¬ 
jective criteria for the optimization scheme, 
machine synchronous frequency, Hz. 
air gap dimension, meters. 


129 





gamti 


atator geometric faotor, non-dimensional. Used to 
find convergence on active length, 
gamr rotor geometric factor, non-dimensional. Used to 

find convergence on active length. 

QMIN minimum machine air gap, 0.002 meters. 

HC1 coercive force, 0.48 oersteds, of Ml9 magnetic 

steel, at BSAT and Tl. 

11 induction machine primary current, amps. 

12 induction machine load current, amps. 

12r copper loss, watts, on stator. 

i2rr copper loss, watts, on rotor, 

jr full load rotor current density, amps/m 1 . 

Jrnl no load rotor current density, used in syn only, 

js full load stator current density, amps/m> . 

JSMAX maximum stator current density, 12e6 amps/meter 1 
ke efficiency weighting factor. 

KM magnet material cost weighting factor, defined as 

$ per pound of magnet material 

KM =-—--—------- 

$ per pound of magnetic iron 

ks[] harmonic winding factors, 

kv volume weighting factor. 

1 machine active length. 

le rotor winding space factor, used in induction 

machine. 

ler combined length of rotor end windings, 

lm radial dimension of permanent magnet, 

loa overall machine length, meters. 

lr ratio of rotor slot width to slot pitch. In per¬ 

manent magnet machines, defined as ratio of magnet 
width to rotor "slot" pitoh. 

lrat length ratio, used in permanent magnet program to 

calculate the effect of magnet overhang. 

Is ratio of stator slot width to 'slot pitch. 
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MAX_TIP_ 

minpwr 

MU 

MUR 

nr 

NU 

overhang 

P 

PC 

pe 

PF 

PI 
ph 
PS I 


rlpl 

r2pl 

rev 

relpi 

RUO 

riv 

RMA7 

rp<n 

RSF 

8CV 

aiv 

slip 


SPEED maximum allowable rotor tangential velocity, 
200 meters per second. 

minimum required power of the machine, horsepower. 
MU of air, 4we-7 Henriea/meter. 

NdFeB relative reversible permeability. The 
figure 1.05 was used. 

number of rotor bars in induction machine, 
anomalous loss factor of M19 magnetic steel. The 
figure of 2 was used, 
percentage of magnet overhang, 
number of machine pole-pairs. 

permanence coefficient of NdFeB. The figure 1.1 
was used. 

eddy current loss, watts. 

power factor of all machines. The figure 0.8 was 
used. 

3.141592654. 
hysteresis loss, watts. 

Trated/Tpullout for induction motor. The figure 

0.55 was used. 

rotor radius, meters. 

per length primary resistance, ohms/meter. 

per length secondary resistance, ohms/meter. 

rotor copper volume, meters*. 

per length Thev6nin equivalent resistance. 

electrical resistivity, 52 micro-ohm-cm, of M19 

magnetic steel. 

rotor iron volume, maters*. 

maximum rotor radius, 2.0 meters. 

machine shaft revolutions per minute, referred to 

as rotor speed. 

rotor slot space factor. The figure 0.35 was 
used. 

stator copper volume, meters* , 
stator iron volume, meters*, 
guessed slip of the induction machine. 
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slipl 

slip2 

smax 

SSF 

T1 

tmaxpl 

vlapl 

va 

vagpl 

vol 

VOLALL 

volmag 

vtpl 

w 

wr 

wm 

wt 

WTALL 

wt_iron 

wtmag 

xlpl 

x2pl 

xbeltpl 

xelpl 

xmpl 

xrdpl 

xrspl 

xa 

xsepl 

xalot 

xsdpl 


larger derived Blip, 
amaller derived alip. 
maximum machine slip. 

stator slot apace factor. The figure 0.35 was 
used. 

thickness. 0.014 inch, of M19 magnetic steel. 

per length maximum torque. 

per length ThevAnin equivalent voltage. 

IP+jQJ, VA rating, volt-amps, 
per length air gap voltage 

machine envelope volume, meters*, with a margin, 
volume allowance for frame and foundation. A ten 
percent allowance was used. 

volume of permanent magnet material, meters*, 
per length terminal voltage, 
electrical frequency, radians per second, 
width of rotor slots. 

mechanical angular velocity of rotor, radians per 
second. 

weight of copper and iron in a machine, plus a 
margin. 

weight margin for frame and foundation. A ten 
percent margin was used, 
iron weight, kg. 

weight of permanent magnet material, kg. 

per length primary impedance. 

per length secondary impedance. 

per length belt impedance. 

per length ThevAvin equivalent impedance. 

per length air gap magnetizing impedance. 

per length rotor differential leakage impedance. 

per length rotor slot leakage impedance. 

p.u. synchronous impedance. 

per length stator end turn impedance. 

slot impedance. 

per length stator differential leakage impedance. 
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xsspl 

XZBpl 


per length stator slot leakage impedance, 
per length zig-zag impedance. 



Table 37. Listing of various functions used in computer 
programs 


/St* a pseudo randon nuatier generator ttt/ 

•define MULTIPLIER 25173 

•define MODULUS 32768 

(define INCREMENT 1384* 

•define MODFLT 32748.0 

double randoif) 

< 

extern long int seed) 

seed* iHULTIPLIERIseed+INCRENENT) X MODULUS; 

return(seed/HODFLT); 


double snfln) It etator Minding factor il 

int n; II haraonic order il 

1 

double kp, kb; 

extern double cosl), sinO; 

kp=cos(0.3142ln); It pitch factor, attunes 0.8 coil pitch 1/ 

kb=(sin(0.5236tn)1/(0.52361ft); It breadth factor, fro* 

Kirtley’s ‘Basic Foreulas ... "and attuies an 
electrical Minding angle of 60' tl 

returnlkptkb); 

1 


double abs(q) 

double q; 

{ 

if (q < 0.0) 

q= (-1.0)tg; 
return(q); 


double sinh(u) 

double u; 

( 

double exp(), ans; 

ant = 0.54lexplu) - exp(-u)>; 
return Ians)| 

1 


double cotb(u) 

double u; 
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double expo, ans; 

ans = t.5t(exp(u) ♦ expl-uH; 
returnlans)j 


double ref(n) /( rotor Minding factor, eaee coeeents at sufO tl 

int n; 
f 

double kp, kb; 
extern double tinl); 
kp=l; 

kb’ (sin (0.5236ti: 1) / <0.5236tn 1; 
return(kptkb); 

> 


(define tol l,0e-6 
(define pi 3.141592654 

double besi (p,x) It bettel Korder, arg) tl 

double x; 
int p; 

{ 

double bi; 
int bial); 

double exp(),sqrt(),bis<); 
if (p<0) abort (‘ besi: negative index'); 
if <x<0) abort (* besii negative arguient’); 
if (xHO) bi=exp(xl/sgrt(2tpilx)i 
if(tx<12):ibia(p,x,lbil) bi*bi*ip,x); 
return (bit; 

} 

double bisfp,x) 
double x; 
int p; 

( 

double fabtd; 
double xx; 
int i,fk,k; 
double bi‘<); 
double t*l; 
xx=x/2.; 

for (i*l;i<prl|ir+) t*tlxx/i; 
if lt<l.0e-M) bi*0; 

el it 

{ 

bi*t; 

xpxxliti 

♦or (kM)(k<1901)ttHfikilt)-)akt(bi>toll>>0)|krr) 
( 

ffc-kt(pek); 
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t-tlxx/fk; 

bi*biH| 

) 

> 

return (bil; 

) 

int bialp,x,pbi) 
double x,tpbi; 
int p; 

I 

int iit,fn,k,}kj 

double **,b!,t; 

double fabsU,sqrtO,*xp(i) 

fnMtpkpj 

t=1.0; 

bi-1.0| 

xx=0.125/x| 

for (k-1 j(k<30)W((fabiltl-fab» , bittolli>0)jk++) 

I 

fk-(2lk-l)l(2»k-l)j 

t*tlxxllfk-fn)/k| 

bi=bi+tj 

i 

if (k==31) isr=lj 
else 

{ 

ii^Oj 

bi=bilexp(x)/sqrt(2.0lpilx)( 

> 

Ipbi=bi; 
return (i&t.-p 
) 

double besk(p,x) It aodified bessel function Kplx) I/ 

double xj 
int p| 

I 

double bk) 

double expo,iqrt(),kOO,kill) 

if (p<0) abort l a bcski negative index *)j 

if (*<0) abort I* bnki negative arjueent *>j 

if lx>40) bk-«xpl-x)/ii,t t(2.Otx/pi)| 

else if (p»M)) bk*kO(i)j 

else if (p«l) bk'kllxlj 

else 

double gO ( gl,gj| 
int j; 
gO»kOlx)| 
gl=klfx)i 

for (j*2)j<pel[j**) 

{ 

gj»2l(J-lHgl/**g#| 

gO«glj 
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gi*gJi 

> 

bk*gjj 

) 

return Ibklj 
) 

double kO(x! 
double x| 

( 

double bk| 

double log(),sqrtU,e»p()j 
if lx<l) 

{ 

double t,b,i,c,p0,«2j,f,hj,rjj 

l«t Ji 

b*0.3tx; 

a*.5772354iflog(b»; 

c*b<ks 

«0=-e; 

*2j a l| 

Ml 

hj=0; 

for (j'lij<7ij**l 
1 

rj L 1.0/j| 

x2j*«2jlc) 

f=flrj»rj; 

hj=hj*rj| 

gO*gO+«2jlftlhj-a)| 

) 

bk*gO; 

) 

else 

1 

double tU2],a,b,c,pa; 

int 1| 

e»exp(-x!) 

b'1.0/xj 

c'sqft(bl) 

tl01*b| 

for <1=1; 1 < 12j 1 ♦■+> tUMU-lltb) 

pa « 1.2533141-.15M442lt[0J| 

pa ♦* . MI811128<t1H-.09139095tt C2)| 

pa +" . 1344596IH31-.2299830lt[41j 

pa .379241W3J-.J24727/ftl4]| 

pa ** .5575348itm-.4262433»tC8)s 

pa ♦« .218431Bltt9]-.06480977ltn0J*,0O91893fl34tIlll| 

bkMlctpai 

) 

return tbk); 
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double kl(,il 
double *j 

( 

double “ikj 

double log(),*«p(),s(jrtllj 
if lk<i/ 

( 

double a,b,c,|l,x2j,f,hj,rj| 

list j; 

b-\/2} 

a*.377215M+loglb>; 

c*blb; 

*2j*b| 

f=l| 

hj'lj 

gl=1.0.'*+«2jl(,5*»-hj)| 
for (j=2jj<»|j*0 
f 

*2.'»«2j;ci 
rj*i.O/j| 
f=flrjlrjj 
tij=hj+rj; 

gl*glex2j4ft(.5el*-MI4j!; 

1 

bk=gl) 

1 

eise 

{ 

double a,b,c,tU2],p.ij 

int l| 

e=«»pl-«){ 

b=l.0/») 

c**qrt(b)( 

tl51*b| 

for (1*1|1<12[1*+1 Ull*tlM]|b» 

pa » 1.233IKU.44W?7l>tOJ-.144WflJlttl3| 
pa ♦« .126042 7 tt(2]-.173i432lt(3]+.2847bl8tt(4]| 
pa -.45945424USl+.‘i283381lt[t]-.4W22WIU7]| 
pa +‘ .50502:i9nt81-.2581J04«ttV3| 
pa ♦* .07880001 ttC 101 -.010324t8W!U| 

bk«alc4pa; 

> 

return (bk); 

) 

double besip <p,argl It deriv of bee Jtorder, arg) »/ 

int pt 
double erg; 

{ 

double x,y t zj 
i 3 ben (p-l,arg)| 
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y * p I b(5i (p,arg)/»rg| 

1 ■* * • y| 
return (z); 

5 

double beskplp,arg) /( deriv of bes kiorder, erg) 1/ 

int p| 
double irg; 

( 

double x,y,z| 

x * -best (p-!,erg)| 

y * - p I beslc (p,»rg)/arg| 

2 » x ♦ y; 
return (i)j 

> 
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Appendix B. Synchronous Machines and General Relations 


An equiva3.ent circuit for a synchronous machine is 


Ra Xs 

+ o-VWNAA-UUU&i 

Ia 



Figure 35. Synchronous machine equivalent circuit 

where Ra is the stator resistance and Xs is the synchronous 
reactance. The internal voltage of the machine, Eaf, is 
developed between the stator and rotor across the air gap. 
The current direction is shown as if the machine is a motor. 


Vta = Rala + JXsIa + Eaf 


Eaf represents a mutual coupling between the stator and 
rotor, and 


Eaf 


« M If 
~ 42 ~~ 


where 


M = 


4_Po ki kf 1 r Ms Nf 

* g P* 


If Is field current, Nf is the number of series field turns 
and kf is the fundamental rotor winding factor. Since it is 
never desired to specify the number of turns on either the 
stator or the rotor, a scheme has been devised so that 
derivations are conducted in volts-per-turn, ampere-turns, 
and ohms-per-tum-squared, which results in power in watts. 
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Figure 36. Phase belt conductor area 

For the stator, the transverse area can be divided into 
phase belts. The area labelled 'A* is the area occupied by 
phase in one direction. Then, Hs la = Aa Ja , where Ja 

is the current density in phase ’a’ conductors. The area 
has some effective conductor area, subject to the need for 
conductor cooling passages and insulation area. Therefore, 
for stator and rotor currents, analysis yields 


It is also desired never to specify the number of slots on 
either the stator and rotor. Accordingly, slot space fac¬ 
tors are defined as 


slo-t pitch _ 


Figure 37. Slot space factors 

stator slot width 

Is =- 

stator slot pitch 
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rotor slot width 
rotor slot pitch 










For a typical turbogenerator conductor bar, the copper area 
is about thirty-five percent of the conductor envelope area. 
Variables titled SSF (stator slot space factor) and RSF 
(rotor slot space factor) embody this thirty-five percent. 
The conductor area of a single stator phase in one direction 


Aa = 


(r+g) 2rr Is ds SSF 


Af = 


r 2 tt lr dr RSF 


The armature resistance is 


and, for the rotor, 


There is only one phase on the 
rotor. 


Ra 


rho It Ns 
WNs~~ 


where rho is the electrical resistivity of copper and It is 
the turn length. If a stator turn can be modeled as 



Figure 38. Stator turn 


and the circumfrence is 2rr(r+g), then 

u - 2i ♦ 2{ *--L**l*lLV*L } 

p 

where 4(4/3) is (l/sin(60 # )) and CP is the coil pitch. Then 

_Ra__ - -L^.L 1 + (2ir(r+g)CP/p)4(4/3) } 

Ns 1 " w"(r+g)~Is"ds'siF".. 
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The stator and rotor copper losses are 


Ral»* = rho It As J»» Rflf» = rho ltf Af Jf» 


The rotor is full pitched 

The vector diagram for the synchronous machine equiv¬ 
alent, circuit is (assuming Ra ==> 0) 



Figure 39. Synchronous machine vector diagram 

By use of the law of cosines, 

Eaf* = Vt* + (Xsla )* - 2VtXsl« cos(r + 90» ) 

where r is the power factor angle. If this equation is put 
into per unit form, with eaf = Eaf/Vt and xs = Xsl»/Vt, then 

eaf* = 1 + xs* + 2 xs sin(T) 

xs can be calculated, allowing the calculation of eaf. 

There is simple linear relationship between eaf, the no-load 
rotor current density (jrnl), and the full-load rotor cur¬ 
rent density (jr). It is Jr = (eaf)(jrnl), because Eaf is 
directly proportional to If and If is directly proportional 
to Jf. Arpere's Law states 




n dA 


If one chooses an integration path around half of the rotor, 


Po H 2g = Po 


4 Nf If kt 
tt p 
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where P is reluctanae, For this conductor configuration, 

„ „ Ho 1 No 1 

Pselt =-ds wn and Pmutual = -ds ws, 

3 2 

with ws equal to slot width. Since ws is not known, we use 
the stator slot spaae faotor, Is, multiplied by the number 
of slots to yield. 

(slots) ws = Is n (r+g) Therefore, 


Lai 


5 Ns* Po 1 ds Is tt (r+g) 
---- 


The real power developed by the machine ia 


Pwr = 3 Vt Ii pf, where pf is the power factor. Through the 
use of Lenz* Law and Ampere's Law, terminal voltage may be 
expressed as 


Vt 


2 r l_w Br k« 
p 


Using our previous relation for la, the expression for Pwr 
is 


Pwr 


2ir r 1 w Br k* Ja SSF (r+g) Is da pf 
P 


Finally, winding factors need to be derived. The wind¬ 
ing breadth factor, kb, is 


sin(m n r/2) 

kb n = - 

m sin(n r/2) 

where r is the electrical angle between adjacent slots, n is 
the harmonic order, and m is the number 'of slots per pole 
per phase. The winding pitch factor, kp, is 
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kpn = sin (n a/2) where a is the eleotrical angle 


between sides of the ooils (pitch angle). For a three phase 
winding, ot = 2r p. If a 0.8 coil pitch is assumed (which 
will rid the machine of certain harmonics during balanced 
operation) then « = 180°(1 - 0.8) = 36® and 
kpn = sin(0.3142 in). Assumptions are needed to calculate 
kbn without specifying the number of stator slots or turns. 
If the winding angle is specified as did Kirtley [32], then 

sin(n 8w/2) 
kb " ~ (n 0w/2) 

A reasonable electrical winding angle is 60° , since most of 
the stator periphery will contain turns. The breadth factor 
devolves to kbn = sin(0.5236 n)/(0.5236 n), for which the 
fundamental harmonic factor equals 0.965. The winding fac¬ 
tor is the product of the breadth and pitch factors. For 
the rotor, a pitch factor of one is assumed. 
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Table 38. Listing of synchronous machine design program 


•include ‘■stdio.h 1 
•include “def.h’ 

It pregraa naaei syn.c tor synchronous, round rotor aachines »/ 

lent int <Kdj It start pi:int for rondo* nuaber generator 1/ 

double b[261I26), htll][2il, ksCB.'l, krlSl; 

It b is “bast* array, h is “hold" array, ks/kr are winding factors tl 

aainl) 

i 

double design_point(), rnd.walkO, swfl), rwfl), ka, kv, ainpwr, 
stapsiza, randoat), abed, frag, rpe; 
int p, iteration, i, j, bast, print.outl), loops; 

FILE If open (), Ifp; 

printf(“\nReading input data froa SYM.DAT . . 


fp=fopen("syn.datVr“); It input seed tar rando* nuibars 1/ 
fscanf(fp,*Zd*,lseed)j 
fscanflfp, “Id*, bp); 
fscanflfp, "Ilf, iainpwr!; 
ainpwrl*74i,0; 
fscanflfp, “Ilf“, bkel; 
fscanflfp, “Ilf“, bkvl; 
fscanflfp, “Ilf, brpa); 
fclasalfpi; 


It input nuabar of pole pairs 1/ 

It input aaihina power, derived fa ASSET i/ 
It convert to watts 1/ 

It CERs for Effective Haight II 

It aachinc aax shaft rp* 1/ 


printf(*\nHow aany loops do you want? *!; 
scanfCZd*, bloops); 

printf(*\n\nDoing prograa calculations . . .in'); 

for (i-i; j < B; i*=2) It haraonic winding factors tl 

C 

ks(il«swfli); 
krtihrwf (1); 

5 


freq c rp*tp/60.0; It aax electrical frequency 1/ 

It HAIN BODY OF THE PR08RAH 1/ 

for (i*l; i <* loops; +«-i) 

( 

stapsize’0.1; 

iteration^; , 

design.point(ainpwr, p, ka, kv, freq); 

It put stuff in the hold array tl 

while (iteration <= 10) 
t 
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rnd.ualkiainpnr, p, stepsize, Ice, kv, freq); 

It stagger around 1/ 

test=0{ It index to best EM of the let 1/ 

for (j=i; j<=10; *+j) 

if (htjHlBl < htbestKlSJ) 

best=j; It find the beet eachine 1/ 

if (absf (htOJriSJ - hibestlllBll/hlOltlBl) < O.OOS) 

ft stall uprovetent in EM 1/ 
f 

steps!ze/=2.0i 
♦♦iteration; 

> 

else It transfers best to 0 position tl 

( 

for (j*l; j <= 25; ++j) 

hfOKjl = NCbestlCj 3; 

> 

1 

for (j«lj j <= 25; ♦♦)) 

biiHjl-hibestlCjlj It keep the best eachine tl 


best-1; 

for li=l; i <* loops; r+il 

if (biiinai < bchestineii 

best=i ; It find and keep the best of ihe best tl 

einpNr/=74di.O; It turn back into hp tl 

printout (best, p, linger, ke, kv, rpe); It output to disk file tl 

fp'fopenCsyn.dat*,VI; /* output seed tl 

fprintf(fp,"Id*, seed); 
fprintf(fp,"\nld', p); 
fprintf (fp/inllf, einpwr); 
fprintfIfp/VnXlf, ke); 
fprintf(fp, a \nUf*, kv); 
fprintf(fp,“\nXlf*, rpe); 
fdose(fp); 

1 

It END OK MAIN PR06RAN; ALL THAT F0LL0M ARE FUNCTIONS tl 

double design.pointieinper, p, ke, kv, freq) 

It detereines a rasdoe design point tl 

double ainpur, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js, Is, lr, dcore, ds, dr, q, e, 1, is, eaf, i2rr, 
va, ph, pe, i2r, vol, «t, effey, in, xsl, xs5, xs7, xsal, 
cm, siv, rev, riv, rev, loa, doa, et.iron, find.iivl); 
extern double sqrtd, randoell; 
int c*0, d>b; 


( 

r*randonO»RHA)(; 

if ((2tRIIrlfrtq/pl < KAl.TIP.SPEEO) It check rotor tip spied «/ 
break; 

> 


w=2IPllfreg; 


lr=randoe()t0.5 a 0,25; 

It rotor slot factor tl 

dr-randoeOlr/5.0; 

It slot no deeper than 201 of rotor radius tl 

dcore=!BPIr)/(8SATIp); 

It back iron depth tl 

ds : randoi()l0.9tdcore; 

It slot depth < 901 of body depth tl 

while Ic ! = 1) 

It gap dieension tl 


( 


g*randoe()»(0.1lr - 6 HIN) a SHIM; 
if (9 > 0 ) 

break; 

> 

ls-randa*()t0.5 ♦ 0.25; It stator slot factor I/ 

js-randoiOUSMAI; It full load stator current density 1/ 

jrnU(BRIglp)/( 8 IHU»RSFIridrllr»kr[UI; It noload rotor current density 1/ 

xsl*kstU»kstllj 
«s5=lks[51tks[5]/25); 
xs7*(k5t71lks[7]/49)| 
xsal=(541slPlids/18); 

xs= <HUIjstSSFtlr+q>t(rkg)tdstPItls t U2l(xsi 
/U2irt8RlksU)); 

It 

if (xs > 2.0) 

continue; It 

else 

♦+d; 

eaf'cqrt.ll a xstxs ♦ 2lxs40,ol; It 0.6 is sin(t), pur factor angle, 
eaf is p.u. internal voltage at full load tl 

jr*eafljrnl; It jr full load, linear with eaf V 

l*(ninpwripl/<2IBRJwiksm»rlSSF»d*«PI»(ragms»jstPFI) It active length tl 

va s 2iPItrillwlBRIksUHj*ISSFt(rag>idstls/p| It va rating tl 

siv • find sivd, r, g, ds, dcore, Is); It stator iron voluoe 1/ 

riv < UPIlrKr - 2tdrtlr); It rotor iron voluo* t/ 

scv > 2IPlt(rkgltdstlst(l ♦ 2.309UPI*tr+g)ICP/p); 

It stator copper voluoe tl 
rev « 2IPlirldrllrlU a 2,30944Pllr/p)| It rotor copper voluee tl 

cw • (rev a scvMDCU; It total copper weight tl 

loa * 1 ♦ 41 frag) ; it length-over-ell tl 

doa • 2*(r+gado+dcorei 5 It over-all-dieeeter tl 

vol > VOLALLI(loalPlidoeldoa/4); It eachine envelope values tl 

wt • MTALLKcw a DHBRNSStriv a siv)); It eachine weight in kg tl 


a xsS a xs7)/(PIIglp) a k sal>) 
p.u. synch iapedance tl 
don’t want xs too big tl 
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et.iron * DKriv ♦ siv); 


It iron Might only tl 




ph * 31.B62:5»8EMfreqiB3IIHCi»«Uron/D; 

It hysteresis loss in Mtts, lies iron Might of aanhine tl 

pe - (106234.9INU!SShTtBSATIT 1 IT I If reqlf raqlut ,iron)/(RHOID) ; 

It eddy curmnt loss in Mtts, uses iron Might oi sachin* V 

i2r = 2.0ISSfldslCRHOIPII (1 ♦ 2.d<M«PH(r*gl»CP/pl IjsIjiNr+gllls; 

It stator copper loss in Mtts 1/ 
It revised t-12-87 il 

i<rr = 2.0IRSFtdrlrllrtjrljrlCRH0tf 1 II (I * 2.3094«Pl*r/p*; 

It rotor excitation losses, 1-12-87 1/ 
i2r+=i2rr; It total copper losses 9/ 


pifcy=(»inper)/l»inper ♦ ph ♦ pe a i2r); 

e*=*t * kel(l-etfcy) * kvtvolj It Effective Might 1/ 


htonil’js; 

h[0H2)=freg; 

htOJ[31=H) 

htOlW-r; 

h(0K5]rg| 

hlOKil^dcore; 

h[0K7hds; 

hI01[81»dr; 

HtOHfMs; 

hlOHlOl-lr; 

hlOJUIl’vol; 

htCMJ21»nt; 

htOH13]=ph; 

htOHMl-pe; 

h[0][15>i2r; 

hlOK161*va; 

ht0U17}»e»fcyi 

htOKlBlMu; 

htOJtl9M; 

ht0K201*jr; 

h(01C2lJ-jrnl; 

hr0H221<xs; 

ht0U231*e»fi 

htOH24]=loa; 

h[0][25)=doa| 


It this section just changed all the variables in the 'hold' array tl 


return; 

)> 


double rndjialMsinper, p, stepsize, ke, kv, freq) 

It Mlks .'bout design.point 19 tins 1/ 
double stepsize, einper, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js, Is, Ir, dcore, ds, dr, g, », 1, ss, »r‘, (2rr, 
va, ph, pe, i2r, vol, tit, efdev, e«, isl, xs3, ss7, xsal, 
cr, siv, sev, riv, rev, loa, doa, nt.iron, find.sivi); 
extern double sqrtil, randoai); 
int iM; 

uhile ii <■ 10) 

( 

It read in the xalk around the design point tl 

js'htOHlllll ♦ stepsizellrandaeO - 0.S11; 
i( (js > JSHAD 

js > JSISAK; It r eset to lioit tl 

»«2IPItfr#qi 

r<it01[41l(l ♦ sttpsizel'.randoal) - 0.5)1; 
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if lUtr/pl ) ItAMIP.SPEE!)) 

continue; It po to next try if violated tl 

q»htQH314U t stepsizet(randoet) - 0.5I)| 
it (g < GIUN) 

g«BHIN| It relit tj the liait 1/ 

dcore^(Bftlr>/(HSATtp); It Mit efficient ui« of iron tl 

ds=hE0K7]|(1 i stepsizellrandoeO * 0.S)>) 

if Ids > dconl It can't have too-deep slots 1/ 

di*dcore; It rcsst to the liait t,' 

dr=hC0K83t(1 ♦ stepsizellrandoeO - 0.3)1; 

1 s a hCOHVlt11 i stepsizellrandoeO - 0.5i); 
if (Is > 0,75) 

ls*0.75; It resat to thr liait tl 

if Its < 0.25) 

ls*0,25; 

Ir-hCOlUOltd ♦ stapsizetlrendoil) - 0.5)); 
if Hr > 0.75) 

lr=0.75| It resat to the liait 1/ 

if Hr < 0.25) 

If0.25t 

It computation section of the talk tl 

jrnl=(BMglp)/(8IKUIASrlrldrtif llrrtl]); It noload rotor current density 11 

isl'ks'lllkstl); 

xs5«(ks[31ikst5)/25l{ 

xs7*(ks[7MsM/4?l; 

xsal=(5llslPlids/16)! 

xs>(HUIjsiSSFt(rrg)i(rtg)idslPltls t (12t(«sl * ssS ♦ xs7)/(Pllgip) * xsal>) 
/(12«rlBK«ksm>; 

It p.u. synch ispedance tl 

if fxs > 2.0) 

continue; I) can’t hare is too big 1/ 

eaf=sqrt(l ♦ xslxs ♦ 2txsl0.6); It 0.0 is sinfT), per factor angle, 

eaf is p.u. internal voltage at full load I/ 

jr*eaftjrnl; It jr full load, linear sith eaf tl 

Msinpnrlp'/UtWIMkimirlSSFIdsJPMIr+glllstjsIPFl; It active length 1/ 

va s 2tPllrtlMBR>k«il]<jsiiSFI(r+g)ldsilsrp; It va rating V 

siv * find s.vll. r, g, ds, dcora, Is); It stater iron volute tl 
riv * llPUrllr - 2ldrtlr); It rotor iron volute tl 

sev = 2IPll(rVg)ldsilsl(l + 2.30fMP!t(rsg)ICP/p); 

il stator copper volute tl 
rev - 2IPIIrtdrtlrl(l ♦ 2.30?4fPIIr/pl; It rotor copper volute tl 

ct = (rev + scvllDCi); It total corpse wight »/ 

loa * 1 ♦ 4Hr+g); It length-over-ill tl 

doa = 2ilrrgvds+dcore); /( over-all-diat*ter tl 

vol * VOlALUilnatPUdoaldoa/4); It tachint envelope volutt tl 

tt = UTALLKct v DKBSifPGIriv t siv)); It aachins wight in kg tl 
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nt.iron > DS(riv + siv); 


It iron diight only 1/ 


ph = 31.B6225IBETflt4req*BRl*HCJ*wt_i ron/Di 

/l hystaretis lots in Mitts, uses iron Might of lachine tl 

p« > ()GB73&.9INUIB5ATIBSATITliTltfreqtfreqtut.iron)/(RH0ID)| 

/I oddy current loti in uatti, him iron Might of uchini tl 

i2r • 2.0ISSFidslCBH0>- Ii II > 2.3094h .'»(r+p)«rP/pl IJiljillrigm*} 

/t ititor coppor loti in iatts 1/ 
It revised 1-12-87 »/ 

i2rr = 2.01flSF»dr>rklrljrljrlCRH0tPI* II * :.3094IPltr/p)| 

It rotor etc i tit ion lotus, 1-12-97 tl 
i2r+-i2rrj /1 total -opptr losses 1/ 


effcyliinpnrl/liinpur + ph + pe + i 2r >; 

eu*ut r Mll-effcyl ♦ kvkvol, 1 It Effective Might 1/ 


himil-js; 

hti H23sfreq; 

h(iK31»Nt 

hf; 3t<l*r( 

hliKSlf) 

httK61=dtore; 

h. r iK71°dt; 

hlillBl'dr) 

hCiH91*ls| 

hlilU01=lr| 

hamn*voi, 

hlilll?]=Nt; 

hmCUl'ph; 

htiltlfl'pe; 

nti)[15)»i2r} 

htiKUl'va; 

ii[iK171*effcyj 

htimaineu) 

hti][19K; 

hli](20]ejr| 

hti]I21]=jriilj 

hliK221=xsj 

hlil f 231*eafj 

hfi][241*l3a| 

htil[251»doa| 


It this secuon just changed all the variables in the 'hold' array tl 


♦+ij It go to the next htilll tl 

1 

return; 

1 


print.out(test, p, einpet, ke, kv, rpi) 
int best, p; 

double ainpw, ke, kv, rpe; 

( 

char outfilelMl; 

FILE Ifpo, IfopenO; 
int i| 

prinif findhat is the naM of the file ehere you want the output? ')) 
scanfl'ls*, outfile); 

fpo*fopen(outfile, Vlj 


fprintf<fpa,'Zd', p)| 
fprintf({po,*\nIlf‘, linpur); 
fprintf(fpo,'\nIlf", ki)| 
fprintf ifpo,'\elir, kv)( 
fprintflfpo,'\nIlf', rpi); 
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(or li=lj i <» 25| *>i) 

(print!l(po,‘\nXlf , ,btbestlfi))| 
(phntdtpo, *\n’I; 
fdoieUpuli 
1 


double find.iivll, r. g, as, dcore, Is) 
double 1, r, g, ds, dcore, !« 

{ 

double one, tuo, three, (our; 

on: * (rrgtde+dcorelMrrgrditdcore) - (r»g)llr*g)i 
teo * 2IP!l(r+g)(d*liSj 

three “ (r+g+d**dcore)l(r+g*ds*dcore) - (rrgedi)l(rrgrde>| 

four - 1 1 (PIPone - tec) *■ PII4»(r*gitthre*j 
returntfcurlj 
1 


153 






Table 39. Listing of synchronous efficiency program 


(include 'stdio.h" 

(include 'def.h' 

It proqrau naiei seff.c to find efficiency of synchronous tachines tl 
It ttorks with only i single eachine I/ 

Mini) 

t 

FILE tfopen(>, tfp; 

double r, jrnl, jr, js, Is, lr, dcure, ds, dr, g, u, 1, xs, eaf, i2rr, 

va, ph, pe, i2r, vol, tit, effcy, eu, xsl, xs5, til, xsal, kstOl, 
siv, riv, et.iron, find_siv(), peeff, psjs, pirpi, dhp, rp», tinpttr, 
ke, kv, freq; 

extern double sufU, sqrtl); 
int e=0, f, p, i; 
char infiletld; 

for <i=l; i < 8; i♦=2) It hanonic winding factors tl 

l 

kstilnuf !i)j 

> 

xsl-ksmiksLll; 

xs5*lksC5Mst5J/25l; 

x»l*(kst7m»l71/4Vl! 

isal=(5«ls(Pllds/l8); 

printfI'NnCalculates efficiency of a single eotor.tn'); 
while le !* 1) 

{ 

f *0; 

printft’Khat is the nam of the input file? *1; 
scanfCXs', infile); 
fp = fopenlinfile, V)| 

fscanflfp, 'Id', Ip); It input nuiber of pole pairs (/ 

fscanflfp, 'Xlf, tainpwl; 

einpnr N 7(6.0; It nc« in uatts it 

fscanflfp, *Xlf“, Ikal; 

fscanflfp, 'Ilf', fckvi; 

fscanflfp, 'Xlf', tirpa); 

fscanflfp, 'Xlf, (js); 

fscanflfp, 'Xlf', (freq); 

fscanflfp, 'Xlf', In); 

fscanflfp, 'Xlf, lr); 

fscanffp, *Xlt‘, Igi; 

fscanflfp, ‘Xlf, Idcore); 

fscanflfp, 'Xlf', Ids); 

fscanflfp, 'Xlf*, Id - ); 

fscanflfp, 'Xlf, lls); 

fscanflfp, *X!f, Ur); 

fscanflfp, 'Xlf*, Ivol); 

fscanflfp, 'Ilf, let); 
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fscanflfp, ‘Ilf*, t,ph)| 
fscanflfp, ’Ilf, Ipd; 
fscanflfp, 'Xlf', li2r)| 
fscanflfp, "Xlf*, tiva); 
fscanflfp, "Xlf', lieffcy)) 
fscanflfp, "XH", ten)) 
fscanflfp, "XH*, till 
fscanflfp, "XH*, Xjrl) 
fscanflfp, "XH 1 , Ijrnlli 
fdoselfpij 

s*v = find.sivll, r, g, ds, dcore, ls)j It stator iron voluee 1/ 

riv * llPItrllr - 2tdrtlr); It rotor iron votu»» 1/ 

etjron * DKriv r siv)t It iron Height only 1/ 

Nhile It !« 1) 

I 

printf(*\nNhat is the sustained speed eachine horsepouer? *>j 
scant CXlf', tdhp); 
dhp 1° 744.0) 

print!I'Hhat is the sustained speed eachine rpt? *>) 
scant(‘Jit*, iperpal; 

treq - perpelp/40,0; It aax electrical frequency tl 

pejs = jsldhp/ainpurj It PH stator current tl 

xs=(HUIjslSSFt(r*g)t(r+g)lditPllls t <12l(xsl ♦ xs5 ♦ xs7)/(Rttgtp) ♦ xsal)) 
/(12trtBRIks(l]); It p.u. synch iepedance 1/ 

eaf=sqrtll ♦ xslxs + 2»xsl0.4); It 0.4 is sin(t), pur factor angle, 
eat is p.u. internal voltage at full load tl 
jr=eaf t jrnl ; It jr full load, linear Kith eat 1/ 

ph = 31.B422518ETAIfreqlBfU»HCl»Ht_iron/Dj 

It hysteresis loss in eatts, uses iron Might of aachine tl 
pe = (104234.9»HUlBSAT»BSAWltTllfreq»frsqtetjron)/<RH04D)| 

It eddy current loss in uatts, uses iron wight of eachine 1/ 
i2r = 2.0*SSF*ds*CRK0»PI* 11 + 2.3094tPHIr+gl4CP/p) Ijstjsllregllls) 

It stator copper loss in uatts 1/ 
i2rr = 2.0*RSFIdr*r»lr*jr»jr*CRHO*Rl* (l ♦ 2.3094IPllr/p)| 

It rotor excitation losses, 1-12-87 tl 
peeff = dhpMdhp ♦ ph ♦ pe + i2r t i2rr); 
printfl'tn Sustained speed efficiency is Xlf, pieffl) 

printfCVnhhat is the endurance speed eachine horsepouer? ')) 
scant <*Xlf\ lidhp); 
dhp t* 744.0) 

printf fHhat is the endurance speed eachine rpe? *)) 
scant CXlf', Iperpel) 

pejs * jsldhp/einpur; It PH stator current tl 

freq * perpetp/40.0) It sax electrical frequency tl 

xsMHUtjstSSftlrsgltiragltdstPItls t I124(xsl ♦ xs5 a xs7)/iPllgtpl ♦ xnalH 
/<12trtBRtksC111; It p.u. synch iepedance 1/ 
eafasqrtll a xslxs * 21xs»0.4l) It 0.4 is sinlrl, per factor angle, 
eat is p.u. internal voltage at full load 1/ 
jr*eaftjrnl; It Jr full load, linear eith eaf 4/ 

ph * 31,B42251BETAlfreqlBRlIHC1let iron/D) 
pe * (104234.9INl)4BSATtBSATtT14Tltfreqtfreqlet.iron)/(RHOfDl| 
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12 r • 2,OI5SFIdslCRHOIPII II ♦ 2. 3094 IP I * (r+Q) ICP/p) Ijiljillngllli 
12rr * 2.0tRSFtdrtrllrljrtjrtCRHOtPH II ♦ 2.3094IPllr/p>| 
peeff * dhp/(dhp ♦ ph ♦ pi * i2r ♦ i2rr> j 
pnntfHn Endurir.ee efficiency it tlf‘, p»*ff)| 

printf CVnSiee eichine? ‘Ij 
scinfCId', U); 
if (I «» 0) 

continue} 
eli* if If » 21 
{ 

* * lj 
brut) 

> 

) It Md of f~loop 1/ 

) It end of *-loop tl 

) It end of Mio progrie 1/ 

double find.sivll, r, g, di, dcort, It) 
double 1, r, g, di, dcort, U| 

( 

double one, tuo, three, four) 

one 1 lrtgtds+dcore)l(r+g+di*dcore) - (ng)llrig)} 
tuo * 2IPH<r*glldslUj 

three = Ir+gidi+dcorelllr+g+di+dcore) - (rig+d*)l(rig+d»)| 

four = UlPIlon* - tuo) ♦ PIMKrigltthreej 
return(four); 

> 


lb6 


Appendix C. Permanent magnet machine 

The equivalent circuit for a permanent, magnet machine 
is almost identical to that of a synchronous machine. The 
only difference is the source of the internal voltage, which 
develops the field flux wave that interacts with the arma¬ 
ture flux wave. The field flux wave is a result of per¬ 
manent magnets built into the rotor to develop magnetic 
polos. 

A typical magnetic circuit, combined with Ampere’s Law, 

shows 


Helm + Hfg = Ni 


Figure 40. Typical magnetic circuit 



We use a constitutive law of B = PH and assume that any 
steel has P = •. If flux leakage is about aero, BmAm = BgAg 
since flux is solenoidal. Then, 


Ni - Hgg . . Po Hg A* 

lm =- and Am =--- 

Hm Bm 


The magnet volume is 


Vm - lm Am 


(Ni - Hgg) Po H* Ag 


Minimum magnet volume occurs when the magnet’s maximum 
energy product (MEP), HmBra, is a maximum. If current is 
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zero and only magnitudes are used, 


y 



_g_Bg» 

Po H», Bm 


The load line of the magnet is developed as 


Im g Ba 
A» ~ Mo ill Ag 

Bm lm Ag 

H* g Aa 


Pc 


and 


The permnnance coefficient, Pc, is the slope of the load 
line. On a magnet diagram 



Figure 41. Magnet operating point diagram 

A good algorithm for machine design is to adjust the 
geometric dimensions to operate at the MEP, on the load 
line. If operation at MEP is assumed, the needed slope is 
determined, the dimensions are randomly generated, P o H* is 
calculated, and the design is maximized for MoHg and mini¬ 
mized for magnet volume, then a search technique has been 
delineated. 

Magnets may "overhang" the active length at either end 
to account for manufacturing imprecision and to permit a 
smaller armature diameter. This overhang affects the 
developed flux. 
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Figure 42. Permanent magnet maohlne diagram 

The flux per pole, 0, with overhang and a given Ja and Jm, 
la the same as would exist in a configuration in which 
Ja' = Jm’ = Ja + ND and there were no overhang. N is a non¬ 
linear parameter promulgated in graph form by Ireland [33]. 

A good power fit for N ia 


N = 0.38558 (- )0. tobu 

D 


= overhang and 0 i overhang 1 0.34. Then, 


Ja S Jm 

5 ' D~ 

Ja 

-+ N 

D 


0with = 0w/o Y 


This flux-with-overhang ia applied to the problem aa would 
be the usual flux. What is the usual flux? A permanent 
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magnet hysteresis diagram shows the residual flux density to 



Figure 43. Permanent magnet hysteresis diagram 


The magnet does not operate at that point, but rather 
on the load line. Also, there is not magnet material at 
every point along the circumference of the rotor. A Fourier 
series is a good way to find the flux-without-overhang. 

Using a developed rotor, 


M 


Ztr r 
P 



Figure 44. Magnet material on a developed rotor 


PM = Bd 

wr p wr 

lr = _ = - 

nr/p rr r 

• n it x 

B(x) = T, An cos(---), 

n = 1 T 


T = 


2irr ' 
P 
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2 

An -- 

T 


nrrx . 

cos (-) dx for n = 1, 2, 3 .. 

T 

B(x) = 

Bd 

for 

0 i x i (nr lr/p) 


0 

for 

(nr lr/p) i x i (nr/p) 


Bd 

for 

(nr/p) i x i (nr/p)(l+lr) 


0 

for 

(nr/p)(1+lr) S x i (2nr/p) 


After integration and manipulation, 


An 


2Bd it lr ir lr 

- [ sin(-----) - cos(---—) + 1] 

tt 2 2 


The equation represents only the fundamental term of the 
Fourier series. The flux-without-overhang is 


2Bd i lr tf lr xp 

Br = - [ sin(-) - cos(-) + 1] coa(-) 

TT 2 2 2r 

As is usual in steady-state analysis, the magnitude is used. 

The next quantity to find is Bd. If the magnet 
operates on the load line, the operating point flux is 

Brem 

Bd =- where Mr is the relative 

1 + Mr/Pc 

reversible permeability of the magnet and Brem and Pc are as 
previously defined. With this relative permeability, the 
magnet length is lm = g Mr . The magnetic machine can now be 
specified. An end view with dimensions is given in Chapter 
Four. 
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Table 40. Listing of permanent magnet machine design program 


tincluds 'stdso.h' 
♦include 'def.h' 


It prograt new: paa.c lor pertinent tagnet aachines 1/ 

long int seed; /J start point for rai.dot nutber generator 1/ 

double bt 101}C273, hCUK27], ksIB); 

It b is 'best* array, h is ‘hold’ array, ks are Minding factors tl 

tain!) 

f 

double design.pointO, rnd_walkl>, seft), ke, kv, ke r linger, 
stepsize, randotO, abs<), freq, rpa; 
int p, iteration, i, j, best, print.outO, loops; 

FILE IfopenO, tfp; 

printf('\nReading input data froa PHH.DAT . . .•>j 


fp=fopen(*pia.dat','r'); It 
fscanf (fp,'Id',listed); 
fscant(fp, 'Id', Ip); 
fscanf(fp, 'Ilf, ttinper); 
•inpur 1*746.0; 
fscanf (fp, "Ilf, Ike); 
fscanf(fp, 'Ilf', Ikv); 
fscanf tf p, 'Ilf, likt); 
fscanf(fp, 'Ilf*, Irpa); 
fclose(fp); 


input seed for randot nutbers tl 

it input nutber of polt pairs 1/ 

It input tachine poser, derived fa ASSET 1/ 
It convert to uatts tl 
It CEfis for Effective Height 1/ 


It taebine tax shaft rpa tl 


printfC\nKo» aany loops do you uant? *); 
scanfv'Id*, Hoops); 

printff\n\nDoing prograa calculations . . .\n’>; 

for (i=i; i < B; it*2) It haraonic winding factors tl 

{ 

kstibsef (i); 

1 

freq=rptlp/40.0; It aax electrical frequency tl 

It MIN BODY OF THE PR06RAH II 

for (i«l; i <= loops; t+i) 

{ 

printf(*\»i»Xd*, i); 
stepsize : 0.1; 
iteration*0; 

design_poi.it (tinpur, p, ke, kv, ka, freq)p 

It put stuff in the hold array 1/ 

while (iteration <« 101 


rnd.ealkltinper, p, stepsize, ke, kv, ka, freq); 
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/I stagger around tl 

best=0; It index to best Ed of the lot tl 

tor (j=lj j<=10; ♦+)) 

if (htjlUBl L htbeotittBl) 

best=jj It find the best uchine tl 

if (absUhiOKiai - htbestm8])/(htO)tl8Jv.001)1 < 0.005) 

It seal 1 iiproveient in Ed tl 
It 0.001 takes care of div by zero t/ 

[ 

stepsize/’2.0; 

♦♦iteration) 

) 

else It transfers best to O position tl 

{ 

for lj=l; j <= 26; ♦♦)) 

htOltjl - hlbestltjl; 

1 


1 

for fj=l; j <= 24; ♦♦j) 

bti][j)=htbest)tj]j ,'t keep the best uchine tl 

) 

best-1; 

for fi=l; i <= loops; ♦♦!) 

if IbtiHlS] < bCbestltiai) 

hest>i; It find and keep the best of the best tl 

»inper/=746.0; It turn back into hp tl 


printout (beet, p, einpur, ke, kv, ki, rp«); It output to disk file tl 


fp=fopent , pM.dat , ,'» , l; It output seed tl 

fprintf(fp,‘Id', seed); 

fprintf(fp, , \nld' l pi; 

fprintf(fp,*\nllf', ainpur); 

fprintf<fp,'\nXlf*, ke); 

fprintf<fp,'\nXlfkv); 

fprintftfp/Vnllf", k»); 

fprintf(fp,*\nXif", rpel; 

fdosetfpl; 

) 

It END OF MAIN fR06R.')H; ALL THAT FOLLOW ARE FUNCTIONS 1/ 


double design.pointleinper, p, ke, kv, kt, freq) 

It detereines a randoe design point tl 

double einpur, ke, kv, ki, freq; 
int p; 

( 

double r, js, Is, lr, Icore, ds, g, u, 1, as, va, ph, uteag, 

pt, i2r vol, wt, effcy, •«, xsl, xs5, «s7, xsal, 1e, voleag, 
bd, find.lrl), cu, siv, scv, riv, loa, doa, find.iivO, et.ironi 
extern double sqrtd, randoel); 
int c=0, d»0; 

ehile !d !* 1) ( 
unde (c !* II 
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{ 

rTandoiOIRRAX; 

if li2>P!trlfreq/p) < HAX_T1P_SPEED) It check rotor tip speed tl 
break; 

> 

<i«2IPHfreq‘, 

bd“BREH/(l + MUfi/PCI| It air qap Hus, operating point, fro* aagnet 

characteristics 1/ 

lr*find_lr(bdl; It find uqnet slot factor t/ 

atule (c !- 1) /I gap diaension tl 

{ 

g=randoa() >10. Hr - GHINI a GRIN; 
if (g > 0) 

break) 

) 

la=g>Rl)R; 

ls=randoa()>0.5 ♦ 0.2S; 
j s-rando*()IJSHAX; 
dcore=(bdt(rala))/fBSATtp); 
ds=randoi()>0.9>dcore; 

xsHksmtksll); 

«s5Hkst5)>ks.l51/25)5 
xs7=(ksl7i>kst71/49); 

»sal=(5>lslP!>ds/18)| 

X5 ! (RU>js>SSI : >(rtgH*)>;rtgtla)>0-lPI>ls 

I il2Hxsl ♦ xs5 ♦ xs7)/(PItqtp) t xsal>>/(12t(r>l■>tbdtksUI>; 

It p.u. synch iapedance tl 

if (xs ) 3.0) 

continue; It don't itant xs too big tl 

else 

++dj It *y escape hatch -- aission coaplete tl 

lMainp*r>p>/(2>bd>«>kstn>(sHa>>SSF>ds>PI*(r*gHOm>js>PF>; 

It active length tl 

va=2tPIt(r+;.a)tltw«bdtks( 1 ItjstSSFtfr>ga-la)tdstls/p; It VA rating 1/ 


It radial length of aagnet, see PR-8, 12/15/86 1/ 

It stator slot factor tl 

It full load stator current density 1/ 

It back iron depth 1/ 

It slot depth < 90X of body depth tl 


siv - find_sivll, r, g, ds, dcore, Is, la); /> stator iron voluae tl 
riv » llPHrlr; It rotor iron voluae tl 
scv x 2>Pt>(r+g+la>>ds>li>(l ♦ 2.3094*P11(r+g+la)>CP/p)j 


c» « scvtDCU) 

loa ■ 1 ♦ 4>lr+gHa); 

doa * 2XrrgH**dsHicore); 

vol • VOLAlLHloalPIldoaldoa/f); 

at » KTALLIlca * SIlBRNSSIriv ♦ siv)); 
»t_iron * SMIriv ♦ siv); 


It stator copper voluae tl 
It total copper Height tl 
It length-over-all tl 
It ovar-all-diaMter tl 
It aachine envelope voluae tl 

It aachine aetght in kg >/ 
It iron aeight only tl 


ph < Jl,B4225l8ETA>fr»q>8Pl>HCllat.iron/0; 
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It hysteresis loss in mitts, uses iron Height at eichine tl 

pe * (106236. OINUtBSATIBSATtTUTMfreqlfraqtmt.ironl/lRHOtD); 

It eddy current loss in mitts, uses iron Might of uchint 1/ 

volng>2tP!lrllllellr; It eignet voluie tl 

nt»«g=volng»DNABi /• eignet Height tl 

12i = 2. OtSSFtdstCRHOtPlI (1 + 2. 3094*PI *(r+g+1■> *CP/p > tjs»jsl(r*geli)lls; 

It ntitor copper lass in mitts tl 
It remised 1-12-87 tl 


effcyMiinpmrl/dinpmr ♦ ph ♦ pe * i2r)| 


em=mt ♦ kitmtug ♦ kejll-effcyi ♦ kvKvol + kelvoleig); 

It Effective Height 1/ 


hlCKil-jt; 

ht0K61=dcore; 

htO]lll]=vol; 

hlflH161=vij 

h[OK20}=H<teig| 

hlOU231=0.0; 

hlOU26]-doi; 

It this 


hlOK2)=freqj 
ht0][71=ds; 
h[0K121>Ht; 
hCO][17J=effeyj 
h[0H211-volng; 


h[0113]=m| h[OHO s rj 

hl0]t81=lej h[0K91»Ui 

hlOHl3)=ph; htOlMhpej 

h[0][ 181=em; hlOlUSl'l j 

ht01C22]««s; 


It this one likes overhing=0 tl ht4H243-bd; 


hlO)t5J*gi 
ht03tl01-lrj 
h[OK151=i2r j 


h[OK25Hoi( 


section just chmged ill the viriibles in the ‘hold* irny tl 


return; 

)) 


double rndjulMtinpmr, p, stepsize, ke, kv, ke, freql 

It milks ibaut design.paint 10 tins 1/ 
double stepsize, einpmr, ke, kv, ke. freq; 
int p; 

{ 

double r, jc, Is, lr, dcore, ds, g, m, 1, xs, ve, ph, pe, i2r, 

val, mt, effey, em, xsl, »s5, ms7, xsil, le, valeig, mtug, bd, 
findjrll, find.brl), hr, overhing, Irit, cm, siv, sev, riv, loi, 
doi, findsivl), mt.ironj 
extern double sqrtl), rindoel); 
int i»i, j y eheck=0, ccjj 


while li <= 10) 

{ 

/I reid in t) mill: tround the design point tl 

ji>h[0Kl]l(! * stepsizellrindoiU - 0.5))) 
if (js > JSHA1) 

js » JSNA1; It rent to lioit tl 

m<2tPItfreg| 

r'htOJHJMl ♦ stepsiiwNrindoeO - 0.5)); 
if ((mlr/p) > MAX TIP SPEED) 
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conUn.it; It go to not try it violated tl 

9*h(0]£5)l(l * stepsiieKranooaO - 0,5)1; 
it <g < GttIN) 

g-GHIN; /I reset to the halt 1/ 

overheng*randoel>l0.34; It aagic nutter '0.3V trot the book 

on pertinent eigntts bv Janes Ireland, 
the effect, of overhang is to concentrate air 
gap flu*, or reduce leakage. XI 

l*htO)Mlt 

la'glNUR; It nee la, based on the neu g V 

bd=PRFN/(l * WIR/PC); It aagnet operating point flux density M 
lr=findjr(bd); It find nett Ir t/ 

ls*hi01[9HI! -r stepsizel IrandoaO - 0.5)1; 
if (Is > 0.75) 

ls-0 75; It reset to tke liait 1/ 

if (Is < 0.2S) 

ls>0.25; 


for !j=l; j <* 10; ♦♦)) It august char convergence loop >I 

( 

lrat>l/(2fir*g*la)l; 
br=bd*findj)r(nverhang, irat), 

It find effect of overhang on air gap flux tl 
dcnre=(brl(rtlel)/(8SVfp); It east efficient use of iron tl 
ds^hiOirntdcore/hiOlU); It keep thee in toe sate ratio 

is upon exit ‘a designjtt tl 


xsl*kst.UCkstlJ; 

xs5*(ks[51tksC51/25l; 

xs7=(k3[7]»kst7)/4V; 

xsal=(5lls»PIIds/l8); 

xs=(HUIjs»SSF»(r+gHi)»(r+g+la)ldslPltls 4 <I2t(xsl ♦ xs5 ♦ xs7)/(PI»glp> 

♦ xsall)/(12i(rela)lbrtksil]); 

It p.u. synch ispedarce tl 

l=(ainparlp) / '21brlslksnil(r<!a)ISSFtdstPII(reg+la)llsljstPF); 

It active length tl 

) It end of aagret char convergence loop 1/ 

if Ixs ) 3.01 
( 

♦♦check; 

if (check > 251 It can't close on decent is 1/ 

i k 

for (ccj*i; ccj <» 10; ee-ccj) 
hlccj1C 181*10000000.0; 

It sake Uis rondos walk undesirable tl 

printfl* burp'); 

break; It go to next design point tl 

) 

cuntinue; It can't have is too big tl 

> 
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va«2IPU(rM»)lllvlbrikstl JlJstSSf't (rtgvletldslls/p; 


It VA rating 1/ 


mv = find.sivU, r, g, ds, dcor*, It, Ul| It stator iron voluse ti 

riv * HPttrlr; It rotor iron voIum tl 

lev -- 2tPH(rtg*l»lldil!slll > 2.30941PIHrvgvl»)tCP/pt) 

It stator copper voIum tl 
c* 3 KvtSCUf It total copper wight 1/ 

lot * 1 ♦ 4»lr*g+l»l| It Ungth-nver-al) tl 

dot * 2ll.'igMitd*>dcor«); It over-all-diasettr tl 

vol * V0lAU.l(loalPIIdoaldoa'4)j It itciiint envelope voIum 1/ 

nt * KTALl.Mce ♦ DtlMKSSIriv ♦ siv)) | It aaebint Might in kg 1/ 

lit.iron = D*Sriv ♦ siv)| It iron wight only tl 

ph * 31.B62?5IBElA4fr»qlBRllHCl»*t.iron,'D| 

It hysteresis lost in wits, uses iron wight, of sachint tl 

pe = (106234.9tNU»8SATtBSATiT14Tlt<r»q»fr«ql»t_iron)/<RH0ID>t 

It eddy current low in wtts, uses iron wight of nachine tl 

vnl«ag°2tPllrlilll.0Mov;rhangi(2l(r«gvlil)))il*ilr| It eegnet voluw tt 

ut»ag=vaUaglDflA 6 ; It aagnet wight tl 

i2r = 2.0kSSFIds«CRH0IPIi (1 ♦ 2.3094«PI»(r+gMe>tCP/p> 4js»Js»(r+gfle)llsj 

/( stator copper loss in eatts 1 / 

It revised 1-12-87 1/ 

effcy s l»inpnr)/t»tnp*r + ph + ps * i 2 r); 

en-Mt + kiinteag + kel(S-efvcy) * kvltvol ♦ kelvoleag); 

It Effective wight t' 

nil K 1 l=j s; hliU 2 ]=freq) MillSl-m; hmk43-rr MiHSl«g| 

hliH61=dcore| htiJI71*iis| hliltSWa; MiHfMtj MiJllOlMr; 

htilUU'vol} MiJU23*»tj htiK131>phi MiH!41=pe; hti H151^i2r; 

htiltUl'vii hliJiniwffcyi hCiU181=e»! hltHlBHf 

MiH20]**tiag; hliH211=vol«agi htiH221>»s| 

nli11231-overhang; htiK247=br; htilt251*loa; hliJ1261*doi! 

It this section just changed all the variables in the "hold* array tl 

++i; It go to the next htilM tl 

check= 0 ; 

> 

return; 

) 


printout (best, p, linpur, ke, kv, k*, rpi) 
int best, p; 

double linpur, ke, kv, ks, rpa; 

t 

char outfilelMl; 
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FILE Mpo, IfopenO; 
int i; 

prints(“NnMhat is the naee of the file xhere you uant the output? ’ 
Mantels*, outfi'e); 

fpo*fnpen(outfile, *w"); 

fprintf (fpo/ld*, p)| 
fprintf(fpo,"\nXlf", ainpur); 
fprintf(fpo, a \nXlf‘, kel; 
fprintf<fpo,'\ftXlf*, kv); 
fprintf (fpo/lnllf*, ki); 
fprintf(fpo,*\nIlf*, rpa); 
for (i=l| i <= 26| e+i) 

fprintflfpo, , \nllf‘ ) btbestHil); 
fprintflfpo, *\n*); 
fclose(fpo)j 
) 


double find_lr(br) It bracketing tl 

double br| 

( 

extern double cosO, siml; 
int c=0; 

double 1r=0.5, top*l , 0 , tbr=0.C, ids'), bta*0.0; 

uhile lc!=l) 
t 

tbr=(2*Bfl£H* (sin(0.5?Plllr) - cos(0.6tPHlr) ♦ l.Oll/PIi 
if (absl(br-tbr)/br) <■= 0.001) It check for convergence I 
break; 
if (hr > tbr) 

( 

bt««lr; 

lr=bt« + (top- btnl/2.0; 

> 

else /I br < tbr »/ 

( 

top»lr; 

lr*bte ♦ (top-bt«)/2.0; 

} 

if (abs((top-lr)/lop) <> 0.0005) It no convergence tl 
( 

printf('\nXlf*, lrl; 
abort(‘no solution for lr'l; 

) 

) 

return(Ir1) 

1 
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double find.br(overhang, lrat) 
double overhang, lrat) 

{ 

extern double po»i)| 

double expon=0.704133501, factor'd.J85576838, enn, yyyj 

/I enn is Ireland's N, yyy is Ireland’s V 1/ 

enn = factor t pon(overhang, expon)| 
yyy * (lrat ♦ overhang)/(lrat ♦ enn); 
return(yyy); 

) 

double find.sivd, r, g, ds, dcore, Is, la) 

double 1, r, g, ds, dcore, ds, Is, la; 

{ 

double one, two, three, four; 

one = (r+g+la+ds+dcoreUlr+gdavds+dcore) - (r+gvle)t(r*gda)| 
tuo * 2«PI*(rtgvle)ldsll5j 

three = (rtg+le+ds+dcorelKr+g+la+dsrdcore) - (rag+la+ds)f(r+g+la+ds); 

four = 11(PItone - teo) * Pll4t(r*gHa)lthree; 
return(four); 

) 



Tabic 41. Listing of permanent magnet efficiency program 


•iDelude "stdio.h" 
♦include "def.h* 


It prograt nate: pefl.c to •find efficiency oi pertinent tagnet tachines 1/ 
It norks with a single eachine 1/ 


lain!) 

{ 

FILE Ifopenl), Ifp; 

double r, js, Is, lr, dcore, ds, dr, g, n, 1, xt, It, kt, sgtffM.O, 

va, ph, pe, i2r, vol, tt, effcy, on, xsl, x$S, xs7, xsal, kstBl, 
siy, riv, nt.iron, find.sivl), pteff, pgeff, teff, pgjs, ptjs, 
pirpt, geff=0.98, pceFf=0.99, gr, dhp, rpt, tinpnr, ke, kv, freq, 
nttag, voltag, overhang, bdi 
extern double shH), sqrtl); 
int e=0, F, p, i; 
char infilelUl; 

siv - find_siv(1, r, g, ds, dcore, Is, It); It stator iron volute tl 

riv = llPItrtr; It rotor iron volute tl 

Nt.iron a Dllriv ♦ siv); It iron Height only tl 


printFI'\nCalculates efficiency For a single tachine.Vn’); 
while (e !* 1) 

{ 

F ■ 0| 

printf(‘Nhat is the nate of the input File? *>; 
scanfl'Xs", infile); 
fp = fopenlinfile, V); 

fscanflfp, 'id*, lip); It input nutber oF pole pairs 1/ 

fscanFIFp, "Xlf, tiinpwr); 

tinpur t= 746.0; It now in tatts tl 

fscanflfp, “Xlf*, Ike); 

fscanflfp, nif", Ikvl; 

fscanflfp, "Xlf*, Ikt); 

fscanflfp, *XIF*, Irpt); 

fscanflfp, a XH a , Ijs); 

fscanflfp, a Xlf a , Ifreq); 

fscanflfp, "Zlf", It); 

fscanflfp, a Xlf a , lr); 

fscanflfp, “Xlf*, tig); 

fscanflfp, ’Xlf", tdcore); 

fscanflfp, "Xlf", Ids); 

fscanFIFp, a Xlf a , lit); 

fscanflfp, a Xlf a , tlsl) 

fscanflfp, -Xlf, Hr); 

fscanflfp, ■Xlf", kvol); 

fscanflfp, ‘Xlf, Int); 

fscanflfp, 'Xlf, lph)| 

fscanflfp, "Xlf 1 , kpe); 

fscanflfp, "XU*, ti2r); 

fscanflfp, a Xlf a , kval| 


170 




fscanflfp, 'XU*, ieffcy); 
fscanf(fp, •Ilf, 4»»l; 
fscanflfp, ’Xlf\ tl); 
fscanflfp, 'Xlf*, inteag); 
fscanflfp, ’Xlf’, Avoleigl; 
fscanflfp, ’Xlf’, ixs); 
fscanflfp, ’Xlf’, ioverhang); 
fscanflfp, ’Xlf’, tbd); 
fclose(fp); 

nhile If !- I) 

{ 

printf C\ndbat is tht sustained speed eachine horsepower? *>| 
scanfl’llf’, (rdhp); 
dhp 1° 746.0; 

printfCHhat ii the sustained speed eachine rp»? •)) 
scanfl’Xlf*, iperpe); 

freq - perpelp/AO.O; It eat electrical frequency tl 

pajs = jsldhp/einpttr; It PH stator current tl 

ph « 31.8A22StBETRIfreqtBRStHClt»t_iror./D; 

It hysteresis loss in Matts, uses iron Might of eachine tl 
pe * (10A23A.9IKl)i8SATIBSATITllTUfreqlfreql*t_iron)/(RHDtD)( 

It eddy current lass in Matts, uses iron Might of eachine tl 
i2r = 2.0tS5FldslCRH0IPII (1 ♦ 2.3094»PI*(r +g+la) *CP/p> »jstjsl<reg*leMs; 

It stator copper loss in eatts tl 

peeff = dhp/(dhp + ph + pe + i2r); 

printf(*\n Sustained efficiency is Xlf*, peeff); 

printf(’\n\nHhat is the endurance speed eachine horsepoeer? '); 
scanfl’Xlf’, Idhpl; 
dhp 1= 746.0; 

printf(‘What is the endurance speed eachine rpe? '); 
scanfCXlf’, iperpe); 

pejs = jsldhp/einpwr; It PH stator current tl 

freq = perpetp/AO.O; It eat electrical frequency tl 

ph * 31.8A225IBETAtfreql8RUHCIIet_iron/D; 

pe = (10A23A.9thUI8SATtBSATITUTllfreqlfreql»t.ironl/(RH0ll))| 

i2r * 2.0ISSFtdslCRH0IPII (1 ♦ 2,3094IP11 (r+Q+l"e) *CP/p) ljs»j»»(r+gtle)»ls; 

peeff = dhp/(dhp * ph ♦ pe ♦ i2r); 

printf("\n Endurance efficiency is Xlf", peeff); 

printf C\nSaee eachine? *1; 
scanfl'Xd’, if); 
if (f == 0) 

continue; 
else if If *= 2) 

C 

e = 1; 
break; 

> 

It end of f-loop tl 
It end of e-loop tl 
It end of eain prograe 1/ 


1 

1 

) 
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double lind.sivd, r, d, ds, dears, Is, !•> 
double 1, r, 9, ds, dcore, It, li; 
l 

double one, t«o, three, lour) 

one = (rt 9 rUtdiedcQr«)llr+QM«*di*dcore) - (r*9+li)l(r+9+l»)| 
tea = 2 IPll(r* 9 +li)«dillt| 

three ’ (r»44i»*dt»dcore)tlf*9»le+di*dcor#) - <r+9+l«td«)»(r+9+l»*-e»), 

lour ' It (PI tone - teec) ♦ Plt4t<r+g+l*) Ithreej 
return(lour)) 

) 
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Appendix D. Induction maohines 


The equivalent circuit for an induction machine is 
shown in the figure. 


R. X, 



Figure 45. Induction machine equivalent circuit 

The stator leakage reactance, XI, is the sum of the 
stator differential leakage reactance (Xsd), the stator slot 
leakage reactance (Xss), and the stator end turn leakage 
reactance (Xse). In turn, Xsd is the sum of the belt and 
zig-zag reactances. 

The rotor leakage reactance, X2, is the sum of the 
rotor differential leakage reactance (Xrd), the rotor slot 
leakage reactance (Xrs), and the skew leakage (Xskew). Xm 
is the magnetizing reactance and He is the core resistance. 
R1 and R2 are the stator and rotor resistances, respec¬ 
tively. 

A derivation of the properties of electric machines, 
using Ampere's Law and the constitutive relation earlier 
postulated, yields the result that 


flux = E I 

n * II kf 1 

k=0 


Mo 6 H« J 1 r 

-15.1 

g tr n* v* 


coji(wt) 
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From this important stator inductances are taken. 

Belt leakage reactance (Xbelt) is the sum of the reac¬ 
tances due to phase-belt harmonics of an "infinite" slot 
winding. In most machines, the most important harmonics 
present are the fifth and seventh, as the third is canceled 
in balanced operation. Then the belt reactance is 


6 Po w N« 2 1 r ka 5 ks7 

Xbelt = --- [(-)* + (-)* ] 

n g p 1 5 7 


This is the harmonic form of the fundamental mutual reac¬ 
tance of Appendix B. The winding factor is again the 
product of the pitch and breadth winding factors. 

Zig-zag reactance is leakage due to all the air gap 
harmonics that would be produced if the winding had one slot 
per pole per phase. For a phase belt of one slot, with each 
slot carrying the same current and equally separated in time 
and space phase, the zig-zag reactance alone would be 
present. Belt leakage occurs because phase belts are ac¬ 
tually several slots wide. Zig-zag reactance has harmonic 
orders higher than seven, with the same form as Xbelt. No 
even or triplen harmonics will be present. 

The fundamental harmonic of the flux yields the mag¬ 
netizing reactance, Xs, which can be viewed as that required 
to "energize" the air gap. 

6 Po wN» 1 rka* 


Figure 46 shows a typical stator slot. The stator slot 
leakage reactance, summing the self and mutual reactances, 
is 
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where ns is the number of stator slots and Np is the coil 
throw, in slots. Np = (CP ns)/(6 p), a result available by 
manipulation. 



ws 

Figure 46. Stator slot diagram 
Then, 


Xss = 


3 H Po 1 Nb* 
ns 


d2 


w2 


ds 
2 ws 


Ml + CP) 


and 


(1 


Is) 2 tt (r+g) 
ws 


where we have traded the need for the knowledge of the num¬ 
ber of stator slots for the need to know the width of an in¬ 
dividual slot. A reasonable relationship between slot 
dimensions is d2 = ds and w2 = ws. Then, 


3 Po w 1 Nb* ds (1 + CP) 

Xss -- 

Is 2ir (r+g) 

The stator end turn leakage reactance may be estimated 
by treating the two end regions as a single helically shaped 
winding. If the active region of the machine is ignored and 
the helix given air core properties, the inductance can be 
found from standard sources. 


-16 w ir Po r* sin* ($w/2) Nr* pirr prrr 

Xse = - - -!p * —)Ki, * (—---) 

®w* le p* le le 


where le is the combined length of both end windings. 
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Presume the winding radius to be r. with helix pitch 
= n/3 and 

le = 2irr tan(ir/3) Np/ns 

Ip’ and Kp’ are the first derivatives with respect to their 
arguments of the hyperbolic Bessel functions Ip and Kp. 

When the three phases are summed, a multiplier of 1.5 will 
be realized. Finally, using the previous results for Np, 

81 w Po Na* r prrr „ . , p*rr v 

Xse =- Ip ’ (-) Kp (- ) 

pw 1 CP tan(ir/3) le le 

Stator resistance is R1 = (rho It Na)/(Aa/Na), as in 
Appendix B, leading to 

6 CRHO Na* (1 + -1(4/3) 2ir(r+g) CP/p) 

R1 = - 

SSF ds tt (r+g) Is 

When the actual calculation is performed, 1 is not known. A 
guos 3 -and-iterate scheme is used. Iteration continues until 
convergence on 1 is achieved. 

Rotor resistance uses a similar scheme, but the 
presence of rotor bars and end rings instead of turns 
changes it somewhat. A model of a rotor bar is below. 



wr 


Figure 47. Typical rotor bar configuration 

Induction motor transformer models provide a way to 
find rotor resistances and inductances. The flux density 
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produced by the stator and rotor is 


B =-(Fa + Ft ) 

g 

3 Ne ka Wa 

The stator mmf is Fa = -j -Ti 

tt p 


where Ti is stator 


current and wa is slip frequency. Rotor mmf is 


Ft 


ws r 1 Br 
Za~ 


where Zn is the rotor surface impedance. The air gap volt¬ 
age, or voltage across Xm, is defined in terms of the flux 
density and rotor mmf as 


Vag = 


w 2 p 1 Nb ka 

-j - Fr 

wa r 


If rotor mmf is now identified with rotor current referred 
to the stator winding, 


Fr 


-j 


3 Nb ka 
W p 


12 , where la is rotor current. 


Then, 


Zb w 6 1 Nb* ka* 

Vag = ---la 

wa tr r 


Separating Zb into it3 real and reactive parts and using a 
rotor surface modal to describe the relation between rotor 
electric field amplitude and rotor surface current yields 


12 1 Na» ka* 

R2 = - ra lot 

nr 

12 1 Nb* ka 1 

Xrs = - xaiot 

nr 

6 Po w 1 Ns* r ka 1 1 1 

Xrd =-E {(—--—)> + 

ir g p+k*nr p-k*nr 
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If magnetic diffusion is ignored, end ring resistance 
can be calculated by comparing losses in the rings and slot 
The ratio of current densities is found by the ratio of the 
areas. This is then squared and multiplied by the ratio of 
volumes. When summed, 


R2 


12 1 N.J 
nr 


niot 


Cl + 


nr r wr 
tt I ler p* 


niot = CRHO/(dr wr), nr is the number of rotor bars 
and ler is the end ring length, approximated as 
ler = 2rr (r - wr/4 - ds/2). The rotor bar width, wr, is 
found by specifying nr and observing that (nr wr) = (2 it lr). 

Rotor skew leakage arises when the rotor slots are 
skewed angularly along the axial length to prevent rotor 
cogging. Then, flux does not fully link the bars. When the 
effect is integrated over the rotor, it is seen that 

Xskew = Xm [1 - (2 sin(skew/2)/skew) i ] with the amount 
of skew measured in radians. When typical values of skew 
are input, we see that (Xskew/Xm) » 0.5%. This is a negli¬ 
gible effect and will be ignored. 

From the previous rotor bar model, it is seen that 
xalot = w Mo (d2/w2 + dr/2t»r). Assume that d2 = wr/4 
and w2 = wr/4. Then xaiot = w Mo (1 + dr/2wr). 

Fitzgerald et al [23] state that only small errors 
result if Rc is omitted. Therefore, the core branch may be 
omitted. 


Once the components of the equivalent circuit have been 
calculated, the designer must turn to power and torque con¬ 
siderations. The internal mechanical power of the machine 

is 


F 


(i_:_siip2 

slip " 


3 l2> R2 


The air gap voltage has been previously defined. ‘The tor¬ 
mina) voltage, Vt, may be found from Vag by means of a volt- 
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age divider and a Thevenin equivalent circuit developed. 



Figure 48. Induction machine Thevenin equivalent circuit 


R1 + j(XI + Xm) 

Rel + jXel = (R1 + jXl) in parallel with jXm. Then operat¬ 
ing point torque is 


_ 1 3 Via* (R2/slip) 

w (Rel + R2/aIIpP ~+"(Xei~+ X2)*~~ 

Torque is a maximum when the power delivered to (R2/slip) is 
a maximum. By matching load ana Thevenin impedances, the 
power is a maximum and a slip-at-maximum-torque is found 

R2 

smai - -(b o1 , + (x e r+ - X2)*~jo. * 

and the corresponding torque is 

1 _ 1 - 6 Y lal 

w Rel + (Rei> _ + (Xel + X2j* jo. s 

Typical induction motors have the ratio between maximum and 
operating point torque as 


T 

T«a x 


« 0.55 
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which is used to find the machine active length to use in 
the circuit component calculations. The equation for rated 
torque can be manipulated to yield a quadratic expression 
for operating point slip, or a Newton's method convergence 
can be used to find operating point slip. Convergence on 
slip and active length through Newton's method is used to 
generate the machines of this thesis. 

Finally, rotor copper losses are (1 - slip)P; stator 
copper losses are found using Ii and Rl. 
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Table 42. Listing of induction machine design program 


•include "stdio.h* 

•include 'def.h' 

It prograa naee: ind.c for induction aachinee, 4/19/B7 tl 

long int seed; It start point lor randoi nuaber generator 1/ 

double b[26][JJ], hUlH331, xsC421; 

It b is 'best' array, h is 'hold* array, xs art winding (actors tl 

eainl) 

{ 

double design joint!), rnd.ealkO, sell), ke, kv, ainpwr, 
stepeire, randoal), absl), (req, r pa, ks£433; 
int p, iteration, i, j, best, print out!), loops, Hag; 

FILE tfopenf), lip; 

print!I'vnReading input data iro* IND.DAT . . 

Fp=fopen('ind.dat*,'r*); It input seed lor rando* nuebsrs tl 

(scan! lip,‘id’,listed); 

iscanilip, 'Id', Ip); It input nuabar oi pole pairs I/ 

iscarllip, 'Ill', (ainpwr); It input aachine power, derived ie ASSET 1/ 

ainp»rt»746.0; It convert to watts 1/ 

iscani lip, 'Hi', tike); It CERs (or Eiiective Height tl 

iscanilip, *X1(*, kkv); 

iscanilip, 'Ill', krpa); It aachine aax shait rpa tl 

idoselipi; 

print!I'VnHow aany loops do you want? '); 
scanil'Id', (loops); 

print!("NnSnDoing prograe calculations . . .\n'l; 

lor (i=l; i < 12; i+=2) It hareonic winding factors 1/ 

i 

ks,i]*swf(i); 

xslil ■ IkstiltksliD/lili); 

) 

(req=rpelp/40.0; It aax electrical frequency 1/ 

It MAIN BODY OF THE PRU6RAH tl 


ior 


li=l; i <= loops; +*i) 

{ 

print!("\nXd*, i); 
stepsize-O.1; 
iteration*0; 

(lag « 0; 

design pointleinpwr, p, ke, kv, ireq, (Flag); 

It put stuff in the hold array tl 


if (Hag » II 


( 


lilOHlBl « 10000000.0; 
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best = Oj 
) 

uhile ((iteration <« 10) il (Hag !* l)i 
( 

rnd.ealklainptiri p, stepsize, ke, kv, freq); 

/t stagger around I/ 

best=0; It index to best Ed of the lot tl 

for lj=l; j<M0; ++j) 

( 

if (hljltlBl < 1.0) 
continue; 

if IhljlUB) < hlbest)£183) 

best=j; It find the beet eachine 1/ 

) 

if (absUluOHlBl - hlbestlllBlI/hlOHlB)) < 0.005) 

It seal! ieproveeent in Ed II 

( 

stepswa/=2.0j 

♦♦iteration; 

) 

else It transfers best to 0 position 1/ 

( 

for (j=l; i <= 32; *+j) 

hiOHjl - hibestliji; 

) 


for (j=t; j<= 32; ♦♦]! 

biilIj]*hibest)Cj); 


It keep the best eachine tl 


best^l; 

for ti=l; i <= loops; sail 
1 

if IblilllBl < 1.0) 
continue; 

if (bli3E183 < btbestHlS): 

best=i; It find and keep the best of the best »/ 

1 

■inpwr/=746.0; It turn back into hp 1/ 

print.out<best, p, einper, ke, kv, rpa); It output to disk file tl 

fpsfopenl'ind.dat',V); It output seed tl 

fprintf(fp.'Id', seed); 

fprintf(fp,*\nZd*, p); 

fprintf<<p, , \nXlf*, oinpur); 

fprintf(fp,’\nZlf*, ke); 

fprintf(fp,*\nXIf a , kv); 

fprintf(fp,“\nXlf", rpa); 

fdoselfp); 

) 

It END Of BAIN PROSRAD; ALL THAT FOLLOd ARE FUNCTIONS 1/ 
double design.pointlainpur, p, ke, kv, freq. flat 
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/• deteraines a randn* design point tl 


double einper, ke, lev, freq; 
int p, Ifla; 

{ 

double r, Is, Ir, dcore, ds, dr, nr, g, *, ut, 1, i2rr, 
rl, r2, xbelt, xn, »ss, xse, xl, xrd, xrs, 
x«, vag, vt, via, rel, xel, taax, radical, sb, ila, ilb, 
x2, va, ph, pe, i2r, vol, ut, effey, ex, la, teinper, 
sa, sc, cu, siv, sev, riv, rev, loa, doa, r2siax, tt«ax; 
double ler, fl, f2, 13, slip, s*ax, jr, is, 
it, i2, xslot, utiron; 

extern double sqrtil, randoal), cosl), sinO, costii), sinhl), vratll, 
floor!), besipl), beskpO, tand, find.sivi), sttfOj 
int cslip=0, nr=71, ccr=t, ccg, ccs=l, k; 

»=2tPIIfreq; It synchronous frequency in rad per sec tl 

w*=w/p; It aech angular velocity in rad per sec tl 

1 * 10.0; It initial guess on length tl 

teax = ainpur/luelPSIl; It pull-out torque tl 

while (cslip 1= 1) 
t 

r=randnel)IRI1AX; 

if !(2IPIlrlfreq/pl < Mfl)t_TIP_SP£ED) It check rotor tip speed tl 
break; 

) 

lr=randoa(U0.5 ♦ 0.25; It rotor slot factor tl 

dr=randoe()tr/3.0; It slot no deeper than 331 of rotor radius 1/ 

dcore=(BRtr)/(BSATtp); It back iron depth tl 

ds*randoa(110.9tdcore; It slot depth < 901 of body depth to start I/ 

g=randoe(H(0.1lr - GRIN) ♦ GRIN; 

if (g < GRIN) 

g - 6RIN; 

ls=randoe()10.5 + 0.25; It stator slot factor 1/ 

f 1=0.; 
f2=0.; 
f3=0.; 

for (k*l; k < 100; +tk) 

1 

fl = (l./(p ♦ klnrll 1 (l./(p ♦ klnrl); 
f2 = U./(p - klnrl) I U./(p - klnrl); 
f3 ♦* fl ♦ f2; 

1 

while (ccr <- 5) It start convergence loop on r tl 

( 

if (ccs > 5) It ue’ve had trouble Mith jr tl 

{ 

r I* jr/JSHAX; 

if (r > (RAJ.TIP SPEEDIp/ul) 
r * RAX TIP SPEEDIp/e; 
if Ir ) RNfiX) 

r - RHAX; 

Ir l» jr/JSRAX; 
if (lr > 0.73) 

Ir * 0.75; 
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> 

»r = (BtPHIrtlr dr/2.>)/(**or ♦ Plllr); It rotor slot width 1/ 

Ur = 2IPIt(r - wr/4.0 - dr/2.0); It end riny length 1/ 

ccs * 1; 

while (ccs <* 5) It start convergence loop on ds and dr 1/ 

( 

ccg = 1| 

while (ccg <° 10) li start convergence loop on length il 

{ 

rl * (htCRHOt'l ♦ 2.3094JP1Kr+g)tCP/p))/(SSfIdsIPl*(r*g)*1s)i 
r2 = (hlxstUICRHOUKl ♦ (2lrlrllr)/!lerlptp))l/(PItrtlrtdr); 
xhelt = (4lwlllHUIrl(xs£51 ♦ xs[71)/(Pltqtplp))| 
tu = litwtltHUtrt(xsllll ♦ xstl31 ♦ xsU7] ♦ xstlOl ♦ xst231 ♦ 
xsl291 ♦ xsl311 v xsf371 4 xsUIDl/IPIlgtptp); 
xss = l3IHUIwtlldst!l4CP))/(lst2tPIIIr4g)); 
le = (PItrlCPttan(PI/3))/(3tp); 

xse - (27(wtlIHUtrtrt<-bcskp(p, ptPltr/leltbesip(p, ptP!tr/le)))/ 
(Pltlelplp); 
xl = xbelt 4 x 2 z 4 xss 4 xse; 
xrd = lhtRUtwtllxs[lllrt<3)/(PKg); 

xrs = (6»xstn»w«UHUMwr>0.5»dr))/(lrlPI*lr - (dr/2.0) - twr/B.O))) 
x2 = xrd 4 xrs; 

xe = (AtHUtwIltrlxstlD/IPlIglpIpl; 

vag = 2lrtBRIwlUswf ll)/p; It air gap voltage tl 

rel = (xitxetrll/trltrl 4 (xl4xe)t(xl4xa)l; 

It thev equiv resistance 1/ 

xel 3 (xattrllrl 4 xilxl 4 xltxa))/(rltrl 4 (xl4ie)t(xl4xe()| 

It thev equiv inductance tl 

r2seax = sqrtlrellrel 4 (xel4x2)l(xel4x2)l; It r2 at saax tl 
vt 3 vao/vrat(r2seax, x2, xl, xe, rl); 

It tereinal voltage at stax tl 
via ^ (vtlxa)/sqrt((xl 4 xelttxl 4 xe) 4 (rltrtl); 

It thevenin equivalent voltage 1/ 
ttaax = (1.3lvlatvla)/(wet(rel4sqrt(rettr»l4|xel4x2)t(xel4x2)))); 

It test eaxieue torque tl 
if (abs((teax-ttiax)/taax) <= 0.005) 

break; It we have convergence tl 

1 I- teax/tteax; It reset 1, reiterate tl 

♦♦ccg; 

if (ccg > 10) 

( 

If la * 1; 

printfl‘\n flag set on length"); 
return; 

) 

) It end convergence loop on length tl 

saax = r2/sqrt(reltrel ♦ (xel4x2)l(xel4x2l); 

slip 3 seax/3.0; It starting point for slip converge l/ # 

cslip 3 1| 

while (cslip < 3 20) 

{ It start convergence loop on slip I/ 
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t«inp*r * <3*vla*vlatr2/slip‘/< <rel+r2/*Up>*(rel+r2/slip> + 
(*elt*2>l(;el+*2))) 

if <abs((*inp*r-teinptir)/ainpttr) <= 0.0051 

break) It m have convergence tl 

slip 1= teinpur/einpurj 
++cslip; 

) It end convergence loop on slip I/ 

i2 = vaqtsqrt<(r2/slip)i(r2/slip) ♦ »2»x2>/«r2/slip)Mr2/slip) ♦ *21*21) 
jr = (3ti2)/(PltlridrlrtRSF)) It rotor current density tl 

if (jr > JSNflK) 

( 

dr U jr/JSHftJ) 
if (dr > (r/JI) 

dr = r/3j It reset to lieit tl 

) 

ila = Ir2/slip)»(xea*2) - *2lr2/slipj 
ilb = *2t(*eax2) + (r2/slip)t(r2/slip); 

il = vagtsqrt(ilatita ♦ ilbtilb)/(*st(*2l*2 ♦ (r2/slip)tlr2/slip)l); 
js = (3.0*il)/(PI*(r+g>lds»lslSSF); 
it (js > JSMAX) 

ds 1= js/JSRAX; 

if ((js <* JSMAXk U (jr <* JSMAX)) 
break) 

♦♦CCS) 

) It end convergence loop on ds and dr tl 

if (ccs <= 5) 

break; 

♦♦ccrj 

if (ccr > 5) 
l 

Ifla = 1; 

printf(’\n flag set on r‘)j 
return) 

) 

) It end convergence loop on r tl 

It calculations 1/ 


sea* = r2/sqrt(rellret ♦ (*ela#2)IUel+x2))| 

siv * find.sivll, r, g, ds, dcore, Is); .'t stator iron voluee 1/ 
riv = llPltrilr - 2$drllrlj It rotor iron voluee tl 
scv = 2IPIt(r+g) Jdsllstd ♦ 2.30MIPII(rtg)tCP/p)j 


rev * srtdrNnrll ♦ 2tlerlj 

cm » (rev ♦ scvllDCU) 

loa * 1 ♦ 4t(r*g)| 

doa = 2l(r+g+ds+dcorel| 

vol = VOLALLt(loatPltdoatdoa/4); 

*t = MTALLt (cm * Dt(BRNSStriv a siv))| 
et.iron - Dtlriv + slvl) 
ph ■ 31.8A225tB£TAIfreqlBftllHCltst_iron/D; 
nachine 1/ 

pe = (104234.9INl)$B3ATiBMTlTMTllfreqlfreql*t 


It stator copper voluee tl 
It rotor copper voluee tl 
It total copper Might tl 
It length-over-all tl 
It over-all-diaeeter tl 
It eachine envelop* voluee 1/ 

II eachine Might in kg tl 
It iron Might on!,-* tl 

It hysteresis loss in eetts, uses iron Might of 
iron)/lRHO»SI) 
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/I eddy current lust in watts, uses iron Height of sachine 1/ 

12r=3.0tittiltr1; It stator copper loss in Matts 1/ 

i2rr = sliptsinpwr; It rotor copper loss in watts 1/ 

vt - sqrt((vag*ilfrl)l(vag+illrl) ♦ (lltxl)l(illil)); It rated voltage tl 

va - 3lvtlil; It VA rating at tereinals tl 

effcy=leinpwr)/(»inpwr ♦ ph ♦ pe ♦ i2r ♦ i2rr 1; ew=wt ♦ kel(l-effcy) r kvlvol; 

(ecti"e weight tl 


h[01M*il; 

h[0H41=dcore; 

htO]tm=vol; 

h(01U61 : va; 

h[0][21]=siaa; 

ht01t261 s rlj 

ht0K31]=loa; 

It this 


hl0]121=freq; 
h[03(7]*ds) 
h [OH 12 l=wt; 


htO][J]=»| 

ht01[81=dr! 

htOH131=ph| 


hiOH41=r; 

h[0H91=lsj 

h[0H14]«pe; 

htOlf193=1; 

hI03C24]=slip; 

hl01E291=x2; 


hC0K5J=gi 

ht0]tl01*lr| 

hl0)tl53=i2r; 

ht03(201=i2rrj 

hC01(251=vag; 

h£03I303=r2i 


hC03[173=edfcy; htOJClfn=e»: 
hC03C22J=t*a*; hC0H233=vt; 
hCO][273=xl; htOK201=xe; 
hC03(323=doa; 

section just changed all the variables in the 'hold' array 1/ 


It Ef- 


return; 

) 


double rnd_walk(einpwr, p, stepsize, ke, kv, Ireql 

It walks about desiqn_point 10 tins 1/ 
double stepsize, winpwr, ke, kv, freq; 
int p; 

( 

double r, Is, lr, dcore, ds, dr, wr, g, w, we, 1, i2rr, ila, ilb, 
rl, r2, xbelt, xzz, iss, use, si, xrd, xrs, r2s«ax, tteax, 
re, vag, vt, via, rel, xel, teax, radical, sb, temper, 
x2, va, ph, pe, i2r, vol, wt, effey, ew, le, cw, siv, sev, 
sa, sc, riv, rev, loa, doa, jr, js, 
ler, f 1, (2, (3, slip, seax, il, i2, xslot, wtjron; 
extern double sqrtl), randoel), cost), sinl), coshl), sinht), vratl), 
besipl), beskpll, tend, find.sivl), swf(); 
int d, i=t, ccg, ccr, ccs, nr*71, k, cslip; 
f1=0.; 

(2=0.i 

f3=0.; 

(or (k=l; k < 100; e+k) 

( 

(1 = (l./(p ♦ klnrl) t 11./(p ♦ fctnr)); 

(2 = (l./(p - klnrl) I <l./(p - klnrl); 

(3 += (1 ♦ (2; 

1 

while (i <=> 10) H ten steps around the design point 1/ 

{ 

It read in the walk around the design point tl 
w*2l?It(reqj 

r»h[01Hll(l ♦ stepsizellrandoel) - 0.5)); 

if Uwlr/p) > HfllJIP.SPEEDI • 

continue; It go to next 1-1oop if violeted tl 

g=ht0HS3l<t e stepsizellrandoiO - 0.3)); 
if tg < BHIN) 
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g*6ftlN; /1 reset to the liuit 1/ 

dcore=(BRIrI/(BSATfp); It lost efficient use of iron tl 

ds=h(0](71l(l a stepsizeKrandoif! - 0.5)1; 

dr=h(0)[81l(l + stepsi 2 elirandoi(! - 0.5)1; 

if (dr > r/3.0) 

dr = r/3.0; It reset to the liiit tl 

!s=h[0H91l(l ♦ stepsizeKrandoiO - 0.5)); 
if (Is > 0.75) 

ls*0.75; It reset to the liiit tl 

if (Is < 0.25) 
ls=0.25; 

1 r=h£03[ 10J*(1 ♦ stepsizeKrandoiO - 0.5)); 
if (lr > 0.75) 

lr=0.75j It reset to the liiit tl 

if (lr < 0.25) 
lr-0.25; 

It lech angular velocity ;n rad per sec tl 
1 = h[0H191; It length starting point tl 

tiax = ht0H221; It pull-out torque tl 

ccs = 1; 

ccr - 1; It first tiie thru, no ad j in r tl 

while (ccr (- 51 It start convergence loop on r tl 

{ 

if (ccs > 5) 

< 

r 1= jr/JSHAX; 

if (r ) (MX.TIP SPEEDIp/u)) 
r = lifts TIP SPEEOlp/tt; 
if (r > fiNflX) 

r = RHftX; 
lr K jr/JSHflSs 
if (lr > 0.75) 

lr = 0.75; 

) 

«r = (BiPIllrtlr - dr/2.))/(4lnr ♦ Pltlr); It rotor slot width tl 
ler * IIPIKr - nr/4.0 - d>72.GI; It end ring length 1/ 

ccs s 1; 

while (ccs <= 5) It start convergence loop on ds and drl/ 

( 

ccg = 1; 

while (ccg <= 10) It start convergence loop on length tl 

rl = (6ICRH0M1 a 2,3094«Pl»(r+g!!CP/p))/(SSFIds»PI»(r+g)lls); 
r2 -- (ilxstlHCRHOWU V (2trtrtlr)/(lertptp)))/!Pltrllrtdr); 

(belt = (itutllHUtrtlxi(5] ♦ *s(71)/(PIIglp»pl); 
xzz ■ (6»emHU4rllxsUU ♦ xs(l31 e xt(17) a xs[19! a xi(231 a 
xs(291 a xsl3ll a xs(371 a xi[41)))/IPIigtpip); 
xss = (3tMJIwtlldsl(laCP)l/(li42tPII(rag))| 
le = (PIIr4CPItan(PI/3)l/(3lp){ 

xse - (27lulllNUtrlrt(-b«slip(p, plPItr/lellbesipIp, pIPItr/le)))/ 
(Pllletplp); 
xl = xbelt ♦ xzz + xss a xse; 
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xrd = (4MlllttiH»xstl]trlf3)/(PHql; 

xrs = (6lxsM 1M1 IHUtlwr+0.Stdr))/UrtPIt(r - (dr/2.C> - («tr/B.O»>); 
x2 = xrd * xrsj 

XI = (btnutwt Mrtxs(ll)/(Pltgtptp); 

vag = 2lr$Bfl»x»UsHf(l)/p| It air gap voltage '> 

rel = (xetxitrD/Irltrl a lxl*x»)l(xl+xi)l; 

It thev equiv resistance 1/ 

xel = (xeMrltrl + xltxl ♦ xltx*))/(rUrl + (xl+xiltlxl+xa)); 

It thev equiv inductance 1/ 

r2s»ax - .qrUrclirel ♦ (xeln2)l<xel+x2)); It r2 at imx 1/ 
vt * vag/vrat(r:!s*ax, x2, xl, xa, rl); 

It terainal voltage at siax 1/ 
via - (vttxel/sqrtUxl + xaltlxl * xi) ♦ <rl»r 1 >)j 

It thevenin equivalent voltage iI 
ttiax = U.5»vlii»vla)/(Nil(rettsqrt(reUreldxel*x2)l(xel+x2))>); 

It test aaxiiui torque 1/ 
if (absdtiax-ttiaxl/tiax) <= 0.005) 

{ 

ccg = lj 

break; It ue have convergence tl 

) 

1 t- taax/ttiax; It reset 1, reiterate tl 

♦+ccg; 

> It end convergence loop on length tl 

if (ccs > 5) 
break; 

siax = r2/sqrt(rellrel ♦ Ixel+x2ll(xel+x2)l; 
slip = saax/3.0; 

It starting point for slip converge tl 

cslip x 1; 
ehile (cslip <* 20) 

{ It start convergence loop on slip tl 

Uinpttr •» (3lvlalvlalr2/slip)/((rel*r2/slip)l(rel+r2/slip) + 
(xel+x2)J(xel*x2)); 

if (absKiinper-tiinpurl/einpur) <= 0.005) 

break; It ee have convergence tl 

slip U tiinpw/iinpnrj 

♦♦cslip; /I no lore than 20 tries 1/ 

) It end convergence loop on slip tl 

il * vagtsqrt((r2/slip)rr2/slip) ♦ x2!x2)/((r2/slip)t(r2/ilip> ♦ x2lx2); 
jr * (3li2)/(PUlrtdrlrlllSF>; It rotor current density tl 

if (jr > JSHhS) 

{ 

dr U jr/JSIIAK; 
if (dr > <r/3>) 

dr * r/3; It reset to liiit tl 

) 

ila » (r2/slip)»(xi+s2) - x2fr2/slip; 
ilb = x2t(xivx2l + (r2/slip)«(r2/slipl; 

il > vagtsqrtlilatila ♦ ilktilfel/(xil(x2tx2 ♦ (r2/slip)»(r2/slip)))| • 

js • <3.0tU>/<PItlr*g)l<iftlitSSF)| 
if (js > JSflAI) 
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ds I- js/JSHAI; 

if tljJ <= JSHAX) U (jr <= JSHA*)) 
break; 

+HC5J 

) It end convergence loop on ds end dr >/ 

if (fees <= 3) !! (ccg > 101) 

break; 

++ccr; 

) It end convergence loop on r tl 

if ((ccg > 10) 1! (ccs > 5) !! (cslip > 20) !! (ccr > 5)1 
i It no design convergence tl 

for (d=ij d <=10; t*dl 

hldlUSl = 10000000.0; 
printfl'burp ■); 
break; 

1 

It calculations tl 

seax = r2/sqrt(rellrel ♦ (xeln2)»(xel+x2l); 

siv = find.iivd, r, g, ds, dcore, Is); It stator iron voluee tl 

riv = UPItrKr - 2ldrtlr); It rotor iron voluee tl 

sev = 2IPUtr*glldsilsl(l ♦ 2.3094IPIt(r+g)»CP/p); 

It stator copper voluee tl 
rev = nrldrllnrll ♦ 2tlerl; It rotor copper voluee 1/ 

cu = (rev t scvltOCU; It total copper Height tl 

loa = 1 ♦ M(r*g); It length-over-all tl 

doa = 2»(r*g+ds+dcorel; It over-all-diaeeter tl 

vol = VOLALLt(loatPltdoatdoa/4); It eachine envelope voluee 1/ 

nt = NTALLIIcn ♦ DKBRNBStriv * sivi); It eachine Height in kg t ! 

ut.iron = Dklriv ♦ siv); It iron Height only tl 

ph = 31.U622SteETAIfreql6RUKCltHt.iron/ll; It hysteresis loss in natts, uses iron Height of 
eachine tl 

pe * (106236.OtNUtBSATtBSATtTltTltfreqtfreqtHt.iron)/(RKDtD); 

It eddy current loss in natts, uses iron Height of eachine tl 
i2r = 3.OtiIti Itr1; It stator copper loss in natts tl 

i2rr = slipleinpnr; It rotor copper loss in natts tl 

vt = sqrt((vagriltrl)t(vag*iltrll ♦ (iltxl)t(iltxl)); It rated voltage tl 
va = 3tvttil; It VA rating at the tereinals tl 

effcy=(einpnrl/(einpnr + ph ♦ pe + i2r ♦ i2rrl; en*nt ♦ ketd-effcyl ♦ kvtvol; It Ef 

fective Height tl 

h[iltll»il{ htiH2]=freq; bliH3)»e| htiH41»r; hli]t5)=g; 

htiU6]=dcore; hliK71=ds; hUHS.'I'dr; hmt9)*!s; htiHtOl-lr? 

h[i)Ul)°vol; h(iU12]=Ht; htiUtSl'ph; hliHlOpe; h[iUI53=i2r; 

hli]U6]=va| hliU171-effeyg htiUKIUen; htiU19]=l; htiJt203-i2rr; 
h(iU2U>ieax| hliK22J*tea»; hlilI23J*vt; hliH241*slip; h[iU2S)=vag; 

hli)1261=rl) hliK271=xl; htiU281*xe; hliU29]**2; hlilt30)=r2; 

h(iH31]=loa; h[ilt32)=doa; 

It this section just changed all the variables in the ‘hold* array tl 

Hi; 

) 


return; 

) 
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printout (bert, p, linpiir, ke, kv, rp«) 
int best, p; 

double ainpur, ke, kv, rpa; 
t 

char outfiletU); 

FILE Ifpo, IfopenOj 
int i; 

printfl‘\nNhat is the naae o< the tile tthere you Hant the output? ') 
scanfCZs - , outfile); 

fpo-Fopen(outfile, V)j 


fprintf(<po,*Xd", plj 
fprint#(fpo,*\nXI+*, ainpur); 
fprintf (fpo,'\nXlf\ ke); 
fprintf(fpo,'\ntlf‘, kv); 
fprintf (fpo/tnllf*, rpa); 
for (i=lj i <= 32; ++i) 

tprintf(fpo,'\nJe",b[best)li)); 
tprintftfpo, *\n’)j 
fcloseIfpo); 

) 


double find.sivd, r, 9 , ds. dcore, Is) 
double 1 , r, 9 , ds, dcore, Is; 

1 

double one, tuo, three, four; 

one * (r* 9 *ds*dcore)*(r+g+eSs+dcore) •• (re-g)ttr-eg); 
tuo • 2tPI< (r+gl tdsds; 

three = (r*g+ds+dcore)l(r+g»ds*dcore> - (r«-g*ds)• (reg+ds); 

four = lKPItone - tuo) * PIHMr+gltthree; 
return(four); 

) 

double vrat(r2saax, »2, ill, xa, cl) It ratio of vag/vt il 
double r 2 saak, x?, xl, xa, r 1 ; 

{ 

Double a, b, c, d, vrati 
extern double sgrtll; 

a « r 2 seaxt(xefxl) * rl»(x 2 ne); 
b * rttr 2 saai - xatx 2 - xll(x 2 *xa)j 
c * atr 2 saax - blx 2 ; 
d « btr 2 saax t alx 2 ; 


vrat * (aatsqrt(etc ♦ dtd))/tata ♦ blb)| 
return(vrat);} 
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Table 43. Lissting o-f induction ef-f iciency program 


(include “stdio.h“ 
lindude 'def.n' 

It prograe naie: ieff.c to find efficiency of induction tachines tl 

It xorks uith a single tachine tl 

lain!) 

( 

FILE tfopenl), ifp; 

double r, Is, lr, dcore, ds, dr, q, *, 1 , i2rr, va, ph, pe, i2r, vol, 

Nt, effcy, e», siv, riv, ut,iron, find_siv(), peeff, pgeff, teff, 
pirpr, geff-0.98, pceff-0.99, gr, dhp, rps, tinpttr, ke, kv, freq, 
sgeff-1.0, tar, *t, ri, r2, x 1, ta, vag, vt, via, ret, xel, radical, 
x2, sa, sb, sc, slip, slipt, slip2, star, il, >2, tear; 
extern double sqrtllj 
int e=0, f, p, i; 
char infileC141; 

printf recalculates efficiency for a single eachine.Vn'); 

xhile le != 1) 

{ 

f = Oj 

printfI'Nhat is the naie of the input file? *1; 
scanfl'Zs", in-file); 
fp = fopenlinfile, *r">; 

fscanflfp, 'id', Ip); It input nueber of pole pairs </ 

fscanfifp, "Ilf*, (tamper); 

eloper 1= 7(6.0; It non in *atts 1/ 

fscanflfp, "XIf*, (ke); 

fscanflfp, *Zlf, (kv); 

fscanflfp, ‘Ilf*, (rpi); 

fscanflfp, nif, (ill; 

fscanflfp, 'Ilf', (freql; 

fscanflfp, 'Ilf', in); 

fscanflfp, “XIf*, (r); 

fscanflfp, "Ilf*, (g); 

fscanflfp, 'Xlf*, (dcore); 

fscanflfp, “XIf*, (ds); 

fscanflfp, ‘Xlf, (dr); 

fscanflfp, *Xlf“, (Is); 

fscanflfp, 'Xlf', (lr); 

fscanflfp, “Xlf*, (vol); 

fscanflfp, *XIP, («t); 

fscanflfp, “Xlf“, (ph); 

fscanflfp, 'Ilf, (pel; 

fscanflfp, “ZIP, (i2r); 

fscanflfp, 'ZIP, (va); 

fscanflfp, "Ilf", (effcy); 

fscanflfp, 'Zlf, (ex); 

fscanflfp, 'Zlf', (I); 

fscanflfp, 'Zlf, (i2rr); 

fscanflfp, ‘Zlf*, (stall; 
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!scan!l(p, *XH“, titaaxI; 
fscantlfp, •X1C, Ivt); 

!scan!(fp, *11!*, tslipl) 
fscan!(!p, "Xl(\ tvag); 
fscandfp, *XU’, trl); 
fscanf(fp, *11!*, 1*11; 
fscantlfp, *Xlf*, lx*); 
fscanf(fp, *XH*, 1*2); 
fscantlfp, ■Xlf, lr2l| 
fdoselfp); 

utiile (1 != II 
{ 

print!(*\r.Wut is the sustained speed tachine horsepower? 
scanfl'llf*, tdhp); 
dhp I* 746.0; 

print!fbhat is the sustained speed lachine rp»? *); 
scanfl'llf*, Iptrptl; 

freq = ptrpelp/60,0; It PH frequency t/ 

via = (vtl*e)/sqrt((»l ♦ xtlilxl ♦ «*l ♦ (rUri)); 

It thevinin equivalent voltage tl 

rel = lx*tx*lrl)/(rltrl + (xl***)dlxUx*>); 

It thev equiv resistance 1/ 

xel - (xtdrllrl r xllxt ♦ xltxi))/(rl»rl ♦ lxl*x*)llxl+xi)); 

It thev equiv inductance tl 

sa * reltrel + Ixel+x2)l(xel+x2); It pieces a! slip quadratic tl 
sb = 2irellr2 - I3»llvlalvlalr2/dhp); 
sc = r2!r?j 

radical = sblsb - Msatsc; 
i! (radical < 0.0) 

abortl‘\nSot a negative radical in the slip eqn.’); 
slipl * (-sb ♦ sqrt(radical))/(2tsal; 
slip2 « l-sb - sqrt(radical))/I2tsa>; 
i! ((slip! < 0.0) M (slip2 < 0.0)) 

abort(*\nT«o negative slips.*); 

It no* we xill use the ssallest positive slip tl 
else i( ((slipl < 0.01 U (shp2 > 0.0)) 
slip * slip2; 

else l! ((slipl > 0.01 ti <$lip2 < 0.0)i 
slip = slipl; 
else 

slip = (slipl > slip2) ? slip2 : slipl; 

It tslip»*in(shpl, slip2) tl 

i! (slip ) siax) 

abortC\n slip is lore than siax.M; 

siv * find.sivll, r, q, ds, dcore, Is); It stator iron volute tl 

riv = ltPlirllr - 2Jdrllr); It rotor iron volute tl 

*t. iron = Dllriv ♦ siv); It iron Might only 1/ 

ph : 31.86225*B£TAftreqtBfl1tHCl*wt_iron/D; 

It hysteresis loss in watts, uses iron Miqht o! tachine tl 
pe * 1106234.9tNUiBSPTI8MTtTUTltfreqlfreqtwt iron)/(RHOtD); 

It eddy current loss in watts, uses iron night o! eachine tl 
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i2 * sqrt((slipldhp)/(3lr2l<I-slip>)l; 
11 = i2 + vag/ie; 


It laid current 1/ 

1 1 i2 plus «agnetizing current 1/ 


t2r=3,0liMillrl; /j stitor copper loss in natts tl 

i2rr = slipidhp; it rotor copper toss in natts tl 

peeff * dhp/(dhp ♦ ph ♦ pe ♦ i2r ♦ i2rr); 
print!C*\n Susta.ned efficiency is Ilf", peeff); 

printf r\n\nHhit is the endurance speed neckine horseponer? *)| 
scentt*Zlf', idhpl; 
dhp »= 746.0; 

printf('Nhat is the endurince speed uchine rpi? *>; 
scenft'Zlf*, Vperpel; 

freq - pnrpetp/60.0; It PR frequency tl 

sb = 2*rel*r2 - (3tltvlatvNtr2/dhp); 
radical * sblsb - flsalsc; 
if (radical < 0.0) 

abort(*Sn6ot a negative radical in the slip eqn.‘); 
slipl = (-sb + sqrt(radical)1/(2tsa); 
slip2 - (-sb - sqrt (radical) )/ (2$sa); 
if ((slipl < 0.0) M (slip2 < 0.01) 

abortl'\nTno negative slips.'); 
else if ((slipl < 0.01 U <slip2 ) 0.0)1 
slip = slip2; 

else if ((slipl > 0.01 U <slip2 < 0.0)) 
slip * slipl; 
else 

slip = (slipl > slip2) ? slip2 : slipl; 
if (slip > snail 

abortC\n slip is lore than siai.'l; 
ph x 31.B622S*BETAtfrcqtBRl>HC1 tut iron/D; 
pe » (106236.9INlUBSATiBSATITUTUfreqlfreql*t_ironl/(RHt!>l»| 

12 x sqrt((sliptdhp)/(3lr2t(l-slip)l); It load current tl 

il ' i2 ♦ vaq/»»; It i2 plus Monetizing current tl 

i2rx3.0liltillrl; It stator copper loss in natts II 

i2rr = slipidhp; It rotor copper loss in natts tl 

pneff • dhp/ldhp ♦ ph ♦ pe ♦ i2r * i2rr); 

printfC\n Endurance efficiency is Ilf', pseff); 

printf('\nSaee eachine? 'I; 
scant('Id', If); 
if (« == 0) 

continue; 
else if (f == 21 
t 

e * l; 

break; 

1 

) It end of f-loop )/ 

1 It end of e-loop 1/ 

) It end of lain prograi tl 


double find siv'l, r, g, ds, dcore. Is) 
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double 1 , r, 9 , ds, dcore, Is; 

t 

double one, two, three, (our; 

one = Irrgtds+dcorelMrtg+dstdcore) - Ir+glllr+g); 
tuo = 2»RIt(r+gltdslls; 

three » lr«g+ds*dcore!t(r*g*ds*dcore) - lr+g*d»INr*g*de)| 

<our = UlPUone - too) ♦ RIHKrtgilthree; 
return(four); 

) 
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Appendix E. Weight and volume algorithms 

The ability to correctly characterize a new naval ship 
technology depends in part on the ability to calculate the 
weight and volume associated with that technology. Algo¬ 
rithms for this purpose were taken from a variety of 
sources, including the Advanced Surface Ship Evaluation Tool 
(ASSET) theory manuals. The ASSET algorithms are the result 
of data analysis for naval ships that have been constructed, 
as well as studies for other ship designs. 

The Ship Work Breakdown Structure (SWBS) categorizes 
all ship weights. The general categories are: 


W100 

Ship structures 

W200 

Propulsion plant 

W300 

Electric generation plant 

W400 

Command and control equipment 

W500 

Auxiliaries and distributed 
systems 

W600 

Outfit and furnishings 

W700 

Ship armament 

WFOO 

Ship Loads 

Within W200 are 

several sub-groups that pertain to 

electric propulsion. 

They are: 

W235 

Electric propulsion devices 

W235.1 

Propulsion motors 

W235.2 

Propulsion generators 

W235.3 

Transmission lines and 
propulsion cables 

W235.4 

Cooling systems 

W235.5 

Switchgear 

W241.1 

Locked-train-double-reduction 
reduction gears 

W242 

Propulsion clutches and 
couplings 

W243 

Shafting 

W244 

Propulsion shaft bearings 

W245 

Propeller weight 
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ASSET allows SWBS groups to be adjusted in weight, 
which allows technology sensitivity analyses, such as this 
thesis, to be performed. Only a few of the W235 sub-groups 
need to be calculated outside of ASSET and adjusted within 
ASSET. These are W235.1, W235.2, W235.4, W235.5, W241.1, 
and W243. The W235.1 and W235.2 weights and volumes were 
calculated as part of the machine design. The rest of the 
needed W235 weights were calculated in "wt.c", a copy of 
which follows. These algorithms are the U. S. Navy 
standard, and were verified against actual ships and com¬ 
ponents . 

Shafting and transmission line weights are dependent on 
motor and generator positions within a ship. Their weights 
and volumes were calculated from ASSET equations, using the 
layouts of the baseline and variant ships. 

Where no volume equation was found in the ASSET 
documentation, or where the result of such an equation was 
unrealistic, a literature search generally found enough ac¬ 
tual equipment to permit a relationship to be empirically 
determined. A linear scaling of such volumes provided 
adequate results. 



Table 44. Listing of off-line weight and volume program 


•include *stdio.h* 

•define K 150.75 It gear hardness factor tI 

It prograi naie: «t.c, to find the weights of SMBS groups 1/ 


lain!) 

< 

FILE (fp, If open I)j 

double xg, xi, zg, zi, pa, ng, na, gra, grg, np, dhp, gn, ds, q, ne, 
*235, *2353, *2354, *2355, (*241=0.0, «243, **298=0.0, xprop, 
wpc, vpc, ppc, wexc, vexc, wbrk, vork, pi, v235, 
v2353=0.0, v2354=0.0, v2335, v24l, v243{ 
extern double sqrtll, powOj 

printfCtnThis prograa calculates SUBS 200 weights for screen output.')) 
printf(*\n\nHhat is the LC6 of the propulsion generator(s>? *lj 
scanfCllf', bxglj 

printfl’Nhat is the LC6 of the propulsion ictorisl? 
scantCZlf■, exi); 

printfl’Nhat is the VCG of the propulsion geneiatorfsl? * 1; 
scant("XIf*, izglj 

printfl'Nhat is the VC6 of the propulsion intorlsl? 
scant(’Ilf*, bzi); 

printfl’Nhat is the nuiber of propulsion generator(s>? ’>) 
scant<*X1f*, Ing); 

printfl’Nhat is the nuiber of propulsion lotorls)? ’)j 
scanfCllf’, bml; 

printfl’Nhat is the nuiber of gas turbines aboard? *); 
scant ("Ilf’, tine); 

printfl’Nhat is the rated horsepower of each gas turbine? ’)j 
scant CZlf •, tpa); 

printfl’Nhat is the propeller rpu? *)| 
scant CZlf’, bnplj 

printfl’How such horsepower is delivered to each propeller? 
scanfCllf’, Idhplj 

printfl’Nhat is the LC6 of the propeller? *1; 
scanfCllf, txprap); 

printfl’Nhat is the gear ratio at the propulsion iotor? ’)) 
scanfCllf’, tgrel; 

printfl’Nhat is the gear ratio at the propulsion generator? *1; 
scanfCllf*, tgrgl| 

It N235.3 Transiission lines 1/ 

w2353 * 0.00000?t<*i-xgtzq-zf27l»(pi»746/30000>| It LT, enhanced •/ 
v2J53 = 0.065451 ixi-xg(zg-zi+27>| /• cubic feet, enhanced »/ 

It LT 1/ 

It swag il 


It N235.4 Cooling systNi tl 
**2354 = 0.26tpalng/2240; 
v2354 = lOO.Otng; 
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It N235.5 Snitch gear 4/ 

**2355 = 0.264<ng ♦ 24n»>; It enhanced ASSET, LI 4/ 

v2355 - 45.04*2353; It switchgear voluie, ft A 3, enhanced tl 


ppc = 0.000939523tpa/ne; 
npc = n«t3.55Ippc/12; 
vpc * ni4340ippc/l2; 


It poner converter rating, HU, enhanced ASSET tl 
It xeight per tonv, LT, enhanced ASSET tl 
It vol per conv, LT, enhanced ASSET tl 


p« = pa4746/1000000; It eotor rating in IW tl 

eexc = (ne ♦ ni)4po*((p*/30), 0.3)j It weight of exciters, LT, enhanced tl 

vexc = 704«e*c; It vol oT exciters, tt A 3 tl 

etrk = 0.264pi4ni; It weight of braking resistors, LT, enhanced tl 

vbrk = 37lpeln*| It vol of resistors, fl A 3, enhanced 1/ 


w235 = *2353 ♦ w2354 ♦ *2353 ♦ *pc * nexc + wbrk ; 
v235 = v2353 + V2354 + v2355 * vpc ♦ vexc * vbrk} 

It 8241 locked-train double reduction gears tl 

if (grg !' 1.0) It there are pg gears, bt* gt and pg tl 

*241 = U.574paipo*!(grg*l), 3.0)4ng)/(36004grg4grg4grg4K); 
if (gr* !* 1.0) It there are pi gears, bt* pi and propeller tl 

*241 ♦» (l.S7ldhplpow((gri*l), 3.0)Ini)/(nplgntgri«K)| It LT 4/ 

v241 = *241134.412; it cubic feet, ratioed fi FF67 1/ 

It d243 Shafting tl 

ds = 2.i524po*((4.223474dhp/np), 0.333); /• shaft dial, sq inches 1/ 

*243 * 1.570Bldslds«2llxprop - xi - 1.0)/2240; It LT tl 

v243 = 0.010914dsids4(xprop - xi - 6.0); It ft A 3 1/ 

It N29B.1 LTDR operating fluids, additon to *t in ASSET 1/ 
if (grg 1= 1.0) 

*298 = 0.27lpatnq/2240; 
if (gri (» 1.0) 

*298 ♦= 0.27lpalm/2240; 

printf(*\n\nM235.3 Transiission lines 
printfI a \n8235.4 Coaling systeis 
printf f\nM235.5 Switchgear 
printfC'n Po*er converters 
printff'\n Exciters 
printfl*\n Braking resistors 
printfC\n\nN233 Electric propulsion 
printfI'\nh241 Reduction gears 
printf C\nN243 Shafting 
printf( a \nN298 Sear operating fluid 
1 


It LT tl 

I7.2H LT I7.21f ft A 3 a , *2333, v2353)| 

X7.2H LT I7.21f ft A 3 a , *2354, v2354)| 

17.214 Li X7.21f ft A 3 a , *2355, v2355); 

X7.21f LT 17.214 ft A 3 a , *pc, vpc); 

Z7.2H LT X7.21 f ft A 3 a , *exc, vexc); 

X7.21f LT I7.21f ft A 3 a , *brk, vbrk); 

17.2H LT X7.214 ft A 3, less P6s and Pits’, *235, v235); 
I7.2H LT X7.21f ft A 3 a , *241, v241); 

X7.21* LT X7.21f ft A 3 a , *243, v243); 

X7.21f LT\n\r, a , *298); 
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Appendix F. Advanced Surface Ship Evaluation Tool output 

The output of ASSET is in text and graphic form. The 
total text output for any particular synthesis run is more 
than thirty pages. Following are several graphic outputs of 
ASSET, showing the mechanical and electrical transmission 
ships used in this thesis. 




Figure 49. Hull isometric view of all ships 
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Figure 51. Plan view of subdivision in mechanical baseline 
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